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Abstract 

Variable-delay  tunable  optical  delay  line  or  optical  buffers  are  critical  for  the  devel¬ 
opment  of  all-optical  networks  components  as  well  as  interferometry  and  analytic 
instruments.  Recent  research  on  slow  light  may  hold  the  key  for  the  development 
of  the  Erst  practical  tunable  optical  delay  device.  In  this  research  an  atomic  vapor 
pulse  delay  model  is  developed  including  hyperhne  structure  and  Voigt  lineshape. 
Frequency  tunable  pulse  delays  of  0-37  ns  are  achieved  in  Cs  D2  at  various  vapor 
pressures  of  0.15-5.28  mTorr  between  78.9  °C  and  137.2  °C  in  agreement  with  model 
prediction.  Furthermore  full- frequency  Cs  D i  line  hyperhne  optical  delay  model  is 
validated  with  the  observation  of  delays  of  1.6  ns  to  24.1  ns.  Additional  optical  control 
of  delays  were  demonstrated  by  pumping  the  Cs  D2  transition  and  observing  resulting 
effects  in  the  Di  delay  spectrum.  For  a  pump  at  four  times  the  saturation  intensity, 
delays  are  reduced  by  a  maximum  of  78%  in  a  narrow  region  of  110  MHz  in  agreement 
with  a  Kramers- Kronig  model  prediction.  Diode-pumped  alkali  laser  (DPAL)  systems 
depend  on  accurate  models  to  scale  to  high  power  and  optimize  performance.  There 
is  currently  no  validated  bleach  wave  model  of  an  operating  DPAL  system.  This  work 
partially  examines  the  bleached  wave  temporal  dynamics  during  a  10  ns  high-power 
D2  pulse  event  with  0-400  Torr  helium  buffer  gas.  The  linear  dispersion  delay  model 
is  a  valuable  improvement  over  currently  used  Lorentzian  approximations  and  pro¬ 
vides  a  significant  addition  to  the  held  tunable  slow  light  delays.  This  work  presents 
the  hrst  reported  spectral  hole-burnt  linear  dispersion  delay  effects  in  an  alkali  vapor. 
The  hyperhne  relaxation  observations  present  insight  into  the  complex  bleach  wave 
dynamics  during  a  high- intensity  pulsed  pump  in  DPAL  systems. 
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TUNABLE  OPTICAL  DELAY  IN  DOPPLER-BROADENED  CESIUM  VAPOR 


I.  Introduction 


Variable-delay  tunable  optical  delay  line  or  optical  buffers  are  critical  for  the  de¬ 
velopment  of  all-optical  networks  and  optical  signal  processing  (Tucker  et  ah,  2005). 
Applications  for  tunable  optical  delay  systems  include  simple  optical  delay  lines,  opti¬ 
cal  switches  and  optical  storage  devices  for  optical  communication  systems  (Camacho 
et  ah,  2007;  Tucker  et  ah,  2005),  as  well  as  interferometry  (Shi  et  ah,  2007;  Purves 
et  al.,  2006)  for  remote  sensing  and  analytic  instruments  (Camacho  et  ah,  2006b). 
Recent  research  in  an  emerging  held  of  optical  physics  known  as  slow  light  may  hold 
the  key  for  the  development  of  the  Erst  practical  tunable  optical  delay  applications. 

Simply  stated,  “group  velocity  gives  the  velocity  with  which  a  pulse  of  light  propa¬ 
gates  through  a  material  system”  (Boyd  and  Gauthier,  2002).  Boyd  defines  the  terms 
“fast”  and  “slow”  light  by  comparing  the  group  velocity  to  the  velocity  of  light  c. 

Recent  research  has  established  that  it  is  possible  to  exercise  extraordinary 
control  of  the  velocity  of  propagation  of  light  pulses  through  a  material 
system.  Both  extremely  slow  propagation  (much  slower  than  the  velocity 
of  light  in  vacuum)  and  fast  propagation  (exceeding  the  velocity  of  light  in 
vacuum)  have  been  observed.  (Boyd  and  Gauthier,  2002) 

A  number  of  methods  have  been  employed  to  achieve  slow  light  including  exploit¬ 
ing  the  natural  resonant  absorption  in  atomic  systems  (Khurgin  and  Tucker,  2009), 
self-induced  transparency  (SIT)  (Allen  and  Eberly,  1987)  and  electronically-induced 
transparency  (EIT)  (Marangos,  1998;  Milonni,  2002;  Boyd  and  Gauthier,  2002;  Flcis- 
chhauer  et  al.,  2005;  Milonni,  2005).  Ideal  optical  system  components  for  optical  delay 
should  have  long  dclay/buffer  times,  no  loss  of  signal  amplitude  and  no  dispersion. 
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For  a  delay  system  to  be  of  practical  use  in  future  applications  it  must  be  capable  of 
storing  large  streams  of  optical  information,  and  allow  for  random  access  of  buffered 
optical  signals  all  in  compact  physical  dimensions  (Tucker  et  ah,  2005).  Fundamental 
understanding  of  the  delay  effects  on  optical  pulse  propagation  and  accurate  delay 
prediction  are  essential  for  any  slow  light  system. 

The  study  of  slow  light  in  atomic  vapors  have  played  an  important  role  in  the 
development  of  tunable  optical  delay  systems.  Many  of  the  preliminary  studies  in  slow 
light  were  accomplished  using  atomic  vapors  (Khurgin  and  Tucker,  2009).  Knowledge 
of  the  physics  of  slow  light  in  atomic  vapors  are  important  in  developing  any  practical 
optical  buffer  system.  Alkali-metal  elements  offer  well-documented  relatively  simple 
atomic  systems  with  physical  properties,  such  as  large  vapor  pressures  at  moderate 
temperatures,  that  make  them  ideal  for  slow  light  laboratory  media. 

It  has  been  theorized  (Brillouin,  1960;  Garrett  and  McCumber,  1970)  and  demon¬ 
strated  (Grischkowsky,  1973;  Tanaka  et  ah,  2003)  to  some  degree  that  optical  pulse 
delay  is  achievable  using  a  single  laser  tuned  in  the  neighborhood  of  resonant  atomic 
transitions.  Full  frequency-dependent  delay  models  have  not  been  developed.  In  this 
work  a  frequency-dependent  Kramers-Kronig  delay  model  is  developed  using  complete 
hyperhne  structure  treatment  with  Voigt  profile  lineshapes  in  a  Doppler-broadened 
alkali-metal  vapor.  Delay  predictions  are  compared  to  observed  pulse  delay  to  validate 
the  delay  model. 

In  EIT  narrow  transmission  windows  generate  areas  of  steep  dispersion.  Similarly 
it  has  been  theorized  that  spectral  holes  in  an  absorption  profile  could  be  used  to 
produce  tunable  delay  effects  (Shakhmuratov  et  al.,  2005;  Agarwal  and  Dey,  2003). 
The  slow  light  hole  burning  technique  has  been  demonstrated  in  solid  photonic  media 
(Bigelow  et  ah,  2003)  and  recently  rubidium  vapor  (Camacho  et  ah,  2006b).  Previous 
work  focused  on  frequency  modulated  pump  beam.  In  this  work  control  of  pulse 
delays  in  cesium  vapor  were  demonstrated  by  pumping  the  D2  (6  2Si/2  —  6  2P  i/2) 
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transition  and  observing  resulting  holes  in  the  Di  delay  spectrum.  Additionally  the 
frequency  dependence  of  the  delays  of  the  probe  laser  in  the  vicinity  of  the  spectral 
holes  is  examined  and  compared  to  the  Kramers- Kronig  model  prediction. 

Tunable  optical  delay  systems  used  in  this  research  are  similar  to  diode-pumped 
alkali  laser  (DPAL)  system  (Rabinowitz  et  al.,  1962;  Krupke  et  ah,  2003;  Ehrenreich 
et  ah,  2005;  Beach  et  ah,  2004)  and  application  of  the  knowledge  gained  in  pulse 
propagation  can  be  useful  in  the  understanding  of  the  temporal  dynamics  of  a  DPAL 
system.  Normally  three-level  lasers  are  not  possible  due  to  the  difficulty  in  achiev¬ 
ing  population  inversion.  The  alkali-laser  is  a  gas-phase  three-level  laser  using  the 
ground  state  (Si/2)  and  the  first  two  excited  states  (P1/2  and  P3/2)  of  the  alkali-metal 
vapor.  Population  from  the  upper  P3/2  state  is  transferred  to  the  lasing  P^  level 
through  a  collisional  gas  species  like  ethane  or  methane.  The  recycle  rate  at  which 
an  individual  photon  excites  the  atom,  collisionally  transferred  to  the  lasing  level 
and  lases  is  very  rapid,  the  quantum  efficiency  of  such  a  three-level  system  is  higher 
than  other  laser  systems,  and  the  gas-phase  laser  does  not  have  the  thermal  manage¬ 
ment  issues  of  solid-state  laser  systems.  While  DPAL  systems  have  been  in  use  for 
several  years,  complete  theoretical  description  of  the  system  is  still  being  developed 
in  order  to  optimize  the  lasers  for  scaling  to  high  output  power.  Optimization  of  a 
DPAL  system  depends  on  accurate  pressure-broadening  models  to  properly  match  the 
atomic  absorption  profile  to  the  broadband  diode  pump  source.  There  is  currently  no 
validated  bleach  wave  model  of  an  operating  DPAL  system.  Transverse  absorption 
spectroscopy  setup  is  used  in  this  work  to  partially  examine  the  temporal  dynamics 
of  the  bleached  wave  of  a  pulsed  alkali-metal  vapor  system. 
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II.  Background 


Theoretical  foundations  of  for  slow  light  have  been  well-known  for  for  over  a 
century.  The  concept  of  group  velocity  originated  in  the  study  of  sound  waves  by 
Sir  William  R.  Hamilton  (1839)  and  Lord  Rayleigh  (1894)  (Brillouin,  1960).  The 
problem  was  of  considerable  interest  in  connection  with  the  velocity  of  light  in  the 
early  1900’s  after  Albert  Einstein’s  Special  Theory  of  Relativity  because  of  the  possi¬ 
bility  of  propagation  faster  than  the  speed  of  light,  c  (Boyd  and  Gauthier,  2002).  In 
1907  Arnold  Sommerfeld  and  later  in  1914,  Sommerfeld  and  Leon  Brillouin,  developed 
the  theory  of  pulse  propagation  through  a  collection  of  Lorentz  oscillators  (treating 
photon-atoms  interactions  as  oscillators).  As  outlined  by  Brillouin  in  “Wave  Propaga¬ 
tion  and  Group  Velocity”  (Brillouin,  1960),  they  concluded  that  a  pulse  front  travels 
at  c  regardless  of  propagation  region  and  that  no  signal  velocity  could  exceed  c  in 
agreement  with  Einstein’s  theory  (Brillouin,  1960).  However,  group  velocity,  often 
confused  with  signal  velocity,  could  be  greater  than  c,  less  than  c,  or  even  take  on 
negative  values. 

Since  the  early  theoretical  work  atomic  vapors  have  played  a  role  in  the  study  of 
slow  light  and  continue  to  be  integral  to  the  fundamental  understanding  of  optical 
delay  research  (Khurgin  and  Tucker,  2009).  Slow  light  research  in  atomic  vapors  is 
broadly  categorized  by  the  method  used  to  achieve  tunable  group  delays. 

In  the  simplest  case  low  intensity  laser  light  near  naturally  occurring  resonant 
atomic  transitions  of  a  two-level  system  provide  the  necessary  dispersion  to  induce 
pulse  delays  (Khurgin  and  Tucker,  2009).  This  dissertation  will  refer  to  this  simple 
linear  case  of  absorption  or  refractive  index  variation  as  dispersion  delay.  Typical 
dispersion  delay  experiments  occur  within  the  linear  intensity  regime  and  can  be 
adequately  described  using  the  semi-classical  approach  where  the  atom  is  treated 
quantum  mechanically  and  the  field  is  treated  classically. 
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Historically  early  studies  of  slow  light  were  performed  as  non-linear  optics  shortly 
after  the  invention  of  the  laser  (Khurgin  and  Tucker,  2009).  Self-induced  transparency 
(SIT)  was  first  discovered  by  McCall  and  Hahn  in  1967  (McCall  and  Hahn,  1967). 
McCall  and  Hahn  passed  pulses  generated  by  a  ruby  laser  through  additional  ruby 
rods  and  observed  significant  group  delay  and  SIT  effect  where  coherent  light  near 
the  resonant  transition  frequency  is  transmitted  with  little  loss  of  pnlse  amplitude. 
Study  of  SIT  was  later  applied  to  gases  such  as  sulphur  hexafluoride  SF6  (Patel  and 
Slusher,  1967)  and  eventually  alkali-metal  vapors  (Bradley  et  ah,  1970;  Slusher  and 
Gibbs,  1972).  The  SIT  effect  is  dependent  on  the  nonlinearity  in  laser  intensity  but 
is  similar  to  dispersion  delay  in  the  wings  of  the  absorption  profile. 

The  third  major  category  of  slow  light  mechanisms  is  electrically-induced  trans¬ 
parency  (EIT)  discovered  by  Harris  et  al.  (Harris  et  ah,  1990).  EIT  is  a  coherent 
optical  effect  which  creates  a  narrow  transparency  window  due  to  quantum  inter¬ 
ference  between  states  to  eliminate  absorption  at  the  resonant  transition  frequency 
(Marangos,  1998;  Fleischhauer  et  ah,  2005). 

This  chapter  begins  with  the  necessary  foundation  in  semi-classical  theory  to 
explain  dispersion  delay  in  the  linear  region,  followed  by  a  discussion  of  SIT  in  terms 
of  non-linear  optics  and  EIT  in  common  quantum  optics  notation.  Pertinent  research 
relevant  to  the  understanding  of  slow  light  and  tunable  delays  will  be  presented. 
Finally  a  brief  description  of  diode  pumped  alkali  lasers  and  their  similarity  to  slow 
light  study  will  be  discussed. 
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Semi-Classical  Wave  Propagation 

From  Maxwell  equations  it  is  understood  that  electromagnetic  fields  have  wavelike 
solutions.  Consider  the  electric  field  in  the  common  monochromatic  plane  wave  form 

E(z,  t)  =  Eq  exp{i  (kz  —  Lot)}  (1) 

where  k  =  nuj/c  =  27t/A,  and  angular  frequency  u  =  2ttu  with  refractive  index  n  and 
the  speed  of  light  in  a  vacuum  c. 

Phase  velocity  v^is  the  speed  with  which  cophasal  surfaces  advance  in  the  medium 
(Born  and  Wolf,  1999).  The  phase  of  the  plane  wave  is  given  by 

(j)  —  kz  —  ut  (2) 

so  that  points  of  constant  phase  move  a  small  distance  Sz  in  a  small  time  St  such  that 

kSz  =  uSt  (3) 


The  phase  velocity  is  the  distance / time  or  v $  =  Sz/5t, 


c 

n 


(4) 


In  the  medium,  the  peak  of  a  pulse,  or  peak  of  the  group,  moves  with  the  same 


phase,  thus  the  phase 

.  ,  nui 

<p  =  kz  —  ut  =  — z  —  ut. 
c 


(5) 


6 


(8) 


where 


(9) 


Starting  the  works  of  Hamilton  and  Raleigh  in  the  held  of  sound  waves,  the  treat¬ 
ment  wave  propagation  was  extended  to  optics  and  the  speed  of  a  pulse  of  light.  Phase 
velocity  v $  and  group  velocity  vg  are  used  to  describe  the  velocity  of  the  phase  front 
and  propagation  a  packet.  Sommerfeld  and  Brillouin  defined  the  velocities  and  clearly 
outlines  how  subluminal  and  superluminal  velocities  do  not  violate  Einstein  causality 
(Brillouin,  1960).  Furthermore  Brillouin  discussed  the  Voigt  profile  and  the  effects  of 
dispersion  near  atomic  resonant  absorption  lines  (Brillouin,  1960).  Subsequent  works 
depend  on  these  fundamental  velocity  definitions. 

An  analytical  discussion  of  the  effect  on  pulse  propagation  near  an  absorption  line 
is  given  by  C.  G.  B.  Garrett  and  D.  E.  McCumber  (Garrett  and  McCumber,  1970). 
In  their  article  they  discuss  the  potential  for  superluminal  or  subluminal  velocities  in 
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various  regions  and  show  that  the  pulse  remains  unchanged  (substantially  Gaussian) 
(Garrett  and  McCumber,  1970). 

Complex  Index  of  Refraction 

Slow  and  fast  light  effects  are  dependent  on  rapid  variations  in  the  index  of  refrac¬ 
tion  as  a  function  of  frequency  (i.e.  dispersion)  that  occurs  in  the  vicinity  of  resonant 
material  properties  of  the  propagation  medium  (Boyd  and  Gauthier,  2002).  Index  of 
refraction  is  a  complex  quantity  given  by  the  equation 

n  =  nRJrinI  (10) 

Based  on  the  macroscopic  Maxwell  equations,  beam  intensity  is  determined  by 
the  Poynting  vector  S  (Griffiths,  1999) 

S  =  E  x  H 

=  — (E  x  B) 
ho 

=  e0c2(E  x  B) 

where  the  magnitude  of  E  and  B  are  related  by 

E  x  B  =  -E  •  E* 

c 

such  that,  the  time-averaged  Poynting  vector  in  terms  of  the  electric  held  yields 

/=  (S)  =e0c(E-E*>  (15) 
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(11) 

(12) 

(13) 

(14) 


Given  the  electric  field  as  in  Equation  (1),  using  complex  index  of  refraction  n 


so  that 


E(z,  t)  =  E0  exp  i(h—z  —  ut) 

L  c 


=  E0  exp 


.  00  .  UJ  ■ 

i  ( nR—z  —  oot)  —  rij—z 
c  c  . 


(16) 

(17) 


r. ,  co 

,  CO  1 

r  . ,  c:  .  a;  I 

Eq  exp 

%  ( nR—z  - 

-  L ot)  ~  Tlj—Z 

exp 

—%  ( nR—z  —  ut)  —  rij—z 

L  c 

c  J 

l  c  c  J 

Eg  exp 

—2  rij—z 

C  - 

El  exp [ 

—a(u)z] 

where 


ot{uj) 


2u 

=  — n, 
c 


or  a(y) 


47 TV 

c 


n,- 


Intensity  becomes 


/  =  gocEq  exp  \—a(uS)z] 


(18) 

(19) 

(20) 


(21) 

(22) 


which  decreases  exponentially  with  z  if  a  is  positive.  Hence  a,  in  cm-1,  is  called  the 
absorption  coefficient  if  a  is  positive  or  gain  if  a  is  negative.  Notation  may  vary  in 
other  sources  if  the  sign  is  included  in  the  coefficient. 

Given  the  intensity  at  the  beginning  of  the  medium,  z  =  0,  is 


I{z  =  0)  =  J0. 


(23) 


The  ratio  of  the  intensity  at  a  given  depth  z  relative  to  the  initial  intensity  Jq  is 


I(z)  =  /oexp  [—a{uS)z\ . 


(24) 
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This  relation  known  as  the  Beer-Lambert  law  (or  simply  Beer’s  law)  is  true  only  in 
the  linear  regime  when  the  intensity  /  is  below  the  saturation  intensity  Isat- 


Kramers-Kronig  Relations 

Both  parts  of  the  complex  index  of  refraction  gives  light  to  important  phenomena- 
the  real  index  of  refraction  effects  group  velocity  and  the  imaginary  index  of  refrac¬ 
tion  accounts  for  observed  absorption  in  a  dispersive  medium.  The  Kramers-Kronig 
relations  connect  the  real  and  imaginary  parts  of  the  complex  index  of  refraction, 
assuming  n(u)  is  regular  at  u  =  0  (Jackson,  1999)  (Hutchings  et  ah,  1992). 

-  1  =  -  P  r  d(  (25) 

7T  J  of-  U 


2  ,  2 co  ^  r  mo  - 1  „ 

ni  M  = - P  /  ~a. - T  (]t 

k  Jo  0  -  u1 


(26) 


where  P  is  the  Cauchy  principal  value.  If  either  term,  nR  or  n,.  is  known,  the  other 
may  be  determined.  From  Equation  (21) 


nM  = 


(27) 


thus  Equation  (25)  can  be  written  in  terms  of  the  absorption  coefficient  as  a  function 
of  linear  frequency  v 


nJv)  -  1  =  -  P 


7 r 


(28) 


2vr2 


<*(0 


(29) 


We  have  developed  an  expression  for  the  real  index  of  refraction  nR{y )  as  a  func- 
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total  power 
P 


cross  section  area 


thickness 


Figure  1.  Absorption  cross  section  construct.  Atoms  within  a  slab  of 
area  A  and  thickness  Sz  interact  with  light  that  strikes  in  cross  sectional 
area  crabs-  Total  power  P  is  evenly  distributed  across  the  area  A. 

tion  of  frequency,  known  as  dispersion  (Menzel,  2007).  Normally  index  of  refraction 
increases  with  increasing  frequency  known  as  normal  dispersion.  When  index  of  re¬ 
fraction  decreases  with  increasing  frequency,  it  is  known  as  anomalous  dispersion. 

In  spectroscopy  it  is  common  to  measure  attenuation  and  determine  the  index 
of  refraction  from  the  Kramers-Kronig  relations.  For  example,  an  observed  discrete 
absorption  line  spectrum,  discrete  step  width  AD  in  wavenumbers,  the  discrete  values 
of  the  dispersion  nR{jp  ■  AD)  at  the  frequencies  p  ■  AD  follow  from  the  Kramers-Kronig 
relation  (Menzel,  2007) 


(30) 


ra=0,7^p 


Absorption  Cross  Section 

A  phenomenological  construct  for  the  medium’s  macroscopic  absorption  is  de¬ 
scribed  using  the  following  model  (Siegman,  1986).  Consider  a  slab  of  material  in 
which  absorbing  particles  are  embedded  as  shown  in  Figure  1.  Number  density  N 
is  absorbers  per  unit  volume  N  =  n/V.  For  simplicity  let  the  particle  cross  sectional 


11 


area  crabs  perpendicular  to  the  incident  radiation  of  total  power  P  represents  an  area 
of  100%  probability  of  interaction.  Probability  of  interaction  is  0%  elsewhere  on  the 
surface  area  of  the  slab  A.  Intensity  is  defined  as  /  =  P/A.  The  relative  change  in 
intensity  of  the  incident  radiation  is  given  by 


out  I in 

1-in 

(31) 

A 

(32) 

-NVGabs 

A 

(33) 

NAcTabs^Z 

(34) 

A 

-Noabs5z. 

(35) 

Integrating  Equation  (35)  with  respect  to 
slab  z  yields  Beer’s  law  attenuation 


In  [I(z)]  -  ln[/0]  = 
I(z)  = 


the  distance  of  propagation  through  the 


1  N a absdt/ 

Jo 

(36) 

-NdabsZ 

(37) 

I0  exp  \—Noabsz\ . 

(38) 

Comparing  Equation  (24)  to  Equation  (38)  reveals  equations  of  similar  form,  such 
that 


a(u)z  =  N a  absz 

(39) 

ot{u)  =  Naabs. 

(40) 

Applying  quantum  mechanical  principles,  we  can  find  the  cross  section  for  absorption 
Pubs-  By  calculating  the  number  density  A  as  a  function  of  the  vapor  pressure,  we 
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can  calculate  the  absorption  coefficient  a(u). 

If  we  consider  cesium  vapor  an  ideal  gas,  we  can  find  the  number  density  N  =  n/V 
by  n/V  =  P/ksT,  where  P  is  the  cesium  vapor  pressure  as  a  function  of  temperature 
T  given  by  Equation  (104)  (Alcock  et  al.,  1984)  as  shown  in  Appendix  B  Figure  42.  At 
this  point  we  seek  crabs  as  a  function  of  photon-atom  interaction  to  End  the  absorption 
coefficient. 


Photon- Atom  Interaction 


Light-matter  interactions  are  described  by  quantum  mechanical  transition  prob¬ 
abilities.  Electric  dipole  interactions  are  the  dominant  factor  for  determining  atomic 
energy  level  therefore  the  matrix  elements  are  based  on  the  electric  dipole  moment 
(Fox,  2006).  We  seek  to  hnd  the  cross  section,  related  to  the  probability  of  interaction, 
in  terms  of  the  eflective  dipole  moment  Hf"f '  (Steck,  2009,  eqn  (43)) 

|  hf"F'  I"  =  ~^Sf"f'  |  (J  |  er  |  J')  |2  (41) 


where  Sf"F'  is  the  relative  strength  factor.  The  radiative  lifetime  rr  is  given  by 
(Loudon,  2000,  eqn  (2.3.21)) 


1  =  a21  = 

rr  37ie0hc3g2 


CO  n 


2J  +  1 


3iie0hc3  \  2  J'  +  1 


<  J I  er  |  J1) 


(42) 


where  e  is  the  unit  charge,  Uq  is  the  transition  frequency,  Di2  is  the  transition  dipole 
moment,  eo  is  vacuum  permittivity,  h  is  the  reduced  Plank  constant,  and  gi/g2  are  the 
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level  degeneracies.  Equation  42  may  be  used  to  solve  for  the  effective  dipole  moment, 


( J  |  er  |  J')  |  =  — 


A  ,2  1  3ne0hc3  (2 J'  +  1 


r 


2J+  1 


,2  c  1  ne0hc3  (2 J'  +  1 

„F,,F,|  -  sF"F7_s—  1 2jyr 


The  absorption  cross  section  is  given  (Bernath,  2005,  eqn  (1.57)) 


@abs 


2i r2  |  fipnpi 

3eo  he 


-v  g(u  -  u0) 


combining  Equations  (44)  and  (45)  yields 


O abs 


&abs 


27 r2  1  ne0hc3  f 2J'  +  1', 

c2  1 


r  ca 


-5 


247T  r 


F"F> 


\  2  J  +  1 
2J'  +  1\  i/ 


2  J  +  1  /  v. 


7^3  9{v  -  Vo) 


Line  Strength 


(43) 

(44) 


(45) 


(46) 

(47) 


The  remaining  pieces  in  our  goal  of  quantifying  the  group  velocity  are  related 
to  cesium  physical  properties  and  Dx  (62Si/2  — »  62Pi/2)  and  D2  (62Si/2  — >  62P3/2) 
D-line  properties  are  cited  from  published  literature,  notably  Steck’s  Cesium.  D-Line 
Data  (Steck,  2009),  Alcock  et  al.  vapor  pressure  equations  (Alcock  et  ah,  1984), 
and  Arimondo  et  al.  experimental  determination  of  alkali-metal  hyperhne  structure 
(Arimondo  et  ah,  1977).  See  Appendix  B  for  additional  cesium  information. 

The  S f"  F'  (6  —j)  symbol  provide  a  measure  of  the  relative  strength  factors  of  the 
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F"  — *  F'  transition,  given  by  (Steck,  2009,  eqn  (41)) 


S 


F"F' 


\  J 

(2F'  +  1)(2J+1)  { 


which  obeys  the  probability  sum  rule 


(48) 


Y,Sf»f>  —  1-  (49) 

F' 

The  values  are  summarized  in  Table  8  and  Table  9. 

For  the  S1/2  ground  state,  J  =  1/2,  and  total  nuclear  angular  momentum  7  =  7/2. 
F  =  J  + 1  so  that  the  magnitude  of  F  takes  on  the  values  of  |  J  —  I  |  ^F^J  +  7  ==>- 
F  =  3  or  4.  For  the  Di  excited  state  P1/2,  J  =  1/2  and  F  =  3  or  4.  See  Figure  44. 
For  the  D2  excited  state  P3/2,  J  =  3/2  and  F  =  2,3,4  or  5.  See  Figure  45.  The 
central  resonant  absorption  frequency  z/0  is  given  in  Table  5  and  Table  6. 


Lineshape  Functions 

Assuming  low  alkali  vapor  pressure  such  that  Cs-Cs  collisions  are  negligible  the 
pressure  broadening  and  spectral  line  shift  are  a  function  of  buffer  gas  pressure  alone. 
The  pressure  effects  are  tabulated  for  several  buffer  gases  in  B  Table  10. 

Pressure  effects  must  be  corrected  for  temperature.  Given  the  data  was  taken  at 
Tdata  =  294.15  K,  the  buffer  gas  pressure  effects  on  system  at  temperature  Tsys  (K) 
and  buffer  gas  pressure  Pj  (Torr)  by 
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where  A///\  is  the  Lorentzian  line  width,  is0  is  the  center  frequency  for  each  hyperhne 
line,  7*  is  the  pressure  broadening  coefficient  for  each  buffer  gas  species  and  Si  is  the 
frequency  shift  caused  by  each  buffer  gas  species. 

The  lineshape  function  g{y  —  u0)  is  a  convolution  of  homogeneous  pressure  broad¬ 
ening  and  inhomogeneous  Doppler  broadening  effects.  The  pressure  Broadening 
(. Lorentzian )  (Bernath,  2005,  eqn  (1.78))  is  given  by 

Az/L/(27r) 

(Az/l/2)2  +  (u-  z/0)2 


where  A is  the  Lorentzian  full  width  at  half  maximum  (FWHM). 


The  Doppler  Broadening  ( Gaussian )  (Bernath,  2005,  eqn  (1.93))  is  given  by 


9d{v  ~  v0)  = 


ln(2) 


Aun  V  71 


exp 


-4  ln(2) 


V—  Vo  ‘ 


A  v 


D 


(53) 


where 


A  vD  =  2z/0 


2kBT\n(2) 


me 


(54) 


The  Voigt  prohle  is  a  convolution  between  the  inhomogeneous  Gaussian  gj  and  ho¬ 
mogeneous  Lorentzian  gB  line  shape  functions  and  is  given  by  (Bernath,  2005,  eqn 
(1.83)) 

g{v  -  v„)  =  i  g,(^- v„)gH(v  -  i)  d(,  (55) 
A  useful  algorithm  for  the  Voigt  prohle  is  found  in  (Thompson,  1993) 


V(u,a )  =  Re  [exp[^2]erfc(^) 


(56) 
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where  z  —  a  +  iu,  a  —  y/In(2)^L,  and  u  =  2^/ln(2) 


,  .  /ln(2)  1  „ 

9v(v,  =  2\  - -r — Re 

7 r  Az// 


exp 


20^(2) 


z/  +  i 


erfc 


Auj 
-i  (  2y^n(2) 


z/  +  i 


Aj^g- 


Auj 


(57) 


Each  of  the  six  allowed  transitions  for  Cs  D2,  are  added  together,  weighted  by 
the  relative  ground  state  population  assuming  Boltzmann  population  distribution 
according  to 

_  grexp\-/3Er\ 

Pr  Ez  3*  exp  [-/?£*] 

resulting  in  the  total  absorption  cross  section. 

Any  effect  that  changes  the  absorption  profile  will  change  the  group  velocity  delay 
through  the  resonant  medium.  For  example,  changing  the  temperature  of  the  alkali 
cell  varies  the  cesium  concentration  in  the  cell.  Attenuation  depends  on  concentration, 
which  directly  effects  the  real  index  of  refraction.  Thus  varying  temperature  becomes 
a  knob  by  which  the  delay  is  tunable,  albeit  a  slow-tunable  method.  Other  possible 
methods  of  tuning  the  group  velocity  in  dispersion  delay  medium  include  buffer  gas 
pressure  or  frequency  tuning  of  the  pulse  across  the  absorption  spectrum.  Saturation 
effects  such  as  SIT  and  hole  burning  may  also  be  alternate  methods  of  optical  delay 
tuning. 


Wave  Propagation  and  Group  Velocity 

The  defining  theoretical  work  in  the  propagation  of  waves  comes  from  the  com¬ 
bined  works  of  A.  Sommcrfeld  and  L.  Brillouin.  The  definition  phase  velocity,  v^, 
and  group  velocity,  vg,  used  to  describe  the  velocity  of  the  phase  front  and  propa¬ 
gation  a  packet  as  given  by  Sommerfeld  from  works  originating  around  1910,  first 
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published  during  World  War  I,  and  reprinted  in  complete  form  in  Wave  Propaga¬ 
tion  and  Group  Velocity  (Brillonin,  1960)  in  1960.  Standard  definitions  for  group 
velocity,  signal  velocity,  and  velocity  of  energy  transfer  are  critically  importance  and 
often  misunderstood  (Born  and  Wolf,  1999;  Jackson,  1999;  Chu  and  Wong,  1982a, b). 
Brillonin  and  Sommerfeld  opened  the  possibilities  that  are  disturbing  or  surprising  to 
many,  such  as  subluminal  (slow  light),  superluminal  (fast  light)  and  negative  group 
velocities  (Brillonin,  1960).  Their  work  carefully  demonstrates  that  (first)  forerunners 
propagate  at  velocity  c  and  the  signal  arrives  at  the  aptly  named  signal  velocity.  Sig¬ 
nal  velocity  does  not  deviate  from  group  velocity  except  in  the  region  of  anomalous 
dispersion (dn/du  <  0)  (Brillonin,  1960,  p.78).  The  signal  velocity  is  always  equal  to 
or  less  than  c  and  Einstein  causality  is  not  violated. 

The  possibility  of  group  velocities  that  are  much  greater  than  c,  “or  even  negative” 
are  outlined  in  the  1970  work  of  Garrett  and  McCumber  (Garrett  and  McCumber, 

1970) .  There  is  a  comment  in  Brillouin,  citing  A.  Schuster’s  Einfuhrung  in  die  theo- 
retische  Optik  (Introduction  to  Theoretical  Optics)  [Leipzig,  1907],  in  which  he  states 
that  “group  velocity  has  meaning  only  so  long  as  it  agrees  with  the  signal  veloc¬ 
ity...  negative  parts  of  the  group  velocity  have  no  physical  meaning”.  Garrett  and 
McCumber  theoretically  demonstrate  through  analytical  approximation  and  numeri¬ 
cal  analysis,  that  there  remains  some  physical  significance  and  that  the  Gaussian  pulse 
is  not  significantly  distorted  if  the  pulse  is  “spectrally  narrow  in  comparison  with  the 
atomic  line“  (Garrett  and  McCumber,  1970).  An  additional  article  by  M.  D.  Crisp 
provides  theoretical  justification  in  the  photon-atom  interaction  as  the  cause  of  the  ob¬ 
served  group  velocity  variations  in  the  neighborhood  of  resonant  atomic  lines  (Crisp, 

1971) .  Crisp  states  that  the  shift  in  the  center  of  gravity  (maximum  pulse  amplitude) 
is  due  to  the  asymmetric  absorption  or  amplification  of  energy  from  the  light  pulse 
and  the  pulse  will  propagate  at  group  velocity  despite  significant  pulse  distortion  in 
a  dispersive  medium  (Chu  and  Wong,  1982a, b). 
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Another  important  theoretical  pulse  propagation  publication  was  printed  in  1993, 
wherein  Bolda,  Chiao  and  Garrision  prove  two  theorems  for  calculating  group  velocity 
in  dispersive  media  in  the  region  of  atomic  resonances  (Bolda  et  ah,  1993).  The  first 
theorem  states  that  regions  of  abnormal  group  velocity,  regions  where  vg  is  superlu¬ 
minal,  infinite  or  negative,  “are  required  by  the  causality  arguments  that  lead  to  the 
Kramers- Kronig  relations”  (Bolda  et  ah,  1993).  The  second  theorem  states  that  in¬ 
formation  about  the  abnormal  group  velocity  regions  is  obtained  from  the  absorption 
(or  gain)  curves  of  the  medium  (Bolda  et  ah,  1993). 

Theoretical  description  of  optical  pumping  mechanisms  and  development  of  the 
quantum  mechanical  formalism  started  with  the  work  of  Barrat  and  Cohen- Tanoudji 
(Barrat  and  Cohen-Tannoudji,  1961b, a),  as  well  as  Harper  and  Mathur  (Harper  and 
Mathur,  1967).  Pulse  propagation  effects  in  terms  of  the  density  matrix  formalism 
were  developed  by  Mathur,  Tang  and  Harper  (Mathur  et  al.,  1970).  These  works  are 
particularly  important  to  the  description  of  coherent  photon-atom  interactions  such 
as  SIT  and  EIT. 

Slow  Light  Mechanisms 

Previous  discussion  describe  pulse  delay  in  the  linear  Beer’s-law  regime.  Here 
pulse  delays  are  the  predicted  consequence  of  Maxwell’s  equations  and  the  Kramers- 
Kronig  relations.  There  are  several  other  related  methods  to  achieve  slow  light  (or 
fast  light)  such  as  saturation  delay,  self-induced  transparency  and  electromagnetically- 
induced  transparency.  It  is  possible  for  these  pulse  delay  mechanisms  to  be  observed 
simultaneously  in  certain  situations  and  determination  of  the  cause  of  the  observed 
delay  is  often  difficult.  This  section  describes  the  history  of  slow  light  research  and 
details  the  various  delay  mechanisms. 
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Linear  Dispersion 


Linear  dispersion  delay  occurs  at  pulse  intensities  well  below  the  saturation  inten¬ 
sity.  In  this  region  the  pulse  amplitude  transmitted  through  the  absorbing  medium 
follows  Beer’s  law  attenuation  and  the  observed  pulse  delay  is  not  intensity  dependent. 
Historically  initial  studies  in  pulse  propagation  appeared  in  the  context  of  nonlinear 
effects  primarily  self-induced  transparency,  discussed  in  a  later  section,  where  the 
observation  of  linear  and  nonlinear  effects  can  be  blended  within  a  single  medium 
(Khurgin  and  Tucker,  2009). 

Demonstration  of  optical  pulse  delay  is  a  direct  result  of  Kramers-Kronig  relations 
and  discussed  in  the  theoretical  predictions  of  Garrett  and  McCumber  (Garrett  and 
McCumber,  1970).  Initial  experiments  sought  to  demonstrate  superluminal  group 
velocity  in  a  gain  medium.  One  of  the  first  demonstrations  of  semi- classical  group 
velocity  and  the  Garret  McCumber  predictions  appear  as  a  measure  of  pulse  velocity 
greater  than  the  speed  of  light  in  vacuum  c  in  a  neon  absorption  cell  within  a  He-Ne 
laser  cavity  (Faxvog  et  ah,  1970).  Theoretical  discussion  of  these  results  is  reported 
in  a  series  of  works  by  Casperson  and  Yariv  (Casperson  and  Yariv,  1970,  1971). 

One  of  the  first  studies  in  linear  slow  light  pulse  propagation  appears  in  the  works 
of  Grischkowsky  in  1973  (Grischkowsky,  1973)  in  relation  to  earlier  reports  of  SIT 
(McCall  and  Hahn,  1967;  Faxvog  et  ah,  1970;  Slusher  and  Gibbs,  1972)  in  that  far  from 
resonance  pulses  can  exhibit  long  linear  delays  outside  the  SIT  region  (Khurgin  and 
Tucker,  2009).  Grischkowsky  reported  difference  in  delays  of  circularly  polarized  light 
in  the  Zeeman-split  rubidium  2P\/2  line  (Grischkowsky,  1973).  Short  pulses  of  narrow 
linewidth  low-intensity  near  the  resonant  line  in  the  dilute  alkali  vapor  were  produced 
group  velocities  as  low  as  l/14c.  Grischkowsky  contends  that  most  of  the  energy  in 
the  propagating  wave  was  contained  in  the  vapor  as  coherent  atomic  excitation  and 
quantitatively  explained  using  adiabatic  following,  in  which  “pseudomoments  of  the 
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atoms  remain  aligned  along  the  effective  field  of  the  laser”  (Grischkowsky,  1973). 

Linear  dispersion  pulse  delays  as  well  as  advancement  (superluminal  group  veloc¬ 
ity)  were  demonstrated  by  Chu  and  Wong  (Chu  and  Wong,  1982a, b)  in  an  epitaxially 
grown  GaP:N  crystal.  Pulse  propagation  in  dispersion  delay  gained  experimental 
notoriety  when  Tanaka  et  al.  (Tanaka  et  ah,  2003)  documented  the  observation  of 
negative  group  delay,  or  superluminal  velocity,  in  rubidium  vapor  Di  transition  using 
“conventional”  group  velocity. 

At  the  same  time  Agarwal  and  Dey  proposed  a  slow  light  method  using  linear 
dispersion  delay  in  a  Doppler-broadened  saturation  spectroscopy  system  burning  a 
spectral  hole  in  the  absorbing  medium  (Agarwal  and  Dey,  2003)  similar  to  the  EIT- 
related  work  in  a  ruby  crystal  by  Bigelow,  Lepeshkin  and  Boyd  (Bigelow  et  ah,  2003). 

The  next  series  of  significant  slow  light  research  came  from  the  Howell  group 
at  the  University  of  Rochester,  NY.  As  a  much  simpler  nearly  passive  alternative  to 
EIT,  using  the  atomic  doublet  structure  (rubidium  D-2  transition)  as  means  to  control 
velocity  of  light  through  the  medium  observed  in  2006  by  Camacho  et  al.  (Camacho 
et  ah,  2006a).  Camacho  outlines  several  desirable  criteria  for  slow  light  to  be  used 
in  communication  applications.  It  can  be  argued  that  dispersion  delay  has  several 
advantages  over  more  complex  EIT  systems.  One  such  potential  application  was 
realized  with  wide-band  tunable  optical  delay  using  slow  light  in  the  two-absorption 
atomic  resonance  setup  (cesium  vapor  D2  transition)  by  Camacho  (Camacho  et  ah, 
2007)  and  all-optical  delay  (Camacho  et  ah,  2007). 

Saturation  Delay 

For  pulse  intensities  above  the  linear  (Beer’s  law)  region  as  the  intensity  is  compa¬ 
rable  to  the  saturation  intensity  a  pulse  distortion  delay  effect  is  apparent.  Saturation 
delay  is  created  as  an  incoherent  delay  effect  which  occurs  when  the  leading  edge  of 
the  pulse  is  sufficiently  intense  such  that  half  of  the  atoms  in  a  two-level  absorber 
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are  excited  to  the  upper  level  and  the  absorption  line  becomes  saturated  (Allen  and 
Eberly,  1987).  The  rest  of  the  pulse  is  then  transmitted  completely  unattenuated  as  if 
the  absorbing  medium  were  transparent.  The  pulse  appears  to  undergo  delay  because 
the  transmitted  “center  of  mass”  of  the  pulse  is  moved  back  in  time  relative  to  the 
initial  pulse.  This  is  a  frequency  dependent  effect  such  that  saturation  intensity  is 
lowest  at  the  absorption  feature  resonance  peak.  Saturation  delay  is  necessarily  inten¬ 
sity  dependent  producing  pulse  distortion  and  delay  that  increases  to  a  maximum  and 
decreases  with  further  increase  in  intensity  (Allen  and  Eberly,  1987).  Pulses  capable 
of  bleaching  the  medium  are  transmitted  with  little  distortion  or  attenuation. 

Self-Induced  Transparency 

Self-induced  transparency  (SIT)  is  a  coherent  delay  effect  described  theoretically 
and  reported  experimentally  by  McCall  and  Hahn  (McCall  and  Hahn,  1967,  1969). 
The  SIT  delay  is  attributed  to  energy  leaving  the  pulse,  coherently  exciting  the  atoms 
of  the  absorbing  medium  and  energy  is  coherently  redeposited  into  the  trailing  edge 
of  the  pulse  with  no  loss  of  signal  amplitude.  Delay  is  produced  by  the  temporary 
storage  of  the  pulse  energy  in  the  absorber  (Allen  and  Eberly,  1987). 

McCall  and  Hahn  first  proposed  the  theory  of  self-induced  transparency  as  the 
coherent  coupling  of  a  coherent  pulse  and  dipole  moment  of  the  atomic  medium  and 
provided  experimental  results  in  a  series  of  cooled  ruby  lasers  (McCall  and  Hahn, 
1967).  Nonlinear  transmission  of  laser  light  is  observed  such  that  low-intensity  pulses 
are  highly  attenuated  and  high-intensity  light  are  transmitted  unattenuated,  which  is 
very  similar  to  the  delay  and  pulse  attenuation  results  observed  with  pure  incoherent 
saturation  delay.  McCall  and  Hahn  introduce  the  pulse  area  theorem  (McCall  and 
Hahn,  1967,  1969)  where  the  pulse  area  0  is  a  dimensionless  parameter  defined  by 
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(Fox,  2006) 


0 


Rl2 

h 


(59) 


For  pulses  of  area  equal  to  n  is  called  a  7r-pulse.  An  atomic  system  coherently  excited 
by  a  7r-pulse  will  have  the  atoms  in  the  excited  state  following  the  pulse,  but  returned 
to  the  ground  state  following  a  27r-pulse  (Fox,  2006;  Allen  and  Eberly,  1987).  The¬ 
oretically  for  pulse  areas  above  37t  the  atoms  are  coherently  excited  and  deexcited 
more  than  once.  The  resulting  pulse  breakup  is  observed  for  such  n7r-pulses  by  Mc¬ 
Call  and  Hahn  (McCall  and  Hahn,  1967)  and  is  presented  as  evidence  of  SIT  (Allen 
and  Eberly,  1987). 

Patel  and  Slusher  were  first  to  present  SIT  in  a  gaseous  medium  where  delays  of 
0.2  fis  were  observed  in  SF6  (Patel  and  Slusher,  1967).  They  also  noted  pulse  shap¬ 
ing  effects  where  the  pulse  become  more  symmetric,  the  pulse  tail  showed  increased 
intensity,  and  delay  was  intensity  dependent.  The  results  of  the  experiment  were  ques¬ 
tioned  regarding  the  degeneracy  of  SF6  by  Rhodes,  Szoke  and  Javan  (Rhodes  et  ah, 
1968).  Courtens  and  Szoke  demonstrated  that  in  the  limit  far  from  the  resonsant 
SIT  peak  the  pulse  delay  follows  linear  dispersion  theory  (Allen  and  Eberly,  1987; 
Courtens  and  Szoke,  1968).  Bradley,  Gale  and  Smith  reported  SIT  in  potassium  va¬ 
por  (Bradley  et  ah,  1970)  using  tunable  cavity  and  a  streak  camera  to  visualize  the 
resulting  pulse  delay.  In  1972  as  a  response  to  previous  degeneracy  criticism  Slusher 
and  Gibbs  reported  SIT  in  a  rubidium  vapor  made  nondegenerate  in  a  tuned  mag¬ 
netic  held  (Slusher  and  Gibbs,  1972)  where  nn  pulse  breakup  appeared  in  agreement 
with  McCall  and  Hahn  pulse  area  theory. 

Results  presented  in  the  present  work  do  not  intentionally  seek  pulse  length- 
intensity  combinations  which  produce  pulse  area  0  =  n2ii  in  which  SIT  is  observed. 
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Electromagnetically-Induced  Transparency 


The  final  slow  light  mechanism  to  be  discussed  is  electromagnetically-induced 
transparency  (EIT).  In  EIT  laser- induced  coherence  between  atomic  states  leads  to 
quantum  interference  between  excitation  pathways  producing  a  narrow  transmission 
window  at  an  absorption  resonance  (Fleischhauer  et  al.,  2005)  which  has  proven  as 
an  important  slow  light  (and  stored  light)  mechanism  dne  to  the  enhanced  nonlinear 
response  in  the  narrow  transparency  region  of  steep  dispersion.  Typical  EIT  schemes 
require  a  three-level  atomic  system  where  electronic  dipole  selection  rules  dictate  that 
two  pairs  of  the  levels  are  dipole  coupled  while  the  third  transition  is  dipole  forbidden 
(Marangos,  1998).  Review  of  EIT  systems  widely  available  in  the  works  of  Marangos 
(Marangos,  1998),  Fleischhauer,  Imamoglu  and  Marangos  (Fleischhauer  et  al.,  2005), 
Boyd  (Boyd,  2008;  Noginov  et  al.,  2009),  Boyd  and  Gauthier  (Boyd  and  Gauthier, 
2002),  Chiao  and  Milonni  (Chiao  and  Milonni,  2002)  as  well  as  Milonni  (Milonni, 
2002,  2005). 

Harris,  Field  and  Imamoglu  coined  the  term  EIT  in  1990  (Harris  et  al.,  1990) 
theoretically  predicting  the  coherent  EIT  effect.  The  following  year  Boiler,  Imamoglu 
and  Harris  demonstrated  EIT  in  a  Sr  heat  pipe.  In  1999  extreme  slow  group  velocities 
of  17 ms-1  were  achieved  using  EIT  in  a  Bose-Einstein  condensate  of  sodium  atoms 
(Han  et  al.,  1999;  Hau,  2001). 

The  next  milestone  experimental  measurement  of  dispersion  delay  comes  in  1995 
from  the  work  of  Xiao  et  al.  (Xiao  et  al.,  1995)  using  a  gas  of  hot  rubidium  atoms 
tuned  near  the  5S'i/2  — >  5P3/2  in  a  Mach-Zehnder  (MZ)  interferometer.  They  did  not 
directly  measure  the  group  velocity,  but  calculated  the  velocity  via  measurement  of 
the  phase  shift  (Boyd  and  Gauthier,  2002)  reaching  a  group  velocity  of  c/13.2  with 
the  coupling  laser  in  a  ladder-type  EIT  system.  Also  that  year  Kasapi  et  al.  reported 
A-type  EIT  group  velocity  of  c/165  ~  1.817  x  106ms-1  in  a  208Pb  vapor(Kasapi 
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et  al.,  1995). 

Alkali-metal  vapors  have  played  a  central  role  in  EIT  slow  light  research.  Zibrov  et 
al.  (Zibrov  et  ah,  1996)  demonstrated  enhanced  index  of  refraction  without  absorption 
in  Rb  D\ . 

Group  velocities  as  low  as  8 ms-1  inferred  in  rubidium  experiment  in  1999  by 
Budker  et  al.  (Budker  et  al.,  1999).  Hot  360  K  87Rb  vapors  produced  group  velocities 
of  90  ms^1  (Kash  et  ah,  1999)  using  a  A-type  EIT  scheme. 

Also  in  2003,  Kim  et  al.  (Kim  et  ah,  2003)  arbitrarily  changed  the  group  velocity 
from  subluminal  to  superluminal  by  changing  the  coupling  laser  intensity,  in  an  EIT- 
like  setup,  with  vg  from  c/3000  to  —c/14400.  At  this  point,  experiments  focused 
primarily  on  the  measurement  of  pulse  propagation  phenomena  in  the  vicinity  of  a 
single  resonant  absorption  line. 

In  a  2007  analytical  paper,  Shi  et  al.  (Shi  et  ah,  2007)  explored  theoretical  use  of 
dispersion  delay  systems  in  MZ  interferometer  system  with  slow  light  cell,  allowing 
for  tunable  optical  path  length  (OPL)  in  MZ  system.  Citing  recent  experimental 
results,  Shi  compared  possible  SL  interferometry  applications  (Shi  et  ah,  2007). 

It  has  been  demonstrated  by  Vudyasetu  (Vudyasetu  et  ah,  2008)  that  transverse 
images  can  be  slowed  and  even  stored  with  little  loss  of  coherence.  The  Howell 
group  again  displayed  tunable  delay  of  the  image  using  temperature  variation  in 
rubidium  vapor  dispersion  delay  cell  as  well  as  using  pump/probe  means  to  control 
the  absorption  and  group  delay. 

Following  the  University  of  Rochester  work,  analytic  study  of  dispersion  delay 
and  limitations  Shakhmuratov  and  Odeurs  (Shakhmuratov  and  Odeurs,  2008)  refer  to 
Camacho  (Camacho  et  ah,  2006a)  for  criteria  to  use  dispersion  delay  in  communication 
application.  Shakhmuratov  cites  low  group  velocity  without  (EIT)  driving  fields  for 
low-frequency  coherence. 

Just  as  linear  dispersion,  saturation  delay  and  SIT  can  appear  in  a  single  sys- 
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tem,  theoretical  description  of  a  combined  SIT  and  EIT  system.  Scully  et  al.  report 
theoretical  analysis  of  SIT  as  a  coupling  field  in  a  V-type  EIT  system  called  “mixed 
induced  transparency”  (MIT)  (Scully  et  ah,  2001). 

Electromagnetically-induced  transparency  is  typically  evident  in  systems  where 
the  Rabi  frequency  Qr  (Loudon,  2000) 

Mr  =  I  Id v2S0h  i  >  7  (60) 

is  much  greater  than  the  homogeneous  broadening  7  (Boyd  and  Gauthier,  2002).  In 
the  present  work  pump  intensities  are  such  that  EIT  was  not  observed. 

Scope  of  research 

This  work  seeks  to  develop  a  linear  dispersion  pulse  delay  prediction  from  a  full- 
frequency  hyperhne  structure  model  with  Doppler-broadened  Voigt  profile  lineshape 
in  cesium  and  rubidium.  This  pulse  delay  model  will  compared  to  observed  pulse 
delay  and  non- Voigt  double- Lorentzian  approximation  model  predictions. 

Once  a  reliable  Doppler-broadened  model  has  been  developed,  modification  of  the 
absorption  profile  will  be  attempted  using  a  second  laser  to  depopulate  the  ground 
state  and  observe  delays  in  a  narrow  spectral  hole  in  the  absorption  profile  in  a 
modified  sub-Doppler  spectroscopy  system. 

Finally  complete  theoretical  description  of  the  system  is  still  being  developed 
in  order  to  optimize  the  lasers  for  scaling  to  high  output  power.  Optimization  of  a 
DPAL  system  depends  on  accurate  pressure-broadening  models  to  properly  match  the 
atomic  absorption  profile  to  the  broadband  diode  pump  source.  There  is  currently  no 
validated  bleach  wave  model  of  an  operating  DPAL  system.  Transverse  absorption 
spectroscopy  setup  is  used  in  this  work  to  partially  examine  the  temporal  dynamics 
of  the  bleached  wave  of  a  pulsed  alkali-metal  vapor  system. 
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III.  Kramers-Kronig  Delay  Calculation 


This  chapter  originally  appeared  in  SPIE  Advances  in  Slow  and  Fast  Light  III 

(Anderson  and  Perram,  2010b). 


In  this  manuscript  predicted  delay  is  developed  for  a  pulse  propagating 
through  an  atomic  vapor.  Theoretical  description  for  pulse  delay  is  devel¬ 
oped  from  atomic  absorption  coefficient  with  hyperfine  structure  and  Voigt 
lineshape.  Predicted  delays  are  in  agreement  with  observed  delays  in  a 
cesium  vapor  cell  at  various  temperatures  between  78.9  °C  and  137.2  °C. 

Introduction 

A  number  of  applications  have  been  conceived  for  light  with  slow  group  velocities, 
including  optical  delay  lines,  interferometry  for  remote  sensing,  and  quantum  informa¬ 
tion.  Electromagnetically  induced  transparency  (EIT)  is  a  well  understood  method 
often  used  to  produce  slow  light  (Marangos,  1998;  Milonni,  2005;  Fleischhauer  et  al., 
2005).  Alternatively  optical  delays  in  the  vicinity  of  atomic  resonance  transitions 
have  also  been  demonstrated  (Grischkowsky,  1973;  Tanaka  et  ah,  2003;  Agarwal  and 
Dey,  2003;  Camacho  et  ah,  2006a).  One  of  the  early  slow  light  experiments  demon¬ 
strated  group  velocities  of  c/14  using  Zeeman  splitting  of  the  D\  line  of  rubidium 
(Grischkowsky,  1973).  Agarwal  and  Dey  used  an  absorption  spectroscopy  setup  to 
observe  delay  (Agarwal  and  Dey,  2003).  More  recently,  all-optical  delays  of  many 
times  the  pulse  width  have  been  demonstrated  in  the  double  Lorentzian  structure  of 
hot  alkali  vapors  (Camacho  et  ah,  2007)  where  delays  of  1.6  GHz-bandwidth  pulse  by 
25  pulse  widths,  or  7  ns,  have  also  been  observed  using  the  hyperfine  structure  of  the 
Cs  D2  line.  Delays  of  up  to  106  ns  have  been  achieved  using  the  hyperfine  splitting  of 
the  D-2  lines  in  85Rb  by  tuning  the  temperature  of  four  10  cm  cells  to  140  °C  (Cama¬ 
cho  et  ah,  2006a).  Indeed,  double  Lorentzian  systems  provide  large  bandwidth,  large 


27 


delays,  and  minimal  absorption  and  dispersion.  Rapid  tuning  of  the  delays  has  also 
been  demonstrated  by  varying  the  modulation  frequency  in  a  two-laser,  pump-probe 
configuration  (Camacho  et  al.,  2006b). 

Precise  absorption  models  which  fully  account  for  the  hyperhne  structure  can  be 
computationally  expensive.  Therefore  it  is  often  advantageous  to  approximate  the 
absorption  coefficient  or  complex  index  of  refraction  as  the  summation  of  Lorentzian 
lineshapes  at  a  pair  of  resonant  frequencies.  Lorentzian  absorption  lineshapes  pro¬ 
vide  adequate  predictions  for  cases  with  homogeneous  broadening  or  when  the  laser 
is  detuned  from  resonance  by  an  order  of  magnitude  greater  then  the  inhomogeneous 
linewidth  (Camacho  et  ah,  2007).  At  moderate  temperatures  (above  100 °C)  the 
Doppler  linewidth  for  the  cesium  D2  line  approaches  0.5  GHz  establishing  a  significant 
limit  on  applicability  such  that  Doppler-broadened  lineshapes  should  be  considered. 
Furthermore,  tuning  of  the  delays  by  varying  laser  frequency  is  largely  unstudied.  In 
the  present  work  we  develop  an  absorption  model,  calculate  predicted  pulse  delay  and 
compare  delay  measurements  to  predicted  delays.  Experimental  results  are  in  agree¬ 
ment  with  prediction  models  using  a  Voigt  lincshape  approximation  of  the  broadened 
absorption  with  full  D2  hyperhne  structure. 

Absorption  Coefficient 

Frequency  dependent  variations  in  the  index  of  refraction  near  resonant  absorp¬ 
tion  transitions  give  rise  to  dramatic  changes  in  group  velocity.  Index  of  refraction 
is  related  to  the  absorption  coefficient.  When  absorption  is  known  through  measure¬ 
ment  or  calculation,  delay  can  be  predicted.  A  rapid  accurate  numerical  absorption 
model  may  prove  useful  for  applications  needing  pulse  delay  predictions  over  a  wide 
frequency  range  in  the  neighborhood  of  resonant  atomic  transitions. 

The  absorption  is  a  function  of  frequency  v  as  well  as  temperature  T .  Frequency 


dependence  is  fundamental  to  photon-atom  interaction  in  that  probability  of  absorp¬ 
tion.  Temperature  dependence  arises  from  both  Doppler  broadening  and  temperature 
dependent  number  density  of  the  atomic  vapor. 

We  seek  the  absorption  coefficient  from  the  usual  definition  of  absorbance 

Ax  =  -  In (J/J0)  =  a(u,  T)l  =  atot(u,  T)n(T)l  (61) 

where  l  is  the  cell  length,  a  is  absorption  coefficient,  a  is  the  cross  section,  and  n 
is  the  number  density.  To  predict  the  pulse  delay  due  to  absorption,  the  absorption 
coefficient  a(y,T)  =  atot{v,  T)n{T)  must  be  determined. 

The  total  nuclear  angular  momentum  /  for  cesium  is  7/2.  For  each  hyperfine 
energy  level,  the  total  atomic  angular  momentum,  F.  is  given  by  F  =  J  +  I.  whose 
magnitude  takes  on  the  values,  \  J  —  I  \  ^  F  ^  |  J  +  I  \ .  Transitions  follow  the 
selection  rule  A F  =  ±1,  0.  In  this  paper  the  notation  used  for  the  D2  absorption 
transitions  follows  spectroscopy  notation  F"  — y  F',  where  F"  refers  a  hyperfine  level 
from  the  ground  state  S\ /2  and  F’  refers  to  a  hyperfine  level  in  the  P3/2  level. 

The  cross  section  a  for  a  single  hyperfine  transition  can  be  calculated  using 
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<7f"f'(v,T)  =  - — —  h2f„f,  v g{y  —  vf"F')  (62) 

o6q/1C 

where  n2F„FI  is  the  transition  dipole  moment  such  that 

/ tpiipr  =  SF"Fi  |  ( J "  |  |  er  |  |  T)  [  =  SF"F '  e2Dj„j,.  (63) 

The  lineshape  function  g{v  —  z/0)  is  discussed  later  in  this  paper. 

Strength  of  the  photon-atom  interaction  is  characterized  by  the  dipole  matrix 
elements.  Relative  strength  of  each  transition  can  be  expressed  as  a  reduced  matrix 
element  and  a  relative  strength  factor,  Wigner  6-j  symbol,  for  each  F"  — y  F'  transition 
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such  that 


2 


SF»F'  =  (2F'  +  1)(2J"  + 1) 


(64) 


The  relative  strength  factor  SpnF,  follows  the  summation  rule  ’YhF,  Sf"f’  =  1-  The 
frequency  for  each  hyperhne  F"  — y  F1  transition  is  calculated  (Arimondo  et  ah,  1977; 
Gerginov  et  al.,  2003).  Cesium  D2  properties  are  summarized  in  Table  1. 


Table  1.  Summary  of  hyperfine  transition  frequencies  and  relative  line 


strengths  cesium  D2. 


Transition 
F"  -a  F' 

^0  ^// 

(MHz) 

Strength  Factor 

S  F"  F' 

3  -A  2 

4831.14 

5/14 

3  — )•  3 

4982.37 

3/8 

3  — )•  4 

5183.65 

15/56 

4  — )•  3 

-4210.26 

7/72 

4  — )■  4 

-4008.98 

7/24 

4  — )•  5 

-3757.89 

11/18 

*  Calculated  from  references  (Arimondo  et  al.,  1977;  Gerginov  et  al.,  2003) 


In  solving  for  the  transition  dipole  moment  (London,  2000), 


1  9J"  n 2 

t  3 «0ftc3  gj,  rJ‘ 


(65) 


where  g%  the  J  degeneracy  such  that  gt  =  2  Jj  +  1,  so  that 


2  _  3e0/tc3  /2J'  +  l\  1 

/W'_  167^  \2J"  +  i)Sf"F't 


(66) 


Combining  Equation  (62)  and  Equation  (66)  the  cross  section  can  be  calculated 


(  rpN,  C2  V  f  2  J'  +  1  ^  c  1  { 

<?F"F'  {y  1  T)  —  - - 3  r>  Jn  i  1  SF"F'-g(v  -  ^o) 


8vr  vl  \  2  J"  +  1 


r 


(67) 


For  the  cesium  D2  transitions,  J"  =  1/2  and  J'  =  3/2  so  that  Equation  (67)  can  be 
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simplified  to 


OF"F'{y ,  T)  —  —  vq) 

r  47t  i/jj 


(68) 


The  relative  ground  state  population  must  also  be  considered  when  calculating  absorp¬ 
tion  for  each  hyperfine  F"  — y  F' .  Assuming  the  S\  /2  population  follows  a  Boltzmann 
distribution  as  a  function  of  temperature  and  the  ground  state  hyperfine  energies,  the 
relative  population  is  given  by 


If"  — 


gFtie~^Ep" 

Ei  SiFE> 


(69) 


where  f3  =  1/ksT  and  kF  is  the  Boltzmann  constant.  The  total  absorption  cross 
section  atot  is  the  sum  each  F"  — >  F'  cross  section. 


atot(u,T)  =  ^  If' 

F",F' 


q  1  C2  V 


(70) 


Lorentzian  lineshapes  are  often  used  to  approximate  atomic  transitions  for  Doppler- 
free  spectroscopic  techniques.  Close  to  the  absorption  peaks  or  when  inhomogeneous 
broadening  effects  are  considerable,  a  Voigt  profile  should  be  used  for  more  accurate 
depiction  of  the  atomic  absorption.  Thus  the  lineshape  itself  is  both  frequency  and 
temperature  dependent.  The  Voigt  profile  is  a  convolution  of  both  homogeneous  and 
inhomogeneous  broadening  lineshapes,  such  that 


gviy  -vo)=  /  gi{v 6  -  D) )gH(v  -  ^o)rfz/o 


(71) 


where  gi(F0  —  //0)  is  the  inhomogeneous  lineshape  and  g h ( v  —  f'0)  is  the  homogeneous 
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lineshape  given  by  the  formulas  (Bernath,  2005) 


9i{v  ~  vo) 

9h(v  ~  v o) 


2 

Aug 


Avh/2h 

(Avh/2)2  +  (u- 


-4  In  2 


v0)2 


f  V  -  Vo 

V  Apg 


2 


(72) 

(73) 


in  terms  of  the  Gaussian  (inhomogeneous)  full  width  at  half  maximum  (FWHM)  AvG 
and  the  Lorentzian  (homogeneous)  FWHM  AvL  . 

A  useful  algorithm  for  the  Voigt  lineshape  (Thompson,  1993)  is 


V(u,a)=  Re  exp(^2)erfc(^) 


(74) 


where 


z  =  a  +  iu 

(75) 

a  = 

Avg 

(76) 

«  =  2Vk2(i/A  "o). 

Avg 

(77) 

Three  broadening  mechanisms  were  included  in  this  model,  Doppler  broadening,  life¬ 
time  broadening  and  buffer  gas  effects.  Lifetime,  or  natural  broadening,  and  buffer 
gas  effects  are  homogeneous  broadening  effects.  Doppler  broadening  comes  from  the 
Maxwell-Boltzmann  velocity  distribution  as  a  function  of  the  temperature  T  and  the 
atomic  mass  M.  The  Doppler  FWHM  is  given  by  the  formula 


AvG  =  Vq 


8/csTln2 

Me2 


1/2 


(78) 


Homogeneous  broadening  from  natural  broadening  and  pressure  broadening  is  given 
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by  the  formula 


**  =  2^)+£ 


T. 


T, 


data 


1/2 


Pi 


(79) 


where  7,  is  the  pressure  broadening  constant  and  Pt  is  the  pressure  for  each  buffer 
gas  species.  The  pressure  broadening  constants  are  temperature  dependent  and  are 
scaled  by  the  temperature  of  the  system  Tsys  and  the  temperature  at  which  the  values 
were  measured  Tdata.  Buffer  gases  can  both  broaden  the  lineshape  profile  and  cause 
a  frequency  shift  as  a  function  of  buffer  gas  pressure.  Included  in  the  alkali  models 
are  helium,  neon,  nitrogen  and  ethane  (Rotondaro  and  Perrarn,  1997).  For  simplicity 
model  predictions  and  data  are  shown  without  pressure  broadened  lineshapes  so  that 
the  only  source  of  homogeneous  broadening  is  natural  broadening. 

Vapor  pressure  equations  for  cesium  used  in  this  model  are  given  by  the  “precise 
equations”  from  Alcock  et  al.  (Alcock  et  ah,  1984).  In  general  the  vapor  pressure  is 
an  exponential  function  Py(T).  Assuming  the  alkali  vapor  acts  as  an  ideal  gas  we 
can  determine  the  number  density 


PV  =  NkBT  ->■  n{T)  = 


PJTj 

kBT 


(80) 


The  model  calculated  absorption  coefficient  plots  are  shown  in  Figure  2  for  78.9  °C, 
102.5  °C  and  137.2  °C. 


Predicted  Delay 

Once  the  absorption  coefficient  is  known  the  predicted  delay  for  a  pulse  propagat¬ 
ing  through  the  alkali  medium  is  determined.  The  absorption  coefficient  is  related  to 
the  imaginary  part  of  the  complex  index  of  refraction  n  =  nR  +  in,,  by  the  equation 

"■ =  (81) 
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Figure  2.  Simulation  output  for  cesium  D2  absorption  coefficient  (cm-1) 
at  various  cell  temperatures  detuned  d=10  GHz  from  the  central  D2  tran¬ 
sition  frequency  351.7257185  THz. 


Using  the  familiar  Kramers-Kronig  relations 

(82) 

(83) 

we  can  connect  the  complex  index  of  refraction  or  the  absorption  coefficient  to  the 
real  index  of  refraction.  The  symbol  P  represents  the  Cauchy  principal  value  of  the 
integral.  Combining  Equation  81  and  Equation  82  yields  an  expression  for  the  real 
index  of  refraction  in  terms  of  the  absorption  coefficient 

(84) 

The  calculated  real  index  of  refraction  for  Cs  D-2  is  shown  in  Figure  3. 

The  group  velocity  shown  in  Figure  4  is  found  from  the  real  index  of  refraction 


Re  [er{y)\  =  1  H —  P 
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Figure  3.  Calculated  real  index  of  refraction  for  cesium  D2  at  vari¬ 
ous  cell  temperatures  detuned  ±10  GHz  from  the  central  D2  transition 
frequency. 


with  the  equation, 


Time  for  the  pulse  to  transit  the 
Figure  5  for  a  5  cm  alkali  cell. 

Delay  of  a  signal  through  the  cell 
signal  frequency  near  the  absorption 
from  the  peaks.  The  lead  FWHM  point  on  the  temporal  pulse  is  used  as  a  reference 
point  for  determining  delay.  Tunable  delays  are  achieved  by  selecting  various  signal 
frequencies  across  the  absorption  spectrum.  Near  the  absorption  peaks  pulses  realize 
the  greatest  group  delay. 


n.  + 


(85) 


alkali  cell  is  calculated  by  t  —  l/vg  and  shown  in 


is  determined  by  comparing  the  transit  time  of  a 
features  with  transit  time  of  a  signal  far  detuned 
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Figure  4.  Calculated  group  velocity  (units  of  c)  for  cesium  D2  at  vari¬ 
ous  cell  temperatures  detuned  ±10  GHz  from  the  central  D2  transition 
frequency. 


Figure  5.  Calculated  delay  (ns)  for  cesium  D2  at  various  cell  tempera¬ 
tures  detuned  ±10  GHz  from  the  central  D2  transition  frequency. 
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Figure  6.  Experimental  setup  for  delay  measurements. 


Experimental  Setup 

The  experimental  apparatus  for  delay  measurements  shown  in  Figure  6  is  similar 
to  the  absorption  spectroscopy  setup  of  Agarwal  and  Dey  (Agarwal  and  Dey,  2003). 

A  25  rnW  tunable  diode  laser  with  less  than  300  kHz  lincwidth  is  tuned  across  the 
Cs  D-2  transition.  The  1.5  mm  diameter  laser  beam  is  directed  through  an  optical 
isolator  and  a  pellicle  beam  splitter  to  a  fiber  coupled  Mach-Zehnder  interferometric 
lithium  niobate  (LiNbOs)  electro-optic  intensity  modulator  (EOM).  A  pulse  signal 
generator  drives  the  EOM  to  create  a  10  ns  FWHM  pulse  with  1  kHz  repetition  rate. 

Part  of  the  beam  is  sent  to  a  Bristol  wavemeter  with  0.003  cm-1  accuracy  for 
frequency  reference.  The  modulated  pulse  propagates  through  a  5  cm  long,  low- 
pressure  alkali  cell  enclosed  in  a  machined  aluminum  heater  block  maintained  at 
a  constant  temperature  between  78.9  °C  to  137.2  °C  using  a  temperature  controller 
system. 

The  pulse  delay  is  monitored  by  a  silicon  photodiode  detector  with  a  1  ns  response 
time.  At  a  fixed  alkali  temperature,  the  laser  source  is  tuned  across  the  cesium 
D2  spectrum.  At  each  sample  frequency  the  recorded  time-dependent  pulse  is  a 
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15  s  average.  This  was  repeated  for  approximately  90  frequencies  ±25  GHz  from  the 
central  D-2  frequency. 

Results 

Pulse  delay  data  were  recorded  in  the  range  from  T  =  78.9  °C  to  137.2  °C.  The 
delay  increases  and  transmitted  intensity  decreases  at  the  resonant  frequencies  and 
the  hyperhne  splitting  of  the  ground  2S'i/2(±//  =  3  and  4)  levels  is  clearly  evident. 
At  the  two  peaks  of  the  absorption  feature,  the  transmitted  pulse  intensities  were 
less  than  five  standard  deviations  (arbitrary  limit)  above  the  background,  making 
accurate  determination  of  the  delay  difficult.  Indeed,  the  optical  depth  at  line  center 
of  the  transitions  originating  from  F"  =  3  at  T  =  114  °C  is  £  =  l/er (z/0)n  =  0.02  mm. 
Delays  were  defined  by  the  leading  edge  at  half  the  peak  intensity. 

Pulse  delay  predicted  with  this  model  and  measured  pulse  delays  are  shown  in 
Figure  7.  Predicted  delays  closely  match  calculations  from  the  absorption  model 
described  previously  in  this  manuscript. 

The  delays  are  tuned  from  0  ns  to  30  ns,  up  to  3  pulse  widths  in  a  single  pass 
through  a  5  cm  cell,  for  the  T  =  114.0  °C  data.  The  present  results  compare  favorably 
with  the  previously  observed  delay  of  about  7  ns  at  zero  detuning  in  a  10  cm  cell  at 
120  °C  (Camacho  et  ah,  2007). 

Conclusion 

Tunable  slow  light  delays  are  demonstrated  through  frequency  modulation  across 
the  hyperhne  doublet  structure  of  an  alkali  metal  vapor.  Delays  of  0  ns  to  37  ns  have 
been  observed  in  a  single  pass  through  a  5  cm  cesium  cell  at  137.2  °C  by  detuning 
from  resonance  ±25  GHz.  Predicted  delays  closely  match  calculations  from  an  absorp¬ 
tion  profile  model  including  Voigt  lineshapes  and  full  hyperhne  structure.  Frequency 
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Figure  7.  Pulse  delay  measurements  in  5-cm  cesium  cell  at  78.9 °C, 
102.5 °C  and  137.2 °C  with  predicted  delay  calculated  from  simulation 
output. 


dependence  of  the  delays  are  sensitive  to  vapor  pressure  and  Doppler  broadening, 
constraining  vapor  temperature  uncertainties  to  about  0.1  °C. 
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IV.  Delay  Prediction 


This  chapter  originally  appeared  in  Physical  Review  A  81(3)033842  (Anderson  and 

Perram,  2010a). 


The  frequency  dependence  of  optical  delays  in  both  the  wings  and  core 
of  the  cesium,  D2  transition  have  been  observed  and  modeled  with  a  Voigt 
lineshape  convolved  with  the  six  hyperfine  components.  Tunable  delays  of 
0-37  ns  are  achieved  by  tuning  the  laser  frequency  through  resonance  at 
various  vapor  pressures  of  0.15-5.28  mTorr. 

Introduction 

A  number  of  applications  have  been  envisioned  for  light  with  slow  group  veloc¬ 
ities,  including  optical  delay  lines,  interferometry  for  remote  sensing,  and  quantum 
information.  Electro  magnetically  induced  transparency  (EIT)  is  often  used  to  pro¬ 
duce  slow  light  (Fleischhauer  et  al.,  2005).  However,  optical  delays  in  the  vicinity 
of  atomic  resonance  transitions  have  also  been  demonstrated  (Tanaka  et  al.,  2003; 
Grischkowsky,  1973;  Camacho  et  al.,  2006a).  One  of  the  early  slow  light  experiments 
demonstrated  group  velocities  of  c/14  using  Zeeman  splitting  of  the  D\  line  of  ru¬ 
bidium  (Grischkowsky,  1973).  More  recently,  all-optical  delays  of  many  times  the 
pulse  width  have  been  demonstrated  in  the  double  Lorentzian  structure  of  hot  alkali 
vapors  (Camacho  et  al.,  2007).  Delays  of  up  to  106  ns  have  been  achieved  using  the 
hyperfine  splitting  of  the  D2  lines  in  S5Rb  by  tuning  the  temperature  of  four  10-cm 
cells  to  140  °C.  Delays  of  1.6-GHz-bandwidth  pulse  by  25  pulse  widths,  or  7  ns,  have 
also  been  observed  using  the  hyperfine  structure  of  the  Cs  D-2  line.  Indeed,  the  double 
Lorentzian  systems  provide  large  bandwidth,  large  delays,  and  minimal  absorption 
and  dispersion.  A  more  rapid  tuning  of  the  delays  has  also  been  demonstrated  by 
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Figure  8.  Experimental  setup  for  delay  measurements. 


varying  the  modulation  frequency  in  a  two  laser,  pump-probe  configuration  (Camacho 
et  ah,  2006b). 

Previous  work  approximates  the  complex  index  of  refraction  as  the  summation 
of  Lorentzian  lineshapes  at  a  pair  of  resonant  frequencies.  Lorentzian  absorption 
lineshapes  provide  adequate  predictions  for  cases  with  homogeneous  broadening  or 
when  the  laser  is  detuned  from  resonance  by  an  order  of  magnitude  greater  then 
the  inhomogeneous  linewidth  (Camacho  et  al.,  2007).  At  moderate  temperatures  the 
Doppler  linewidth  for  the  cesium  D2  line  is  almost  0.5  GHz  establishing  a  significant 
limit  on  applicability.  Furthermore,  tuning  of  the  delays  by  varying  laser  frequency  is 
largely  unstudied.  In  the  present  work  we  report  on  delay  measurements  across  the 
resonances  and  show  the  results  are  in  agreement  with  Voigt  lineshapes  including  the 
full  hyperhne  structure. 

Experimental  Setup 

The  experimental  apparatus  for  delay  measurements  is  shown  in  Figure  6. 

A  20-mW  tunable  diode  laser  with  less  than  300  kHz  linewidth  is  tuned  across  the 
Cs  D2  transition.  The  1.5-nnn  diameter  laser  beam  is  directed  through  an  optical 
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1.0  - 


Figure  9.  Amplitude  modulated  10-ns  FWHM  laser  pulses  transmitted 
through  a  5-cm  Cs  cell  at  114.0 °C  (a)  detuned  by  24.9  GHz  and  (b) 
detuned  by  6.3  GHz.  The  FWHM  for  each  pulse  (solid  lines)  and  pre¬ 
trigger  noise  levels  (dashed  line)  are  illustrated. 


isolator  and  a  pellicle  beam  splitter  to  a  fiber  coupled  Mach-Zehnder  interferometric 
lithium  niobate  (LiNbOs)  electro-optic  intensity  modulator.  A  pulse  signal  generator 
drives  the  amplitude  modulator  to  create  a  10  ns  full-width  half  maximum  (FWHM) 
pulse  with  1kHz  repetition  rate  as  shown  in  Figure  9. 

Part  of  the  beam  is  sent  to  a  wavemeter  with  0.003  cm-1  accuracy  for  frequency 
reference.  The  modulated  pulse  propagates  through  a  5-cm-long,  low-pressure  al¬ 
kali  cell  enclosed  in  a  machined  aluminum  heater  block  maintained  at  a  constant 
temperature  between  78.9  °C  to  137.2  °C  using  a  temperature  controller  system. 

The  pulse  delay  is  monitored  by  a  silicon  photodiode  detector  with  a  1  ns  response 
time.  At  a  fixed  alkali  temperature,  the  laser  source  is  tuned  across  the  cesium  D2 
spectrum.  At  each  sample  frequency  the  recorded  time-dependent  pulse  is  a  15  s 
average.  This  was  repeated  for  approximately  90  frequencies  ±25GHz  from  the  central 
Do  frequency. 
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Figure  10.  Frequency  dependence  of  transmitted  pulses  delays  through 
114.0  °C  5-cm  cesium  cell. 

Results 

The  observed  transmitted  pulses  for  a  series  of  laser  wavelengths  at  T  =  114.0  °C 
are  illustrated  in  Figure  10.  Similar  data  were  recorded  at  T  =  78.9  °C  to  137.2  °C. 
The  delay  increases  and  transmitted  intensity  decreases  at  the  resonant  frequencies 
and  the  hyperfine  splitting  of  the  ground  2Si/2(F"  =  3  and  4)  levels  is  clearly  evident. 
At  the  two  peaks  of  the  absorption  feature,  the  transmitted  pulse  intensities  were  less 
than  five  standard  deviations  above  the  background,  making  accurate  determination 
of  the  delay  difficult.  Indeed,  the  optical  depth  at  line  center  of  the  transitions 
originating  from  F"  =  3  at  T  =  114  °C  is  £  =  l/a(v0)n  =  0.02  mm.  Delays  were 
defined  by  the  leading  edge  at  half  the  peak  intensity  and  are  illustrated  for  T  = 
78.9  °C  and  114.0  °C  in  Figure  11. 

The  delays  are  tuned  from  0  ns  to  30  ns,  up  to  3  pulse  widths  in  a  single  pass 
through  a  5  cm  cell,  for  the  T  =  114.0  °C  data.  The  present  results  compare  favorably 
with  the  previously  observed  delay  of  about  7  ns  at  zero  detuning  in  a  10-crn  cell  at 
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Figure  11.  Cesium  £>2  delay  data  for  78.9°C  (open  circles)  and  114. 0°C 
(open  triangles).  Predicted  delays  (solid  lines)  are  calculated  using  the 
full  spectrum  hyperfine  structure  absorption  model.  Maximum  error  in 
the  delay  measurement  is  ±0.15  ns  while  the  average  error  is  ±0.03  ns. 


120  °C  (Camacho  et  ah,  2007). 

The  pulse  delays  are  directly  related  to  the  absorption  spectrum.  Figure  12  il¬ 
lustrates  a  diode  laser  absorption  spectrum  of  the  Cs  D2  hyperfine  split  line  at  low 
pressure. 

The  absorption  is  readily  computed  from  the  hyperfine  splittings,  line  strengths, 
and  broadened  lineshape  (Pitz  and  Perram,  2008).  The  absorbance  A  at  a  given  cell 
temperature  T  and  length  l  is  related  to  the  effective  absorption  cross  section  o  and 
number  density  n  by 

71  =  -  ln(///0)  =  a(u,  T)l  =  atot(v,  T)n(T)l  (86) 

where  each  of  the  six  hyperfine  components,  F"  =  3  —  4  and  F'  =  2  —  5,  contributes 
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Figure  12.  Cesium  £>2  absorbance  for  a  2-cm  cesium  cell  at  25 °C  at  low 
pressure  (gray  line),  model  prediction  for  total  absorbance  (solid  black 
line),  and  model-simulated  absorbance  for  individual  hyperfine  transi¬ 
tions  (dashed  black  lines). 


to  the  cross  section  such  that 


&tot 


V  fF„S. 

F",F' 


1  C2 
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F"F' 


r  47tz/0f//f/3 


gv(y-v o  ) 


(87) 


The  frequencies  and  the  line  strengths  Sp"p'  of  each  hyperhne  component  are 

summarized  in  Table  2. 


Radiative  lifetime  for  the  Cs  D2  line  is  r  =  30.32(5)  ns  (Amini  and  Gould,  2003). 


The  statistical  distribution  of  the  populations  in  the  two  hyperhne  split  62S'i/2  levels 


is 


If" 


(2  F"  +  l)e"w 
k,  ab 

7  +  9e  kT 


(88) 


where  k  is  Boltzmann’s  constant  and  A E  =  9.192  631  770  GHz  is  the  energy  difference 
of  the  hyperhne  levels  of  the  Cs  625’i/2  ground  state. 

The  Voigt  prohle  gy(v  —  z/0)  is  a  convolution  of  both  homogeneous  and  inhomo- 
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Table  2.  Summary  of  cesium  D2  hyperfine  transition  frequencies  and 
relative  line  strengths. 


Transition 
F"  F' 

^0  pH  pr 

(MHz) 

Strength  Factor 

S  F"  F' 

3  — )■  2 

4831.14 

5/14 

3  — )■  3 

4982.37 

3/8 

3  — )■  4 

5183.65 

15/56 

4  — )■  3 

-4210.26 

7/72 

4  — )■  4 

-4008.98 

7/24 

4  — )■  5 

-3757.89 

11/18 

1  Calculated  from  references  (Arimondo  et  al.,  1977;  Gerginov  et  al.,  2003) 


geneous  broadening  lineshapes  and  is  approximated  by  (Thompson,  1993) 


gv(u,a )  =  Re[exp(^2)erfc(^) 


(89) 


where 


a  +  iu 

(90) 

A  vD 

(91) 

2Cta2  (\  "»)2 

Aud 

(92) 

The  inhomogeneous  Doppler  width,  Az/p  «  0.43  GHz  at  114  °C,  is  relatively  large 
compared  to  the  homogeneous  broadening.  The  cell  contains  no  buffer  gas  and  the 
alkali  pressure  is  low,  so  that  the  homogeneous  broadening  approaches  the  natural 
limit  of  A uh  =  5.2  MHz. 

The  temperature  dependence  of  the  cesium  vapor  pressures  P  are  derived  from 
Alcock  et  al.  (Alcock  et  ah,  1984)  and  the  corresponding  concentrations  n  =  P/kT 
range  from  4.08  x  1012  to  1.24  x  1014  atoms/cm3  for  T  =  78.9  °C  to  137.2  °C.  The 
modeled  absorption  from  Equations  (86)  and  (87)  compares  very  well  with  the  ob- 
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served  spectrum,  as  illustrated  in  Figure  12.  The  only  adjustable  parameter  employed 
is  a  small  correction  to  the  thermocouple  temperatures.  The  reported  temperatures 
are  determined  by  matching  the  observed  and  predicted  absorption  or  delays.  For 
example,  the  thermocouple  temperatures  were  taken  at  a  position  not  touching  the 
heater  element,  in  contact  with  only  the  alkali  cell  and  exposed  to  room  temperature 
air  within  the  laboratory.  For  data  reported  at  114  °C,  thermocouples  read  113.7  °C 
and  hyperfine  model  temperature  was  adjusted  until  the  observed  delay  best  fit  model 
predicted  values  within  0.1  °C.  The  hyperfine  splitting  of  the  excited  state  is  smaller 
than  the  Doppler  width  as  shown  in  Figure  12. 

The  delays  are  readily  predicted  from  the  absorption  coefficient  by  employing  the 
Kramers- Kronig  relationship  between  the  real  and  imaginary  portions  of  the  index  of 
refraction  n  =  nR  +  ini , 


nR  =  1  +  —  P 

Z7T 


£M0 


di 


where  P  represents  the  Cauchy  principal  value  of  the  integral, 
related  to  the  imaginary  part  of  the  index 


a(u)  =  - — unAu). 
c 


(93) 

The  absorption  is 


(94) 


Thus  the  real  portion  of  the  index  is  related  to  Equations  (86)  and  (87)  by 


nR  =  1  + 


<*(0 

e-»2 


dA 


(95) 


Finally,  the  observed  delays  are  predicted  as  td  =  l/vg ,  where  the  group  velocity  is 
defined  as 


v 


9 


C 

nR  +  v 


dnR  ' 
dis 


(96) 
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Figure  13.  Comparison  of  group  velocity  at  114 °C  computed  from  the 
hyperfine  model  of  Equation  (87)  (solid  line)  and  double  Lorentzian  index 
of  refraction  approximation  (dashed  line).  Percent  difference  between 
models  is  shown  for  comparison. 


These  predicted  delays  are  overlaid  with  the  experimental  observations  in  Figure  11. 
The  agreement  is  excellent  for  all  frequencies  and  temperatures. 

The  delays  are  often  computed  assuming  two  Lorentzian  profiles  split  by  the 
ground  state  hyper  hue  structure  (Camacho  et  ah,  2007,  2006a,  2007;  Shi  et  ah,  2007). 
This  approximation  is  adequate  when  the  laser  is  detuned  from  the  resonances  by  an 
order  of  magnitude  greater  than  the  inhomogeneous  linewidth  (Camacho  et  ah,  2007). 
However,  significant  differences  in  the  frequency  dependence  of  the  predicted  group 
velocities  are  observed  when  the  two  approaches  are  compared,  as  seen  in  Figure  13. 
For  the  two  Lorentzian  model  of  the  index  of  refraction  (Camacho  et  ah,  2007) 


n(S) 


9 1 


OJ  i  +  i'y 
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5  —  cu2  T  i'j 
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where  N  is  number  density,  \n\2  =  2.1947  x  10-29Cm  is  the  far  detuned  dipole 
moment  (Steck,  2009),  g\  =  7/16,  r/2  =  9/16,  uq  =  —3880MHz,  uj2  =  4980MHz, 
7  =  2.617 45  MHz  (27  is  the  Lorentzian  FWHM  «  natural  linewidth),  and  A  = 
2.984  65  MHz  at  114  °C. 

Across  the  vapor  pressure  equation’s  temperature  range  (Alcock  et  ah,  1984), 
300  K  to  550  K,  used  in  the  hyperhne  Voigt  profile  model,  the  group  velocity  deviation 
due  to  Doppler  broadening  remains  approximately  as  shown  in  Figure  13.  Indeed  the 
best  case  for  model  agreement  occurs  with  no  pressure  broadening.  Between  the 
absorption  peaks  at  zero  detuning  a  7.3%  difference  in  the  predicted  group  velocity 
is  observed  for  the  model  baseline  conditions  of  Figure  13.  For  comparison  at  higher 
pressure  broadening,  with  1/4  atmospheres  helium  buffer  gas  and  the  alkali  cell  at 
114  °C,  the  Lorentzian  FWHM  Avu  is  approximately  1/3  of  the  hyperhne  splitting. 
In  this  case,  where  the  two  absorption  peaks  remain  relatively  well  defined  as  the 
absorption  coefficient  at  zero  detuning  is  approximately  half  the  peak  values,  the 
two-Lorentzian  approximation  group  velocity  deviates  49.6%  from  the  full  hyperhne 
structure  Voigt  prohle  group  velocity. 

Pulse  Delay  Images 

Pulse  delays  were  observed  using  a  setup  similar  to  the  aforementioned  single- 
beam  pulse  delay  measurements  using  a  Princeton  Instruments  PI-MAX  intensified 
charge-coupled  device  (ICCD)  gated  camera  system,  shown  in  14.  A  beam  expander 
was  placed  before  the  alkali  cell  so  the  expanded  pulse  hooded  the  alkali  cell  and 
surrounding  space  inside  the  heater  block.  The  PI-MAX  was  used  in  place  of  the 
silicon  detector  without  the  compound  lens  attached  so  the  expanded  beam  shown 
directly  on  the  PI-MAX  ICCD  array. 

The  PI-MAX  camera  used  is  capable  of  imaging  frame  gates  as  small  as  0.5  ns. 
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Oscilloscope 


Figure  14.  Experimental  setup  for  pulse  delay  observation  with  Prince¬ 
ton  Instruments  PI-MAX  image  array. 

The  camera  used  contains  256  x  1024  spectrometer  ICCD  array.  The  camera  was 
turned  on  its  side  to  align  the  ICCD  array  vertically  with  the  expanded  beam  and 
the  vapor  cell  such  that  pulses  arriving  both  inside  and  outside  the  cell  could  be 
imaged  simultaneously  thus  providing  temporally  imaged  evolution  of  delayed  and 
non-delayed  pulses.  Close-up  detail  of  the  image  array  relative  to  the  alkali-metal 
vapor  cell  are  shown  in  15. 

Four  0.5-ns  frames  of  a  multi-frame  TIFF  imaged  with  the  PI-MAX  ICCD  are 
shown  in  16.  Each  image  of  the  256  x  1024  array  captures  laser  pulses  outside  the 
alkali-metal  vapor  cell  shown  at  the  top  of  each  image  and  laser  pulses  imaged  through 
the  reference  cell.  The  edge  of  the  Pyrex®reference  cell  is  highlighted  with  dashed 
curves.  The  series  of  images  details  the  time  evolution  of  a  10-ns  FWHM  pulse  delayed 
20  ns  where  frame  (a)  shows  the  cell  prior  to  the  pulse  arrival  at  0  ns,  (b)  shows  peak 
non-delayed  pulse  intensity  outside  the  cell  at  23  ns,  (c)  captures  the  peak  delayed 
pulse  intensity  at  43  ns,  and  (d)  shows  the  vapor  cell  after  the  pulse  passage  at  55  ns. 
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Figure  15.  Close  up  detail  showing  half  heater  block  with  alkali-metal 
vapor  cell.  The  PI-MAX  intensified  CCD  (ICCD)  image  array  is  shown 
as  the  extruded  rectangle  intersecting  the  vapor  cell.  Notice  the  image 
array  captures  non-delayed  pulses  outside  the  alkali-metal  vapor  cell  with 
the  expected  delayed  pulses  imaged  through  the  reference  cell. 


Figure  16.  Four  frames  taken  from  PI-MAX  intensified  CCD  (ICCD)  im¬ 
age  array.  Each  256  x  1024  image  array  frame  captures  laser  pulses  outside 
the  alkali-metal  vapor  cell  shown  at  the  top  of  each  image  and  laser  pulses 
imaged  through  the  reference  cell.  The  edge  of  the  Pyrex®  reference  cell 
is  highlighted  with  dashed  curves.  The  series  of  images  details  the  time 
evolution  of  a  10-ns  pulse  delayed  20  ns  where  (a)  shows  the  cell  prior 
to  the  pulse  arrival  at  0  ns,  (b)  shows  peak  non-delayed  pulse  intensity 
outside  the  cell  at  23  ns,  (c)  captures  the  peak  delayed  pulse  intensity 
at  43  ns,  and  (d)  shows  the  vapor  cell  after  the  pulse  passage  at  55  ns. 
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Conclusion 


Tuning  of  optical  pulse  delays  obtained  in  the  hyperfine  doublet  structure  for 
the  alkali  atoms  may  be  achieved  through  frequency  modulation.  Indeed,  delays  of 
Ons  to  37ns  have  been  observed  in  a  single  pass  through  a  5-cm  Cs  cell  at  137.2  °C 
by  detuning  from  resonance  25  GHz.  Delay  predictions  are  nearly  exact  when  the 
modeling  precisely  matches  the  absorption  profile  including  Voigt  lineshapes  and  full 
hyperfine  structure.  The  frequency  dependence  of  the  delays  are  sensitive  to  va¬ 
por  pressure,  constraining  cell  temperature  uncertainties  to  about  0.1  °C.  Using  the 
Voigt  profile  and  hyperfine  structure  allows  for  accurate  prediction  of  the  pulse  delay 
achieved  by  frequency  tuning. 
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V.  Hole  Burning 


The  full  frequency  dependence  of  the  optical  delay  in  the  Cs  D\  (62<Si/2  — 
62-Pi/2^  line  has  been  observed,  including  all  four  hyperfine  split  compo¬ 
nents.  Pulse  delays  of  1.6  ns  to  2Jhl  ns  are  obtained  by  scanning  across 
the  hyperfine  splitting  associated  with  the  lower  2Si/2  state.  Optical  con¬ 
trol  of  pulse  delays  in  cesium  vapor  were  demonstrated  by  pumping  the  D2 
(62S1/2  —  Q2P\/2 )  transition  and  observing  resulting  holes  in  the  D\  delay 
spectrum.  For  a  pump  at  four  times  the  saturation  intensity,  the  pulse 
delays  are  reduced  by  a  maximum  of  18%  in  a  narrow  region  of  110  MHz. 
The  frequency  dependence  of  the  delays  of  the  probe  laser  in  the  vicinity 
of  the  spectral  holes  agrees  with  a  Kramers- Kronig  model  prediction. 


Introduction 

Numerous  applications  from  communications  to  interferometry  for  remote  sensing 
stand  to  benefit  from  the  tunable  delay  of  optical  pulses.  Several  methods  have  been 
used  seeking  to  achieve  a  rapidly  tunable  optical  delay  system  capable  of  delaying 
pulsed  information  multiple  pulse  widths  with  little  loss  of  signal  amplitude  and  no 
pulse  dispersion.  It  is  well  understood  that  coherently  coupled  fields  can  produce  slow 
light  through  electronically  induced  transparency  (EIT)  (Fleischhauer  et  al.,  2005; 
Boiler  et  ah,  1991;  Kasapi  et  ah,  1995).  Researchers  have  exploited  an  alternate 
method  of  controlling  pulse  delays  through  the  rapidly  varying  frequency-dependent 
absorption  near  atomic  transitions  (Grischkowsky,  1973;  Tanaka  et  ah,  2003;  Camacho 
et  ah,  2006a).  We  previously  reported  controlling  optical  delay  by  tuning  a  single  laser 
across  the  D2  (6  2Si/2  —  6  2Pi/2)  absorption  spectrum  of  cesium  vapor  (Anderson  and 
Perram,  2010a).  By  extending  the  spectral  analysis  to  include  Voigt  profiles  with 
the  Doppler  component,  the  full  frequency  dependence  of  the  observed  delays  was 
adequately  described.  In  this  paper  we  extend  this  analysis  to  the  Cs  Di  (62Si/2  — 
62Pi/2)  spectrum  where  the  hyperfine  splitting  the  large  enough  to  resolve  all  four 
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components.  Agarwal  and  Dey  proposed  using  hole  burning  in  a  Doppler-broadened 
atomic  vapor  using  a  counter  propagating  pump  beam  to  saturate  the  media  and 
achieve  group  indices  of  about  103  (Agarwal  and  Dey,  2003).  Narrower  resonances  may 
be  achieved  using  quantum  coherence  effects  (Bigelow  et  al.,  2003).  Shakhmuratov 
et  al.  extended  the  technique  to  persistent  spectral  hole  burning  and  slow  light 
in  inhomogeneously  broadened  solid  materials  (Shakhmuratov  et  ah,  2005).  Most 
recently,  Camacho  et  al.  experimentally  realized  these  concepts  in  hot  rubidium 
vapor  (Camacho  et  ah,  2006b).  The  pump  laser  induces  transparency  by  moving 
population  from  the  probed  hyperfine  state  to  the  other  hyperhne  component,  not  by 
coherent  effects.  This  work  focused  on  frequency  modulation  of  the  pump  beam  to 
tune  delays  of  up  to  45  ns  (Camacho  et  ah,  2006b). 

In  this  experiment  we  use  a  continuous  wave  (CW)  D2  pump  laser  to  depopulate  a 
given  velocity  group  in  the  Cs  ground  state  and  probe  the  full  frequency  dependence 
of  the  delays  observed  in  the  Di  absorption  spectrum.  Monitoring  pulse  delays  as  a 
function  of  probe  frequency  clearly  revealed  a  pump  intensity  dependent  hole. 

Experimental  Setup 

Experimental  setup  for  hole  burning  delay  measurements  shown  in  Fig.  17  is  a 
modified  sub-Doppler  absorption  spectroscopy  experiment  where  counter  propagating 
beams  intersect  in  the  hot  cesium  vapor  cell.  A  New  Focus  Velocity™  model  6316  20- 
mW  tunable  diode  laser  with  a  linewidth  less  than  300  kHz  serves  as  the  pump  laser 
on  the  D2  Cs  transition  (852  nm).  The  pump  laser  is  directed  through  a  pellicle  beam 
splitter  where  a  beam  sample  is  sent  to  a  Bristol  model  621 A  wavelength  meter  with 
0.075  GHz  (0.0025  cm-1)  accuracy  for  frequency  reference.  The  main  pump  beam  is 
directed  through  a  5-cm-long,  low-pressure  cesium  vapor  cell  enclosed  in  a  machined 
aluminum  heater  block  maintained  at  a  constant  temperature  ±2  °C  between  40  °C 
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and  120  °C  using  a  temperature  controller  system. 


Detector 


Pump  Beam 
Laser  Splitter 


Beam  Splitter 


,  ,  .  Beam 
Laser  Isdat0r  Splitter 


Mirror 

Figure  17.  Experimental  setup  for  hole  burning  delay  measurements. 


A  New  Focus  Velocity™  model  6318  6-mW  tunable  diode  laser  with  less  than  300- 
kHz  linewidth  is  used  as  the  D\  probe  laser  at  894  nm.  The  probe  beam  is  directed 
through  an  optical  isolator  and  a  pellicle  beam  splitter.  A  probe  beam  sample  is 
sent  to  a  second  Bristol  model  621A  wavelength  meter  while  the  main  probe  beam 
continues  to  a  EOSPACE  fiber-coupled  Mach-Zehnder  interferometric  lithium  niobate 
(LiNbOs)  electro-optic  intensity  modulator  (EOM).  A  Picosecond  Pulse  Labs  1.6- 
GHz  pulse/pattern  generator  drives  the  EOM  to  create  a  Gaussian  7.5-ns  full- width 
at  half-maximum  (FWHM)  pulse  with  133  Hz  repetition  rate.  The  modulated  probe 
is  split  exiting  the  EOM  fiber  to  a  silicon  photodiode  detector  with  2.3-ns  response 
time  providing  initial  intensity  reference  before  the  alkali-metal  vapor  cell.  The  probe 
beam  is  aligned  overlapping  but  counter  propagating  to  the  pump  beam  in  the  Cs 
cell  enclosed  in  the  heater  block.  Probe  pulse  delay  is  measured  after  the  alkali-metal 
vapor  cell  by  a  silicon  photodiode  detector  with  a  1-ns  response  time. 

Pulse  delay  temporal  waveforms  were  2000  pulse  averages  recorded  on  a  3  GHz 
oscilloscope  triggered  by  the  pulse  generator  signal.  The  GW  D2  pump  is  set  at  a 
fixed  power  and  frequency  detuned  from  an  absorption  peak.  Neutral  density  Liters 
placed  in  the  D2  pump  beam  path  to  vary  pump  power  at  0  mW  (blocked),  5  mW 
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and  16  mW.  Probe  laser  is  scanned  across  the  Di  absorption  spectrum  near  F"  =  4 
to  F'  =  3,4,5  as  depicted  on  the  Cs  energy  level  diagram  in  Fig.  18.  Simulated 
absorption  coefficients  for  Cs  D\  and  D2  at  50  °C  are  shown  in  Fig.  19.  The  x-axis 
is  labeled  using  both  nanometer  (nm)  and  detuning  frequency  v&i  and  z/A2  (GHz) 
scales.  The  detuning  frequencies  and  z/A2  represent  detuning  from  the  central 
atomic  transition  frequency  335.116049  THz  for  Di  and  351.725718  THz  for  D-2. 
Regions  of  the  absorption  spectra  used  in  this  experiment  are  highlighted  with  gray 
squares. 


Figure  18.  Energy  level  diagram  depicting  the  pump  D 2  transitions  and 
probe  D 1  transitions  Energy  level  separations  are  not  drawn  to  scale. 


Optical  Delay 


Previous  work  estimated  pulse  delay  using  a  group  velocity  through  the  Kramers- 
Kronig  relations  to  find  the  real  index  of  refraction  rif{  from  the  absorption  coefficient 
a  (Anderson  and  Perram,  2010a, b)  such  that 


nR  =  1  + 


<*(0 

e-v2 


dt 


(99) 
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852.360  nm  852.355  nm  852.336  nm 


Figure  19.  Simulated  absorption  coefficients  for  the  (a)  D\  and  (b)  £>2 
transitions  at  50 °C.  Single  hyperfine  transition  absorption  lineshapes 
are  shown  in  the  D2  spectrum  as  thin  black  lines.  Regions  of  the  D- 
line  absorption  spectra  used  in  this  experiment  are  highlighted  by  gray 
squares. 


Optical  pulse  delays  are  calculated  as  td  =  l/vg  for  cell  length  l  where  the  group 
velocity  is  defined  as 


V9  = 


(100) 


Observed  D\  pulse  delays  compared  favorably  to  previous  work  in  Cs  D2  at  similar 
temperatures.  Kramers- Kronig  delay  model  prediction  and  observed  delays  for  a  Cs 
vapor  cell  at  109.5  °C  are  shown  in  Fig.  20  where  the  solid  line  represents  the  delay 
model  prediction.  The  Cs  D 1  on  resonance  optical  thickness  £  =  l/a(u0)n  ~  0.05  mm 
for  the  F"  =  4  — »  F'  =  3  transition  z/ai  =  —4.665  GHz  at  109.5  °C. 

The  pulse/pattern  generator  produced  a  nearly  Gaussian  7.5-ns  pulse  as  shown 
in  Fig.  21.  Sample  pulses  at  three  detuning  frequencies  z^ai  ate  depicted  as  a 
series  of  gray  shapes  where  (1)  is  farthest  from  the  absorption  peak  while  pulse  (3) 
is  closest  to  the  resonant  transition  frequency.  Solid  black  lines  show  Gaussian  fit 
curves,  and  the  peak  of  the  fit  curve  is  marked  for  reference.  As  expected  the  pulse 
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Detuning  Frequency,  vA]  (GFIz) 

Figure  20.  Pulse  delay  for  full  spectrum  Cs  D\  scan  from  i/^i  =  -8  GHz 
to  +8  GHz.  Cesium  cell  was  maintained  at  109.5  °C.  No  £>2  pump  was 
used. 

amplitude  decreases  for  pulses  closer  to  the  peak  absorption  while  the  pulse  width 
increases  due  to  dispersion.  Peak  time  for  each  Gaussian  fit  is  recorded  as  the  pulse 
arrival  time.  Zero-delay  time  is  determined  when  the  probe  is  detuned  far  off  resonant 
frequency  and  the  value  is  subtracted  from  the  arrival  time  for  each  pulse  to  calculate 
the  pulse  delay.  Maximum  error  in  delay  time  measurement  is  ±0.15  ns  with  average 
error  ±0.03  ns.  Near  the  resonant  transitions  resolving  the  peak  is  difficult  as  the 
pulse  amplitude  is  attenuated  approaching  detector  noise  level.  Pulse  peaks  lost  in 
the  noise  appear  in  Fig.  20  as  ‘stray’  data  points  below  the  model  prediction  curve. 

As  demonstrated  in  Fig.  21,  pulses  generally  decrease  in  amplitude  near  absorp¬ 
tion  peaks  and  the  pulse  FWHM  increase.  Dispersion  of  the  Gaussian  fit  is  shown 
in  Fig.  22.  As  the  probe  is  tuned  toward  the  transition  frequency  with  no  pump 
laser,  pulses  initially  appear  to  broaden  before  pulse  narrowing  is  briefly  observed. 
At  higher  pump  intensities,  pulse  dispersion  appears  more  evident  near  resonant  fre¬ 
quency  peaks,  however  narrowing  is  not  observed  before  pulse  attenuation  prevents 
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Figure  21.  Pulse  shape  for  selected  pulses  from  data  depicted  in  Fig.  20 
at  109.5  °C.  Pulse  data  for  pulse  at  three  D\  detuning  frequencies  uai  are 
depicted  using  gray  lines.  Solid  black  lines  show  Gaussian  fit  curves,  and 
the  peak  of  the  fit  curve  is  marked  for  reference.  (1)  =  7.44  GHz, 

peak  at  4.66  ns,  width  7.37  ns  (2)  i/^i  —  6.73  GHz,  peak  at  12.04  ns, 
width  8.60  ns  (3)  i/ai  =  6.5  GHz,  peak  at  19.84  ns,  width  11.83  ns. 


accurate  resolution  of  the  pulse. 

High-temperature  delay,  where  Doppler-broadened  absorption  masks  the  individ¬ 
ual  hyperhne  D\  peaks  as  in  Fig.  20,  compares  favorably  to  previous  work(Anderson 
and  Perrarn,  2010a).  This  work  reports  higher  frequency  sampling  with  200  frequency 
samples  over  20  GHz  compared  to  previous  work  in  Cs  D2  90  samples  over  20  GHz 
(Anderson  and  Perram,  2010a).  The  full  frequency  dependence  of  the  optical  delays 
for  the  Cs  D2  line,  similar  to  that  now  reported  for  the  D\  line  in  Fig  20,  have  been 
successfully  modeled  as  expressed  by  equations  99  and  100.  However,  the  hyperhne 
structure  of  the  excited  2P3/2  state  was  unresolved  and  delays  at  frequencies  between 
upper  state  hyperhne  components  has  not  previously  been  investigated.  The  spectral 
simulation  illustrated  in  20  agrees  very  well  with  the  observed  pulse  delays  in  this 
region. 
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Figure  22.  Pulse  width  in  the  neighborhood  of  D i  absorption  peaks. 
Gaussian  pulse  fit  FWHM  is  shown  for  0  mW  (no  pump)  and  5.6  mW 
across  the  D\  spectrum  between  i/^i  =  -6  GHz  and  -2.5  GHz. 


Intensity-Dependent  Delay 

Crossing  laser  beams  in  a  modified  saturation  spectroscopy  setup  is  a  simple  way 
to  reduce  Doppler  broadening  of  spectral  lines  and  examine  hyperfine  structure  near 
the  natural  broadening  limit.  Given  a  pump  at  frequency  vpurnp  and  D2  hyperfine 
transition  at  vm-,  the  beams  interact  at  similar  velocity  groups  determined  by  the 
frequency  difference  5vpump  =  VpUrnp  —  uori-  Spectral  holes  are  burnt  at  frequency 
Vhoie  such  that  Vhoie  =  iyDi  —  fivpump  for  each  hyperfine  transition.  For  beams  that 
are  not  perfectly  aligned  to  be  counter  propagating  and  intersecting,  there  will  be 
some  residual  Doppler  component.  The  depth  and  width  of  the  holes  are  intensity 
dependent  where  hole  width  A Vh0ie  is  a  function  of  intensity  and  saturation  intensity 
I sat  such  that 

/  /  \  V2 

A  Vhoie  0C  (  1  +  —  .  (101) 

Overall  D\  absorption  decreases  as  D2  pump  laser  power  increases  due  to  ground 
state  depletion  and  spectral  holes  are  evident  in  the  Di  absorption.  The  absorption 
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Figure  23.  Absorption  coefficients  for  Cs  D i  with  Cs  D 2  pump  at  1^2  = 

-3.84  GHz  (852.3566  nm,  11732.1790cm-1)  various  pump  powers. 

spectra  for  a  low-pressure  Cs  cell  at  44  °C  with  a  D2  pump  at  a  fixed  frequency  v\2 
=  -3.84  GHz  (852.3566  nm,  11732.1790  cm-1)  at  0  mW  (no  pump),  5.7  mW  and 
16.4  111W  are  shown  in  Fig.  23.  At  the  resonant  frequency,  the  D\  laser  is  2.3% 
I sat  while  D‘2  pump  intensity  is  approximately  4  x  Isat  at  5.6  111W  and  12  x  Isat  at 
16.4  mW. 

In  the  presence  of  holes  burnt  in  the  D 1  absorption  spectra,  predicted  delay  is 
calculated  from  observed  absorption  spectra  using  the  Kramers-Kronig  model  to  cal¬ 
culate  group  velocity  and  pulse  delay.  Observed  delay  and  the  predicted  delay  for  Cs 
at  52.5  °C,  z^A2  =-3.84  GHz  at  three  D2  pump  powers  are  shown  in  Fig.  24.  Pump 
power  is  adjusted  using  neutral  density  filters  as  changing  diode  current  causes  a 
subsequent  shift  in  pump  laser  frequency. 

Rapid  variations  in  predicted  delay  near  absorption  peaks  in  Fig. 24  are  due  to 
frequency  instability  of  the  pump/probe  lasers  in  a  region  of  highly  dispersive  ab¬ 
sorption.  In  general  the  D\  optical  delay  decreased  with  increasing  D2  pump  power. 
It  is  speculated  that  this  may  provide  a  rapidly  tunable  optical  delay  mechanism 
by  varying  pump  power  though  it  has  not  been  verified  how  quickly  the  system  can 
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respond  to  changes  in  pump  intensity. 


Hole  Burning  Effects 

Sub-Doppler  spectroscopy  removes  some  of  the  velocity-induced  broadening  using 
counter  propagating  beams  interacting  with  atoms  moving  within  a  narrow  velocity 
class.  For  a  single  resonant  transition  at  z/0  a  pump  beam  detuned  +5v  interacts 
with  atoms  at  vo  +  5v  effecting  the  observed  spectrum  of  a  counter  propagating  probe 
beam  at  u0  —  5u.  In  cesium  D2  there  are  three  closely  spaced  hyperhne  components 
corresponding  to  the  allowed  transitions  of  the  2P3/2  energy  splitting  as  shown  in 
Fig.  18  and  Fig.  19.  A  D2  pump  interacts  with  velocity  groups  detuned  +<5zy  from 
each  hyperhne  component.  The  counter  propagating  D1  probe  sees  spectral  holes  at 
corresponding  velocity  groups  —  from  each  D\  hyperhne  component. 

In  the  neighborhood  of  the  spectral  hole  absorption  changes  rapidly  across  a  nar¬ 
row  frequency  range.  The  Kramers-Kronig  delay  model  predicts  a  significant  change 
in  delay  in  the  hole  region.  To  observe  the  effect  of  a  spectral  hole  on  the  pulse  delay, 
the  D2  pump  laser  was  set  at  z/A2  =-3.86  GHz  for  16.4  mW  and  52.5  °C  while  the  D\ 
probe  was  scanned  across  the  absorption  spectrum  from  =-6  GHz  to  -2  GHz. 

The  observed  transmittance  is  shown  in  Fig.  25  with  a  no  pump  reference  absorp¬ 
tion  spectrum  (dashed  lines)  to  visualize  the  spectral  holes.  A  simulated  absorption 
spectrum  was  then  created  using  the  observed  spectrum  as  a  reference  by  subtracting 
a  series  of  Lorentzian  line  shapes  from  the  reference  absorption  spectrum.  The  re¬ 
sulting  difference  is  in  the  absorption  spectrum  is  shown  in  Fig.  25  as  the  solid  lines. 
Lorentzian  line  shapes  used  in  the  prediction  have  frequency  differences  of  192  MHz 
and  285  MHz  compared  to  accepted  hyperhne  splitting  of  approximately  201  MHz 
and  251  MHz.  The  Lorentzian  FWHM  for  the  absorption  corrections  is  110  MHz 
which  is  more  than  20  times  the  natural  linewidth  for  Cs  Di  and  D2  and  less  than 
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Detuning  Frequency,  VA1  (GHz) 


Figure  24.  Observed  delay  (black  dots)  and  the  predicted  delay  (gray 
lines)  for  Cs  at  52.5  °C,  ^2  =  -3.84  GHz  at  three  D?  pump  powers  of 
(a)  0  mW  (no  pump),  (b)  5.7  mW  and  (c)  16.4  mW.  Predicted  delay  at 
each  power  is  calculated  using  the  observed  absorption  spectra  and  the 
Kramers-Kronig  relations  model. 
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Figure  25.  Observed  Cs  D\  transmittance  with  reference  absorption 
spectrum  (dashed  lines)  to  visualize  the  spectral  holes.  Difference  in  the 
absorption  spectrum  is  shown  as  solid  lines. 

the  Doppler  width  (~400  MHz).  Frequencies  used  for  each  Lorentzian  line  shape  used 
to  approximate  the  reference  absorption  spectrum  are  on  the  order  of  the  hyperhne 
separation  of  the  Cs  D-2  F"  =  4  — >•  F'  =  3, 4,  5  transitions.  The  modified  absorption 
curve  similar  to  the  observed  absorption  spectrum  is  then  used  to  predict  the  optical 
delay  using  the  Kramer-Kronig  model  and  the  resulting  delay  prediction  and  observed 
delay  is  given  in  Fig.  26.  The  observed  delay  in  the  region  of  the  hole  match  the 
Kramer-Kronig  model  prediction  in  general  shape  and  magnitude. 

Conclusion 

The  frequency  dependence  of  the  optical  pulse  delay  in  cesium  D \  closely  match 
the  previously  developed  model  predictions  using  Voigt  lineshapes  and  the  full  hyper- 
fine  structure.  The  observed  D\  optical  pulse  delays  compared  favorably  to  similar 
results  in  Cs  D2.  The  delay  prediction  Kramers-Kronig  model  accurately  predicts 
the  delays  between  Di  upper  state  hyperhne  components  at  lower  temperatures  than 
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Figure  26.  Pulse  delay  in  the  narrow  spectral  region  surrounding  a  hole 
burnt  in  the  Cs  D\  absorption  spectrum.  Solid  line  depicts  the  model 
delay  prediction. 

previously  reported  with  smaller  frequency  sampling  separation.  Delays  up  to  24  ns 
where  achieved  by  detuning  from  resonance  by  10  GHz. 

Tunable  optical  delays  were  demonstrated  with  a  Cs  D i  pulse  by  varying  D2  pump 
laser  power.  Increasing  Do  pump  power  reduced  D\  absorption  overall.  Observation 
of  hole  burning  effects  on  delay  close  to  the  spectral  hole  were  noted.  Location  of 
the  holes  and  predicted  delay  from  a  modified  absorption  spectrum  matched  model 
predictions.  For  a  pump  intensity  of  four  times  saturation  intensity  Isat  observed 
delay  is  reduced  by  an  average  of  50%  across  the  spectrum  and  reaches  a  maximum 
of  78%  in  a  narrow  110  MHz  spectral  hole  region. 
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VI.  Hyperfine  Relaxation 


A  Cesium.  D2  (62Si/2  —  6 2P^/2)  probe  laser  is  used  to  examine  temporal 
pulse  dynamics  during  a  10 -ns  high-power  D2  pulse  event.  The  probe  laser 
is  set  perpendicular  to  the  pump  pulse  in  a  heated  Cs  optical  cell  with  he¬ 
lium  buffer  gas.  The  probe  frequency  is  scanned  across  the  D2  transitions 
and  the  dynamic  excitation  and  relaxation  of  the  hyperfine  components 
is  studied  relating  to  temperature,  buffer  gas  pressure,  and  pump/probe 
intensity. 


Introduction 

Alkali-metal  vapor  laser  systems  were  available  shortly  after  the  laser’s  invention 
using  cesium  vapor  6Si/2  — »  8Pi/2  pumped  at  388.8  nm  (Rabinowitz  et  ah,  1962). 
Alkali-metal  elements  have  a  simple  atomic  structure  with  a  single  valence  electron 
and  are  relatively  simple  to  model.  The  quantum  defect  can  be  quite  small  leading 
to  highly  efficient  laser  systems  such  as  diode-pumped  alkali  laser  (DPAL)  systems 
(Krupke  et  ah,  2003).  Diode-pumped  alkali  lasers  are  excellent  three- level  lasers 
provided  the  Di  and  D2  absorption  profiles  are  sufficiently  broadened  to  match  pump 
linewidth  and  population  excited  to  the  2P3/2  level  can  be  quickly  relaxed  to  the  upper 
lasing  2P 1/2  state.  Broadening  is  typically  achieved  using  a  homogeneous  pressure 
broadening  species  such  as  helium  and  a  collisional  P3/2  — >  P1/2  relaxation  gas  such 
as  ethane  (Krupke  et  ah,  2003).  DPAL  systems  have  been  successfully  demonstrated 
in  Rb  (Krupke  et  ah,  2003),  Cs  (Beach  et  ah,  2004;  Ehrenreich  et  ah,  2005)  and  other 
alkali-metal  vapors  building  hopes  that  DPAL  systems  offer  a  scaling  path  to  high 
power  laser  systems  without  the  thermal  management  issues  of  solid  state  lasers  or 
the  hazards  of  chemical  laser  systems. 

Modeling  of  the  intra-cavity  dynamics,  such  as  thermal  effects  and  bleached  am¬ 
plification  medium,  of  a  DPAL  system  and  subsequent  validation  of  the  theoretical 
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description  are  critical  to  power  scaling  and  maximizing  laser  efficiency.  Optical  mea¬ 
surement  of  the  bleach  wave  dynamics  have  not  been  previously  accomplished.  Except 
for  a  very  small  increase  in  probe  intensity  at  the  absorption  peak,  the  methodology 
described  in  this  section  did  not  prove  successful  in  imaging  the  bleach  wave,  however 
other  significant  pulsed-pump  hyperffiie  dynamics  were  observed. 

Direct  observation  of  hyperffiie  relaxation  is  difficult  due  to  the  hyperffiie  filters  re¬ 
quired  to  observe  resonance  lines  (Beverini  et  al.,  1971).  Typical  hyperffiie  relaxation 
rates  are  measured  for  alkali- foreign-gas  collision-induced  relaxation  of  the  various 
ground  state  levels  (Beverini  et  al.,  1971;  Bernabeu  and  Alvarez,  1980;  Tornos  and 
Amare,  1986),  normally  accomplished  at  low  temperatures  to  reduce  Cs-Cs  collisions. 
Wall  collisions  and  diffusion  are  predominant  factors  besides  the  intended  foreign- 
gas  interactions.  In  previous  work  hyperffiie  relaxation  of  the  moderately  spaced  Di 
levels  with  a  collisional  gas  were  studied  (Morgus  et  ah,  2008)  using  counter  prop¬ 
agating  pump  beams  at  Cs  Di  and  to  excite  6Pi/2  atoms  to  the  8Si/2  level  with  a 
761.1  nm  pump  laser.  Relaxation  was  observed  measuring  side  fluorescence  of  the 
8Si/2  — »  6P3/2  photons. 

To  the  best  of  my  knowledge  pump-probe  study  of  the  hyperffiie  relaxation  related 
to  pulsed  pump  rates  for  DPAL  systems  has  not  been  accomplished.  This  work 
examines  pulsed  pump  dynamics  of  a  DPAh-like  cesium  vapor  cell  related  to  the 
excitation  and  relaxation  to  pre-pulse  equilibrium  of  hyperffiie  components. 

Experimental  Setup 

The  experimental  setup  to  measure  pulse  dynamics  is  shown  in  Figure  27.  A 
Spectra  Physics™  Quanta  Ray  PRO-290  Nd:YAG  laser  is  used  as  a  pump  laser  for 
the  Sirah™  Precision  Scan  PRSC-D-18  pulse  dye  laser  (LDS  867  dye  mixture)  which 
produced  a  10  ns  pulse  with  800  KJ  per  pulse  at  10  Hz  and  a  spectral  linewidth 
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Figure  27.  Setup  for  hyperfine  relaxation  experiments.  A  10  Hz  Sirah 
dye  laser  generates  high-power  pump  pulses  sent  through  the  Cs  vapor 
cell.  A  low-power  D2  probe  passes  perpendicular  to  the  pump  pulse. 

of  1.4  GHz  at  852.34  nm.  For  comparison  closely  spaced  group  of  three  hyperfine 
lines  at  Cs  D2  v&2  ~  —4  GHz  has  a  Doppler-broadened  FWHM  of  0.8  GHz  so  that 
the  pump  pulses  cover  the  full  set  of  three  Cs  D2  upper  state  hyperfine  split  energy 
levels.  The  pulsed  pump  beam  is  directed  through  a  series  of  mirrors  and  an  aperture 
to  limit  the  spatial  pump  beam  width  to  the  diameter  of  the  cesium  cell.  The  2.5  cm 
beam  passes  through  the  heated  alkali  cell  and  into  a  beam  block.  Care  was  taken  to 
eliminate  reflections  and  scattered  pump  light  with  the  extensive  use  of  black  foam 
board  and  blackout  curtains. 

A  New  Focus  Velocity™  model  6316  20-mW  tunable  diode  laser  with  a  linewidth 
less  than  300  kHz  serves  as  the  continuous  wave  probe  laser  on  the  D2  Cs  transi¬ 
tion  (852  nm).  The  probe  laser  is  directed  through  a  pellicle  beam  splitter  where 
a  beam  sample  is  sent  to  a  Bristol  model  621A  wavelength  meter  with  0.075  GHz 
(0.0025  cm-1)  accuracy  for  frequency  reference.  The  probe  beam  is  directed  through 
the  alkali-metal  vapor  cell  perpendicular  to  the  pump  beam  and  is  then  directed  along 
a  5  m  path  through  two  apertures  to  eliminate  pump  beam  reflections  and  side  flu¬ 
orescence  on  the  way  to  the  detector.  A  1-ns  rise  time  silicon  photodiode  is  used  to 
collect  probe  beam  intensity  measurements  recorded  on  a  3  GHz  oscilloscope. 


Results 


The  pump  laser  output  50  mW  pulsed  power  at  the  cesium  cell  and  was  set  to 
pump  the  D2  absorption  feature  peaks  at  (1)  852.357  nm  (ua2  =-4.  GHz)  and  (2) 
852.335  nm  (va2  =+5  GHz)  as  shown  in  Figure  28.  Points  A,  B,  C  and  D  refer  to 
frequencies  detuned  from  the  absorption  peaks  until  the  probe  beam  transmission 
is  approximately  50%.  Hyperfine  F"  — >•  F'  transitions  are  labeled  in  Figure  28  for 
reference. 


Detuning  Frequency,  vA2  (GHz) 

Figure  28.  Cesium  D2  absorption  with  hyperfine  structure.  Reference 
locations  A,  B,  C  and  D  are  set  at  50%  probe  transmission  and  pump 
locations  are  labeled  1  and  2  at  the  absorption  peaks. 


It  was  hoped  that  the  pump  beam  would  bleach  the  Cs  vapor  and  the  probe  beam 
would  pass  unattenuated  around  the  pulse  arrival  time  so  that  temporal  and  spatial 
observation  of  the  DPAL  bleach  wave  could  verify  theoretical  results.  While  this 
was  not  the  case,  other  interesting  frequency-dependent  hyperfine  relaxation-related 
phenomena  were  observed. 

There  are  distinct  short-term  and  long-term  effects  following  a  pump  pulse  event. 
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Sample  pump  and  probe  intensities  are  shown  in  Figure  29.  Pulse  data  is  shown 


Figure  29.  Sample  pump  and  probe  intensities  at  the  detector.  These 
sample  pulses  taken  for  probe  Point  A  at  80  °C  illustrate  (a)  the  short¬ 
term  nanosecond  timescale  effects  in  the  neighborhood  of  the  pump  pulse 
and  (b)  long-term  microsecond  hyperflne  effects  and  relaxation  to  pre¬ 
pump  equilibrium. 


for  probe  frequency  at  point  A  (852.357  nm),  pump  at  frequency  (2)  852.335  nm 
(ua 2  =+5  GHz),  Cs  cell  temperature  T  =  80 °C.  The  time  axis  is  synchronized  to 
the  q-switch  signal  from  the  Quanta  Ray  Nd:YAG  laser.  Arrival  time  for  the  pump 
at  140  ns  is  shown  as  the  gray  line.  Figure  29  (a)  shows  short-term  effects  on  probe 
intensity,  dropping  initially  during  pulse  arrival,  partially  returning  before  decreasing 
again  past  200  ns.  Figure  29  (b)  shows  long-term  probe  intensity  decay  to  return  to 
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pre-pulse  equilibrium  around  50  ps. 

Full  frequency  spectra  from  v&2  =  ±7  GHz  illustrating  short-term  and  long¬ 
term  effects  are  shown  in  Figure  30  through  Figure  33  for  a  pump  at  frequency 
(2)  852.335  nm  (u^2  =+5  GHz).  In  these  figures,  the  absorbance  spectrum  plot  is 
interpreted  such  that  lower  values  indicate  little  absorption  so  that  high  probe  signal 
intensity  is  transmitted. 

Figure  30  shows  nanosecond-scale  pulsed  pump  effects  as  the  probe  is  scanned  in 
frequency  across  the  Cs  D2  absorption  spectrum.  Figure  31  shows  the  absorbance  at 
four  distinct  frequencies  (A,  B,  C  and  D)  from  Figure  30  (see  Figure  28  for  frequency 
reference  points). 


Detuning  Frequency  (GHz) 


Figure  30.  Nanosecond-scale  frequency-dependent  probe  intensity  after 
a  pump  pulse  event  at  +5  GHz. 


Figure  32  shows  microsecond-scale  pulsed  pump  effects  as  the  probe  is  scanned  in 
frequency  across  the  Cs  D2  absorption  features.  Figure  33  shows  the  absorbance  at 
four  distinct  frequencies  (A,  B,  C  and  D)  from  Figure  32  (see  Figure  28  for  frequency 
reference  points). 
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Figure  31.  Nanosecond-scale  frequency-dependent  probe  intensity  after 
a  pump  pulse  event  at  +5  GHz  at  Points  A,  B,  C,  and  D  (See  Figure  28 
for  graphical  depiction  of  reference  points). 


De<“™ngFrequency(GH;;) 


Figure  32.  Microsecond-scale  frequency-dependent  probe  intensity  after 
a  pump  pulse  event  at  +5  GHz. 
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Figure  33.  Microsecond-scale  frequency-dependent  probe  intensity  after 
a  pump  pulse  event  at  +5  GHz  at  Points  A,  B,  C,  and  D  (See  Figure  28 
for  graphical  depiction  of  reference  points). 
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For  the  long-term  effect  each  pulse  is  fit  to  a  double  exponential  of  the  form 


y  =  a  +  6exp(— cx)  +  dexp(— ex). 


(102) 


Detector  intensity  data  is  transformed  using 


V  =  ~  In 


I 

T0 


(103) 


where  Jo  is  the  pre-pulse  intensity.  A  sample  pulse  data  plot  with  double  exponential 
fit  is  shown  in  Figure  34  for  a  probe  at  point  B  with  21.4  mW  probe  power.  Data 


0  50  100  150 

Time  (ps) 


Figure  34.  Sample  data  fit  is  shown  for  a  probe  at  point  B  with  21.4  mW 
probe  power.  Data  before  the  peak  is  excluded  from  the  fit  data  and  the 
fit  function  is  shown  as  a  solid  black  line..  Data  for  all  samples  is  fit  to 
a  double  exponential  function  of  the  form  y  =  a  +  bexp(-cx)  +  dexp(— ex). 
Residuals  from  the  fit  are  shown  below  the  data  plot. 


before  the  peak  is  excluded  from  the  fit  data,  the  resulting  fit  function  is  shown  as  a 
solid  black  line  over  the  pulse  data  and  residuals  are  shown  below  the  data  plot. 
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Data  Trends 


It  should  be  noted  in  that  in  Figure  30  the  short-term  timescale  pnlse  effects  at  the 
time  of  the  pnlse  are  observed  far  from  the  absorption  peak,  especially  near  point  A, 
the  Cs  D-2  F"  =  4  — y  F'  =  3, 4  transitions.  Little  pnlse  effect  is  noted  around  point  B 
in  Figure  30  and  the  shape  of  the  pump  pulse  is  dramatically  different  compared  the 
response  at  other  frequencies.  There  is  a  small  region  of  increased  intensity  (decreased 
absorption)  at  the  —4  GHz  absorption  peak  with  the  pump  tuned  to  peak  2  at  5  GHz. 

Similar  hyperhne  frequency-dependent  pulse  response  and  double-exponential  re¬ 
turn  to  equilibrium  is  evident  in  the  long-term  microsecond  scale  Figures  32  and  33. 
Decay  curves  in  the  A,  C,  and  D  regions  are  similar  while  the  decay  around  B  is 
markedly  longer  lived  at  80  °C. 

Recall  some  of  the  general  trends  in  the  absorption  coefficient.  As  temperature 
increases,  alkali  concentration  and  inhomogeneous  Doppler-broadening  increase.  Ho¬ 
mogeneous  broadening  increases  with  buffer  gas  pressure.  Pressured  broadened  Cs 
D2  with  He  is  broadened  sufficiently  that  the  50%  transmission  “saddle”  between  the 
ground  state  splitting  is  no  longer  distinguishable  before  15  Torr  He  at  90  °C.  There¬ 
fore  at  higher  temperatures  and  buffer  gas  pressures  all  four  frequency  points  are  no 
longer  distinguishable.  Therefore  it  is  only  possible  to  compare  points  A  and  B  at 
low  pressures  and  temperatures.  Figure  35  shows  the  double  exponential  fit  decay 
constants  c  and  e  for  points  A  and  B  at  various  HE  pressures  up  to  12  Torr  at  93.5  °C. 
It  it  evident  that  the  constants  rapidly  approach  constant  values  with  increasing  He 
pressure. 

Temperature  dependence  for  points  A  and  B  data  in  Figure  36  shows  that  the  c 
and  e  values  converge  to  a  single  value  indicating  that  the  double-exponential  decay 
curve  converges  to  a  single-exponential  decay. 

Closer  examination  of  the  He  pressure  dependence  is  shown  in  Figure  37.  As  He 
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Figure  35.  Decay  rates  c  and  e  (from  the  double  exponential  fit  function 
y  =  a  +  bexp(-cx)  +  dexp(-ex))  for  probe  frequencies  A  and  B  are  shown 
at  various  He  pressures  (Torr). 
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Figure  36.  Decay  rates  c  and  e  (from  the  double  exponential  fit  function 
y  =  a  +  bexp(-cx)  +  dexp(-ex))  for  probe  frequencies  A  and  B  are  shown 
at  low  He  pressures  (0.06  Torr). 
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Figure  37.  Decay  rates  c  and  e  (from  the  double  exponential  fit  function 
y  =  a  +  6exp(— cx)  +  dexp(— ex))  for  probe  frequency  A  are  shown  at  low 
He  pressures  (0-1.2  Torr). 

pressure  increases  the  c  decay  rate  peaks  around  1  Torr  He  for  point  A  at  93.5  °C  and 
begins  to  decrease  to  the  high-pressure  value.  The  e  decay  constant  remains  relatively 
constant. 

The  decay  constants  quickly  approach  an  asymptote  at  higher  pressures  as  shown 
in  Figure  38. 

Decay  rate  as  a  function  of  pump  and  probe  laser  power  is  also  examined  in 
Figures  39  and  40  with  a  low-pressure  alkali  sample  at  95.5  °C  and  pump  laser  at 
peak  2  (5  GHz).  Decay  rate  as  a  function  of  probe  power  is  shown  in  Figure  40. 

Discussion 

Several  pump  and  probe  characteristics  were  also  examined.  Various  linear  and  cir¬ 
cular  pump  polarization  combinations  were  used  and  no  polarization-induced  changes 
could  be  identified.  Pump  intensity  has  the  intuitive  effect  in  that  excitation  increases 
with  pump  power  (no  pump  shows  no  excitation  or  relaxation).  Probe  power  depen- 
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Figure  38.  Decay  rates  c  and  e  (from  the  double  exponential  fit  func¬ 
tion  y  =  a  +  bexp(-cx)  +  dexp(-ex))  for  probe  frequency  A  are  shown  at 
moderate  He  pressures  (0-400  Torr). 
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Figure  39.  Decay  rates  c  and  e  (from  the  double  exponential  fit  function 
y  =  a  +  bexp(—cx)  +  dex p(— ex))  for  probe  frequency  B  are  shown  at  low 
He  pressures  (0.06  Torr)  at  various  pump  powers. 
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Figure  40.  Decay  rates  c  and  e  (from  the  double  exponential  fit  function 
y  =  a  +  bexp(-cx)  +  dexp(-ex))  for  probe  frequency  B  are  shown  at  low 
He  pressures  (0.05  Torr)  at  various  probe  powers. 

dence  was  not  evident.  Decay  rates  were  constant  over  useful  range  of  probe  powers 
sufficient  to  penetrate  the  optically  thick  alkali-metal  vapor  and  produce  detector 
signal  outside  noise  limitations. 

The  full- frequency  spectra  of  Figures  30  through  33  illustrate  the  pump  pulse 
induced  effect  on  hyperhne  state  excitation  and  relaxation  return  to  equilibrium. 

For  short-term  nanosecond  scale  effects,  there  are  several  potential  explanations 
for  the  observed  decrease  in  probe  intensity  at  pump  arrival  including  spatial  vari¬ 
ations  in  temperature,  pressure  caused  by  the  pump  laser  inducing  changes  in  the 
index  of  refraction  reducing  the  probe  beam  incident  on  the  detector,  indicated  by  a 
rise  to  the  long-term  maximum  on  the  time  scale  of  the  Cs  D2  lifetime  sfter  the  pump 
pulse.  To  test  this  theory  a  fast  frame  rate  (200k  frames/second)  camera  was  setup  to 
capture  images  of  probe  spot  movement  as  the  30-mW  pump  laser  is  focused  through 
the  alkali  vapor  cell  at  80  °C  as  shown  in  Figure  41.  Beam  deflection  was  observed 
up  to  1  mm  at  60  cm  after  the  alkali  vapor  cell,  which  would  translate  to  over  3  mm 
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Figure  41.  Experimental  setup  for  beam  steering  observation.  Move¬ 
ment  of  the  probe  beam  is  captured  using  a  high  speed  (200k  frame/sec¬ 
ond)  camera  as  the  pulsed  pump  beam  is  focused  through  the  alkali-metal 
vapor  cell. 


deflection  at  the  approximately  2-nt  path  length  used  in  the  experimental  setup  de¬ 
picted  in  Figure  27.  In  such  cases  the  probe  beam  would  be  steered  completely  away 
from  the  1  mm  detector  for  deflections  of  this  magnitude.  Observed  deflections  lasted 
up  to  100  ps.  Deflections  were  most  evident  when  the  probe  beam,  detuned  from 
cesium  D2  peak  (-4  GHz),  was  closely  aligned  spatially  to  the  tightly  focused  pump 
beam  tuned  to  the  D2  peak  absorption  (+5  GHz).  Repetition  of  the  probe  beam 
deflection  roughly  followed  the  10-Hz  pulsed  pump  rate,  though  deflection  was  not 
consistent  in  magnitude  or  direction  from  pulse  to  pulse.  This  may  have  resulted 
from  inconsistent  pulse-to-pulse  pump  intensity  or  frame  synchronization  artifacts  as 
well  as  physical  variations  within  the  vapor  cell. 

Frequency  dependence  of  the  absorption  shows  there  is  a  hyperhne  component 
to  the  decay  rate’s  return  to  equilibrium.  In  the  cases  illustrated  the  F'  —  5  upper 
state,  the  primary  transition  for  frequency  point  B,  is  electronically  forbidden  from 
the  F"  =  3  pump  tuning.  The  deviation  of  the  point  B  excitation  and  decay  from 
the  other  frequency  regions  is  expected. 
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Temperature  and  He  pressure  dependence  is  predicted  as  the  total  Voigt  lineshape 
broadening  increases  transfer  between  upper  state  components.  The  F'  —  4  — >  F'  —  5 
energy  separation  is  only  251  MHz.  At  90  °C  the  Doppler-dominated  Voigt  FWHM  is 
420  MHz.  Over  the  range  of  temperatures  used  in  this  experiment  the  Voigt  FWHM 
ranges  from  410  MHz  to  440  MHz  at  0  Torr  He. 

Pressure  broadening  with  He  has  a  pronounced  effect  on  the  total  broadening  of 
the  hyperhne  lines.  At  low  pressures  the  homogeneous  broadening  is  natural  broad¬ 
ening  limited  at  5.25  MHz.  As  pressure  increases  to  0.24  Torr  the  pressure  induced 
broadening  is  equal  to  the  natural  linewidth.  By  18.6  Torr  He  at  90  °C  the  total  homo¬ 
geneous  broadening  exceeds  the  Doppler-broadening  and  the  total  Voigt  broadening 
is  680  MHz.  The  total  separation  of  the  Cs  2P3/2  upper  state  is  603  MHz.  Laser  exci¬ 
tation  should  be  well  mixed  and  the  frequency  dependence  is  reduced.  Therefore  the 
asymptotic  relaxation  rates  at  high  He  pressures  is  expected  as  depicted  in  Figure  35 
where  the  c  decay  rates  for  points  A  and  B  quickly  reach  a  common  constant  value. 

Low  pressure  response  shows  decay  rates  initially  increase  due  to  He  pressure 
before  attaining  the  pressure-broadened  asymptotic  constant  showing  that  He  helps 
the  energy  transfer  return  to  equilibrium  before  the  states  are  sufficiently  pressure 
broadened. 

It  is  therefore  beneficial  to  examine  the  effect  at  low  He  pressures  at  various 
temperatures.  Temperature  has  little  effect  on  Doppler  broadening  over  the  range 
used  in  this  experiment.  However  the  concentration  of  alkali-metal  vapor  atoms  varies 
exponentially  with  temperature  above  the  room  temperature  melting  point.  Collisions 
between  Cs-Cs  are  rare  at  93  °C  where  the  number  density  is  1.0  x  1012  atoms/cm3 
and  2.2  x  1013  atoms/cm3  by  105  °C.  Cesium-cesium  collisional  broadening  is  less 
then  natural  broadening  until  approximately  120  °C.  As  evidenced  by  the  change  in 
decay  rate  as  a  function  of  temperature  as  shown  in  Figure  36  temperature  effects 
the  data  fit  decay  rates  as  the  double-exponential  decay  becomes  a  single-exponential 
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decay.  This  suggests  that  decay  rates  trend  toward  a  single  dominant  hyperfine 
relaxation  process.  At  low  temperatures  at  least  two  processes  are  present  in  the 
decay  which  could  be  representative  of  Cs-wall  collisions  or  Cs-Cs  interactions. 


82 


VII.  Conclusion 


A  hyperfine  structure  cross  section  model  was  created  to  generate  simulated  ab¬ 
solute  absorption  coefficients  for  cesium  and  rubidium.  The  absorption  coefficients 
were  used  to  predict  group  velocity  and  group  delay  used  a  Kramers-Kronig  delay 
model  and  the  predicted  delays  were  validated  by  measuring  pulse  delays  across  the 
Cs  D\  and  Cs  D2  spectra. 

Tunable  slow  light  delays  are  demonstrated  through  frequency  modulation  of  a 
single  probe  laser  across  the  hyperfine  structure  of  an  alkali  metal  vapor.  Delays 
of  0  ns  to  37  ns  have  been  observed  in  a  single  pass  through  a  5  cm  cesium  cell  at 
137.2  °C  by  detuning  from  resonance  ±25  GHz.  Predicted  delays  closely  match  calcu¬ 
lations  from  an  absorption  profile  model  including  Voigt  lineshapes  and  full  hyperfine 
structure.  Frequency  dependence  of  the  delays  are  sensitive  to  vapor  pressure  and 
Doppler  broadening,  constraining  vapor  temperature  uncertainties  to  about  0.1  °C. 
Additionally  temporal  delay  of  pulses  using  a  gated  camera  demonstrated  the  pulse 
delay  effect  in  an  image. 

Tuning  of  optical  pulse  delays  obtained  in  the  hyperfine  doublet  structure  for 
the  alkali  atoms  may  be  achieved  through  frequency  modulation.  Indeed,  delays  of 
0  ns  to  37ns  have  been  observed  in  a  single  pass  through  a  5-crn  Cs  cell  at  137.2  °C 
by  detuning  from  resonance  25  GHz.  Delay  predictions  are  nearly  exact  when  the 
modeling  precisely  matches  the  absorption  profile  including  Voigt  lineshapes  and  full 
hyperfine  structure.  The  frequency  dependence  of  the  delays  are  sensitive  to  va¬ 
por  pressure,  constraining  cell  temperature  uncertainties  to  about  0.1  °C.  Using  the 
Voigt  profile  and  hyperfine  structure  allows  for  accurate  prediction  of  the  pulse  delay 
achieved  by  frequency  tuning. 

The  frequency  dependence  of  the  optical  pulse  delay  in  cesium  D\  closely  match 
the  previously  developed  model  predictions  using  Voigt  lineshapes  and  the  full  hyper- 
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fine  structure.  The  observed  D\  optical  pulse  delays  compared  favorably  to  similar 
results  in  Cs  D2.  The  delay  prediction  Kramers-Kronig  model  accurately  predicts 
the  delays  between  D\  upper  state  hyperfine  components  at  lower  temperatures  than 
previously  reported  with  smaller  frequency  sampling  separation.  Delays  up  to  24  ns 
where  achieved  by  detuning  from  resonance  by  10  GHz. 

Spectral  hole  burning  produced  optical  delays  in  Cs  D\  pulses  by  varying  D2  pump 
laser  power.  Increasing  D2  pump  power  reduced  D\  absorption  overall.  Observation 
of  hole  burning  effects  on  delay  close  to  the  spectral  hole  were  noted.  Location  of 
the  holes  and  predicted  delay  from  a  modified  absorption  spectrum  matched  model 
predictions.  For  a  pump  intensity  of  four  times  saturation  intensity  Isat  observed 
delay  is  reduced  by  an  average  of  50%  across  the  spectrum  and  reaches  a  maximum 
of  78%  in  a  narrow  110  MHz  spectral  hole  region. 

Hyperfine  relaxation  and  temporal  dynamics  were  examined  in  a  cesium  vapor  cell. 
Short-term  nanosecond  scale  pump  effects  and  long-term  microsecond  scale  relaxation 
of  hyperfine  populations  were  observed.  Temperature  and  helium  pressure-dependent 
absorption  features  were  noted. 

Fundamental  understanding  of  pulse  delays  may  lead  to  the  development  of  prac¬ 
tical  all-optical  delays  for  optical  communication  systems,  improved  interferometry 
instruments,  remote  sensing  techniques  and  analytic  devices.  Exploration  of  the  tem¬ 
poral  dynamics  in  alkali-metal  vapors  is  crucial  to  the  optimization  of  scalable  pumped 
DPAL  systems. 
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Appendix  A.  Physical  Constants 


Table  3.  Physical  constants^ 


Symbol 

Description 

Value/Units 

CXfs 

Fine-structure  constant 

1/137.035  999  679(94) 

c 

Speed  of  light  in  vacuum 

299  792  458  ms-1 

h 

Planck  constant 

6.626  068  96(33)  x  10“34  Js 
4.135  66733(10)  x  10”15eVs 

h 

Planck  constant  /  2n 

1.054  571628(53)  x  10~34Js 
6.582118  99(16)  x  10-16eVs 

ks 

Boltzmann  constant 

1.380  6504(24)  x  10”23  JK"1 
8.617343(15)  x  lO^eVK"1 

Co 

Permittivity  (Electric  constant) 

8.854187817  x  lO-^Fm"1 

ho 

Permeability  (Magnetic  constant) 

12.566  370  614  x  lO^NA"2 

e 

Elementary  charge 

1.602176487(40)  x  10~19C 

me 

Electron  mass 

9.109  38215(45)  x  10"31  kg 

rrip 

Proton  mass 

1.672  621637(83)  x  10~27kg 

mu 

Atomic  mass  constant 

1.660  538  782(83)  x  10~27kg 

atm 

Standard  atmosphere 

101  325  Pa  =  760  Torr 

t  NIST  2006  CODATA  recommended  values 
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Appendix  B.  Cesium  Data 


Physical  Properties 


Atomic  numbed 
Total  nucleons' 

Ground  stated 

Relative  natural  abundance* 
Atomic  mass' 

Melting  point* 

Boiling  point* 

Density  (solid)* 

Nuclear  spin* 

Magnetic  moment 
Ionization  energy* 


Z 

Z  +  N  133 

6s2Ai/2 

■q(133Cs)  100% 

m  132.905  452(1)  u 

Tm  28.44  °C 

Tb  671.0  °C 

p  1879  kg  m~3 

I  7/2 


/i  +2.579 

Ei  3.893  90(2)  eV 


Table  4.  Cesium  (Cs)  physical  properties 
Property  Symbol  Value 

55 


*  NIST  “Chemistry  Webbook.”  2005 

*  Winter,  M.  “WebElements.”  2008 


Cesium  Vapor  Pressure 

Cesium  vapor  pressure,  in  atmospheres,  and  Figure  42  are  derived  from  Equation  (104) 
(Alcock  et  al.,  1984)  where  temperature  T  (K)  is  between  melting  point  and  550  K. 
Number  density  derived  from  Equation  (104)  and  the  ideal  gas  las  is  shown  in  Figure  43. 

log10  Pv  =  8.232  -  4062T"1  -  1.3359  log10  T  (104) 


Cesium  D  Line  Optical  Properties 

D\  optical  properties  are  listed  in  Table  5  and  illustrated  in  Figure  44. 
D-2  optical  properties  are  listed  in  Table  6  and  illustrated  in  Figure  45. 
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Figure  42.  Cesium  vapor  pressure  as  given  by  equation  (104)  Alcock  et  al. 
(1984). 


Figure  43.  Cesium  number  density  derived  from  equation  (104)  Alcock 
et  al.  (1984). 
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F  =  4 


62Pi/2 
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Figure  44.  Cesium  D±  (62S1;/2  — >  62P!/2)  hyperflne  energy  level  split¬ 
ting.  Hyperflne  splitting  not  drawn  to  scale.  Relative  splitting  is  only 
accurate  for  each  hyperflne  level. 


F  = 


62P3/2 


852.34727582(27)  nm 
351.72571850(11)  THz 
11732.3071049(37)  cm1 
1.454620563(53)  eV 


62Si/2 


263.8906(24)  MHz 


12.78951(82)  MHz 


251.0916(20)  MHz 


188.4885(13)  MHz 


T 


201.2871(11)  MHz 


339.7128(39)  MHz  151.2247(16)  MHz 

i  A 


4.021  776  399  375  GHz  (exact) 


9.192  631  770  GHz 
(exact) 


5.170  855  370  625  GHz  (exact) 


F  =  4 

F  =  3 

F  =  2 

F  =  4 


F  =  3 


Figure  45.  Cesium  £)2  (62S!/2  — >  62P3/2)  hyperflne  energy  level  splitting. 
Hyperflne  splitting  not  drawn  to  scale.  Relative  splitting  is  only  accurate 
for  each  hyperflne  level. 


Table  5.  Cesium  D\  (62S!/2  — >  62Pi/2)  transition  optical  properties(Steck, 
2009) 


Property 

Symbol 

Value 

Frequency 

D) 

335.116  048  807(41)  THz 

Energy 

hv o 

1.385  928  495(34)  eV 

Wavelength 

A 

894.592  959  86(10)  nm 

Lifetime 

T 

34.791(90)  ns 

Decay  Rate 

r 

28.743(75)  x  HRs”1 

Natural  Line  Width  (FWHM) 

4.575(12)  MHz 

Absorption  Oscillator  Strength 

/ 

0.3449(26) 

Table  6.  Cesium  Z)2  (62S!/2  — >  62P3/2)  transition  optical  properties (Steck, 
2009) 


Property 

Symbol 

Value 

Frequency 

D) 

351.725  718  50(11)  THz 

Energy 

hiA) 

1.454  620  563(35)  eV 

Wavelength 

A 

852.347275  82(27)  nm 

Lifetime 

T 

30.405(77)  ns 

Decay  Rate 

r 

32.889(84)  x  lOS"1 

Natural  Line  Width  (FWHM) 

5.234(13)  MHz 

Absorption  Oscillator  Strength 

/ 

0.7164(25) 

Cesium  D  Line  Hyperfine  Energy  Levels 


A Ehfs 


1  4  TS  ,  D  §A’(A'+l)-2/(J+l)J(J+l) 
T  -D/i/s  47(27— 1)J(2J—1) 


,  r  5772(7\/4+l)+7C[7(7+l)+J(J+l)+3— 37(7+1) J(J+1)]  — 57(7+1) J(J+1) 
0/l/s  7(7— 1)(27— 1)  J(J— 1)(2J— 1) 


(105) 


where 


K  =  F(F  +  1)—  /(/  +  1)  -  J(J  +  1)  (106) 

Ahfs  =  magnetic  dipole  constant 
Bhfs  =  electric  dipole  constant 
Chfs  =  magnetic  octupole  constant 


Table  7.  Cesium  D  line  magnetic  and  electric  interaction  con- 


stants(Steck,  2009). 

62S1/2 

62P 1/2 

62P  3/2 

Ahfs  h-  2.298  157942  5  GHz 

h  •  291.9201(75)  MHz 

h  ■  50.288  27(23)  MHz 

Bhfs 

— 

h  ■  —0.4934(17)  MHz 

Chfs 

— 

h  ■  0.560(70)  kHz 

Cesium  D  Line  Strength 


Line  strengths  are  determined  using  the  Sp"  f>  from  (Steck,  2009,  eqn  (41),  Tab 


8). 


S, 


F' 


(2F'  +  1)(2J+1) 


F" 


(107) 


5V»  f'  obeys  the  rule(Steck,  2009,  eqn  (42)) 


^  Sf"  f'  =  1  (108) 

F' 


Strength  factor  values  can  be  determined  using  Mathematica's  SixJSymbol  func¬ 
tion  and  Equation  (107). 

Tables  8  and  9  catalog  the  transition  line  strength  values  for  the  cesium  D\  and 
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D2  hyperfine  lines. 


Table  8.  Cs  Di  (62S1/2  — >  62P!/2)  hyperflne  transition  strength  fac- 
tors(Steck,  2009). 


Transition 

Strength 

S F"  F' 

(unitless) 

544 

5/12 

S4  3 

7/12 

S3  4 

3/4 

S3  3 

1/4 

Table  9.  Cs  D2  (62S1/2  — »  62P3/2)  hyperflne  transition  strength  fac- 
tors(Steck,  2009). 


Transition  Strength 

Sp"  f'  (unitless) 


54  5 

11/18 

5*4  4 

7/24 

5*4  3 

7/72 

S3  4 

15/56 

S3  3 

3/8 

S3  2 

5/14 

Buffer  Gas  Pressure  Effects 

Pressure  effects  vary  by  buffer  gas  type  and  pressure.  The  line  is  broadened  as 
well  as  shifted  in  frequency.  Table  10  shows  the  frequency  shift  and  line  width  effects 
on  the  cesium  Di  and  D2  lines  due  to  nitrogen  and  helium  buffer  gases. 

Pressure  effects  must  also  be  corrected  for  temperature.  Given  the  data  was  taken 
at  Tdat.a  =  294.15  K,  the  buffer  gas  pressure  effects  on  system  at  temperature  Tsys,  in 
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Table  10.  Pressure  broadening  and  shift  of  cesium  D  line  due  to  buffer 
gases(Andalkar  and  Warrington,  2002).  Data  was  taken  at  294.15 K 


Buffer 

Gas 

Measured 

Quantity 

Value 

(MHz  Torr-1) 

n2 

D\  width 

19.51 

n2 

Di  shift 

-8.23 

n2 

D2  width 

22.68 

n2 

D2  shift 

-6.73 

He 

D i  width 

26.21 

He 

D\  shift 

4.46 

He 

D2  width 

23.50 

He 

D2  shift 

0.75 

K,  and  buffer  gas  pressure  Pt.  in  Torr  by 


A  UL  =  ^~  - 


2n  \  Tr 


T. 


sys 


T, 


data 


1/2 


Pi 


rri  \  1/2 

u0  S >  H)  +  ^  Si  f  S  'JS  j 

'  4  data  / 


(109) 

(110) 


where  is  the  Lorentzian  line  width,  and  u0  is  the  center  frequency  for  each 
hyperhne  line. 
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Appendix  C.  Matlab®  Code 


Cesium  D2  Absorption  Coefficient 


The  calculated  absorption  coefficient  compares  favorably  to  experimentally  deter¬ 
mined  Cs  absorption  as  shown  in  Figure  12  in  Chapter  IV  Delay  Prediction. 


function  [abs]  =  fnAbsCoef fCsD2 (freq, T, PN2, PHe, PCH4) 

o. 

O - 


%  Monte  D.  Anderson 

%  Department  of  Engineering  Physics,  Air  Force  Institute  of  Technology 
%  19  Mar  2010 

o, 

o _ 

%  Calculates  cesium  D2  absorption  coefficent  (1/cm)  at  a  given  frequency 
%  at  temperature  T  (K)  and  buffer  gas  pressures  PN2,  PHe,  and  PCH4  (Torr) 

o, 

o _ 


%%  CONSTANTS 

global  c  eo  h  hbar  kB  AMU; 
global  AmCs  MCs  Inam  vo  tau; 
global  TransM  numTrans; 
global  dvS12  sO  vD2; 


%  Physical  Constants  (NIST  2006  CODATA  Internationally 
%  http : / / physics . nist . gov/cuu/ Constants/ index . html 


c  =  2.997  92458e8 ; 
eo  =  8 . 85418781762e-12; 
h  =  6. 62606896e-34; 
hbar  =  1 . 054571628e-34; 
kB  =  1 ,3806504e-23; 

AMU  =  1 . 660538782e-27; 


vacuum  speed  of  light  (m/ 
vacuum  permittivity  (F/m) 
Planck  const  (Js) 

Planck  const/2pi  (Js) 
Boltzmann  const  (J/K) 
Atomic  Mass  Unit  (kg) 


recommended  values) 
s) 


%  Cesium  Atomic  Data 
AmCs  =  132.905451931; 
MCs  =  AmCs * AMU; 

Inam  =  7/2; 


%  Atomic  mass  (AMU) 

%  Mass  of  Cs  atom  (kg) 

%  Total  nuclear  angular  momentum 


%  Cs  D2  Data 

vo  =  351 . 72571850el2; 

tau  =  30.32e-9; 


%  Frequency  (Hz) 
%  Lifetime  (s) 


REF  [1] 


REFERENCE : 

[1]  J.  M.  Amini  and  H.  Gould 

High  Precision  Measurement  of  the  Static  Dipole  Polarizability  of  Cs 
Phys .  Rev.  Lett.,  American  Physical  Society,  2003,  91,  153001 


%  Transition  Matrix,  contains  allowed  transitions  and  pertenent  quantum 
%  numbers  in  the  format  (Inam,  Jupper,  Flower,  Fupper) 
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59 
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65 

66 

67 

68 
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70 

71 

72 

73 

74 
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%  NOTE:  Jlower  is  1/2  (ground  state) 

%  NOTE:  Inam  is  the  same  for  all  transitions  in  Cs .  This  function  was 
%  built  to  accomdate  elements  with  multiple  isotopes  like  Rubidium. 
TransM  =  [7/2  3/2  3  2; 

7/2  3/2  3  3; 

7/2  3/2  3  4; 

7/2  3/2  4  3; 

7/2  3/2  4  4; 

7/2  3/2  4  5] ; 

numTrans  =  length (TransM) ; 

%%  Hyperfine  Structure  Energy  Levels 
vD2  =  zeros ( 1 , numTrans ) ; 
for  j  =  1: numTrans 

vD2 ( j )  =  (hf sEnergy (TransM ( j , 1 ) , 1, TransM ( j ,  2 )  ,  TransM (j,4))/h)-... 

(hf sEnergy (TransM (j,l), 0,1/2, TransM ( j , 3) ) /h) ; 

end 

dvS12  =  hf sEnergy ( Inam, 0 , 1/2 , 4 ) /h  -  hf sEnergy ( Inam, 0 , 1 /2 , 3 ) /h; 

%%  Line  Strength  Factors,  SFF 
SFF  =  zeros ( 1 , numTrans ) ; 
for  n  =  1: numTrans 

SFF (n)  =  S_f un (TransM (n, 3 ) , TransM (n, 4 ) , 0 . 5, TransM (n, 2 ) , TransM (n, 1 )  )  ; 

end 

%%  Cross  Section  at  Line  Center  (sigma_0) 
sO  =  zeros ( 1 , numTrans ) ; 
for  f  =  1: numTrans 

sO(f)  =  sigma_0 (SFF,  vD2, f ) ; 

end 

%%  main  function  output 

abs  =  NumDensity (T) *CrossSection ( f req, T, PN2 , PHe, PCH4 ) ; 

end  %end  of  fnAbsCoef fCsD2  main  function  - 


o, 

o _ 

%%  f nAbsCoef f CsD2  SUBFUNCTIONS 

%%  HYPERFINE  ENERGY  LEVELS  - 

function  [dEhfs]  =  hf sEnergy ( Inam, L, J, F) 

%  Calculates  hyperfine  structure  energy 

o, 

o _ 

%  REFERENCES: 

%  [2]  E.Arimondo,  M.Inguscio,  and  P.Violino 

%  "Experimental  determinations  of  the  hyperfine  structure... 

%  ...in  the  alkali  atoms" 

%  Rev.  Mod.  Phys . ,  American  Physical  Society,  1977,  49,  31-75 

o. 

o 

%  [3]  V.Gerginov,  A. Derevianko,  and  C.E. Tanner 

%  "Observation  of  the  Nuclear  Magnetic  Octupole  Moment  of  Csl33" 


94 


94 

95 

96 

97 

98 

99 

100 

101 

102 

103 

104 

105 

106 

107 

108 

109 

110 

111 

112 

113 

114 

115 

116 

117 

118 

119 

120 
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123 

124 
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Phys .  Rev.  Lett.,  American  Physical  Society,  2003,  91,  072501 


o, 

o _ 

global  h; 


%Constants 

AS 12  =  h*2 ,2981579425*1E9; 

API 2  =  h*291 . 9201*1E6; 

AP32  =  h*50.28827*lE6; 

BP 3 2  =  h*-0 . 4934*1E6; 

CP32  =  h*0 . 560*1E3 ; 

K  =  F*  (F  +  1)  -  Inam* (Inam  +  1) 

if (L  ==  0) 

dEhfs  =  . 5*AS12  *K; 
elseif (L  ==  1  &&  J  ==  .5) 
dEhfs  =  . 5*AP12*K; 
elseif (L  ==  1  &&  J  ==  1.5) 

dEhfs  =  ( . 5*AP32*K)  +  BP32* ( (1 . 5*K* (K+l) 

2*Inam* (Inam+1) *J* (J+l) ) / (4*Inam* (2*Inam-l) *J* (2*J-1) ) )  +. . . 

CP 32* (5* (K~2) * ( (K/4 ) +1 )  +... 

K* ( Inam* ( Inam+1 )  +  J*(J+1)  +  3  -  3*Inam* (Inam+1) *J*  (J+l) )  -... 
5*Inam* (Inam+1) *J* (J+l) ) / (Inam* (Inam-1) * (2*Inam-l) *J* ( J-l) * (2*J-1) ) ; 

else 

fprintf ( ' \nIMPROPER  INPUTS\n\n') 
dEhfs  =  0; 

end 


%  (Hz) 

%  (Hz) 

%  (Hz) 

%  (Hz) 

%  (Hz) 

-  J* ( J  +  1)  ; 


end 


%%  LINE  STRENGTH  FACTORS  - 

function  [S]  =  S_fun (F, Fprime, J, Jprime, I) 

SI  =  Wigner 6 j (J, Jprime, 1 , Fprime,  F,  I )  ; 

S  =  (2*Fprime  +  l) * (2*J+1) * (Sl~2)  ; 
end 

%%  WIGNER  6 j  FACTORS - 

function  Wigner  =  Wigner6 j ( jl,  j2,  j3,  Jl,  J2,  J3) 

%  Wigner  6j-symbol  calculator.  Written  by  Amita  B  Deb,  Clarendon  Lab.  2007. 
%  Improved  by  Richard  A.  Holt,  Univ.  of  Western  Ontario,  2009. 

%  Further  improved  to  deal  with  large  input  arguments  by  Lee  Harper, 

%  CRL  Oxford,  2009. 

%  Builds  on  the  m-file  Wigner6 jcoef f  posted  by  Richard  A.  Holt. 

%  Calculates  {  jl,  j2  , j 3 }  using  Racah  formula. 

%  Jl  J2  J3 

o. 

O - 

%  REFERENCE: 

%  [4]  I.I.Sobelman 

%  Atomic  Spectra  and  Radiative  Transitions 

%  Springer  Series  in  Chemical  Physics,  Vol .  1,  Springer  (June  1979) 
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o. 

o _ 

%  Finding  Triangular  coefficients 
tril  =  triangle_coef f ( jl, j2, j3) ; 
tri2  =  triangle.coef f ( jl, J2, J3) ; 
tri3  =  triangle.coef  f  ( Jl,  j2,  J3)  ; 
tri4  =  triangle.coef f ( Jl, J2, j3) ; 

if  (tril  =  =  0  1 1  tri2  =  =  0  1 1  tri3  =  =  0  1 1  tri4==0) 

Wigner=0 ; 
return 

end 

%  Finding  the  range  of  summation  in  the  Racah  formula. 


a  ( 1 )  = 

jl 

+ 

j2 

+ 

j3; 

a  (2)  = 

jl 

+ 

J2 

+ 

J3 ; 

a  (3)  = 

Jl 

+ 

j2 

+ 

J3 ; 

a  (4)  = 

Jl 

+ 

J2 

+ 

j3; 

rangei 

max (a) 

r 

Ml)  = 

jl 

+ 

j2 

+ 

Jl 

+ 

J2 ; 

M2)  = 

j2 

+ 

j  3 

+ 

J2 

+ 

J3 ; 

M3)  = 

j3 

+ 

jl 

+ 

J3 

+ 

Jl; 

rangef 

min (k) 

r 

Wigner  =  0; 

for  t=rangei : rangef 

Wigner  =  Wigner  +  ( (-1 ) "t ) *exp (gammaln (t+2 ) 
fung (t, jl, j2, j3, Jl, J2,  J3) ) ; 

end 

Wigner  =  (tril*tri2*tri3*tri4) " (0 . 5) *Wigner; 


end 

% - Wigner 6 j  subfunctions - 

function  r  =  fung(t,  j 1 ,  j2 ,  j 3 ,  Jl ,  J2 ,  J3 ) 

%  Calculating  the  logarithm  of  the  denominator  in  Racah  Formula,  using 
%  the  gamma  function  in  place  of  the  factorial. 

r  -  sum (gammaln ( [ (t-jl-j2-j3) ; (t- jl-J2-J3) ; (t-Jl- j2-J3) ; (t-Jl-J2- j3) ; . . . 

( jl+j2+Jl+J2-t) ; ( j2  + j3+J2+J3-t) ; ( j3+ jl+J3+Jl-t ) ]  +  1) ) ; 

end 

o. _ 

o 

function  tri  =  triangle_coef f (a,  b,  c) 

%  Calculates  triangle  coefficients  for  angular  momenta. 

%  This  version  returns  0  if  the  triangle  inequalities  are  violated.  (RAH) 
if  (a<0  |  |  b<0  |  |  c<0) 
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tri=0 ; 
return 

end 

for  xa  =  abs  (a-b)  : 1 : a+b 
if  c==xa 

tri  =  factorial (a+b-c) ^factorial (a-b+c) . 

factorial (-a+b+c) / (factorial (a+b+c+1) )  ; 
return 

end 

end 


tri=0 ; 


end 

%%  CROSS  SECTION  AT  LINE  CENTER  - 

function  [sO]  =  sigma_0 (SFF, vD,  x) 

%  Calculates  cross  section  at  line  center  for  each  hyperfine  level, 
global  c  vo  tau  TransM; 

dvnat  =  1/  (2*pi*tau) ;  %  natural  linewidth 

sO  =  SFF(x) * (1/tau) * ( ( (c*100) ~2) / (8*pi) ) * (1/ ( (vo+vD(x) ) ~2) ) *. . . 

( (2*TransM (x, 2) +1) / (2*. 5+1) ) * (2/ (pi * dvnat ) ) ; 

end 

%%  NUMBER  DENSITY  - 

function  [ndensity]  =  NumDensity (T) 

%  calculates  number  density  (cm'' -3)  given  Temperature  (K)  from  ideal  gas  eq 
global  kB; 

ndensity  =  (CsVaporPressure (T) )/ (kB*T) *101325* ( le-6)  ; 
end 

%%  VAPOR  PRESSURE  - 

function  [Pv]  =  CsVaporPressure (T) 

%  Provides  Cesium  vapor  pressure  in  atmospheres  given  Temperature  (K) 

o. 

o _ 

%  REFERENCE: 

%  [5]  C.B.Alcock,  V.P.Itkin,  and  M.K.Horrigan 

%  Vapour  pressure  equations  for  the  metallic  elements:  298-2500K 

%  Canadian  Metallurgical  Quarterly,  1984,  23,  309-313 

O, 

O _ _ 

%  Uses  "Precise"  equation  for  Cs  vapor  pressure. 

%  Uses  melting  point  for  Cs  at  301. 59K  (NIST  value) 

%  Valid  for  temperatures  from  melting  point  to  550K  (per  Alcock  et  al) . 

if (T  >  301.59  &  &  T  <  550) 
a  =  8.232; 
b  =  -4062; 
c  =  -1.3359; 
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elseif  (T  <  301.59  ||  T  >  550) 

fprint f ( ' \nTemperature  is  outside  of  allowed  range. \n\n') 
return 

end 

Pv  =  10.  ~  (a  +  b*(T~-l)  +  c*loglO (T) ) ; 
end 

%%  CROSS  SECTION  - 

function  [cs]  =  CrossSection (x, T, PN2, PHe, PCH4 ) 

%  calculates  cross  section  as  function  of  frequency,  x 
global  tau  sO  vo  TransM  vD2; 

dvnat  =  1/ (2*pi*tau)  ; 

[gamma, shift, dvD, dvL]  =  Buf ferGas (T, PN2 , PHe, PCH4 )  ; 

cs  =  (dvnat*pi/2 ) * ( . . . 

Bolt zRatio (TransM (l,3),T)*sO(l)*... 

( (x+vD2 (1) +shift+vo) / (vD2 (1) +shift+vo) ) * . . . 
gVoigt ( (x-vD2 ( 1 ) ) , dvL, dvD) + . . . 

Bolt zRatio (TransM (2,3),T)*s0(2)*... 

( (x+vD2 (2) +shift+vo) / (vD2 (2) +shift+vo) ) * . . . 
gVoigt  (  (x-vD2  (2 )  )  ,  dvL,  dvD)  +  .  .  . 

Bolt zRatio (TransM (3,3),T)*sO(3)*... 

( (x+vD2 (3) +shift+vo) / (vD2 (3) +shift+vo) ) * . . . 
gVoigt  (  (x-vD2  (3)  )  ,  dvL,  dvD)  + .  .  . 

Bolt zRatio (TransM (4,3),T)*sO(4)*... 

( (x+vD2 (4) +shift+vo) / (vD2 (4) +shift+vo) ) * . . . 
gVoigt  (  (x-vD2  (4  )  )  ,  dvL,  dvD)  +  .  .  . 

Bolt zRatio (TransM (5,3),T)*sO(5)*... 

( (x+vD2 (5) +shift+vo) / (vD2 (5) +shift+vo) ) * . . . 
gVoigt  (  (x-vD2  (5)  )  ,  dvL,  dvD)  +  .  .  . 

Bolt zRatio (TransM (6,3),T)*sO(6)*... 

( (x+vD2 ( 6) +shift+vo) / (vD2 ( 6 ) +shift+vo) ) * . . . 
gVoigt  (  (x-vD2  (6)  )  ,  dvL,  dvD)  )  ; 

end 


%%  BUFFER  GAS  EFFECTS  - 

function  [gamma, shift, dvD,  dvL]  =  Buf ferGas (T, PN2 , PHe, PCH4 ) 

%  Cs  D2  buffer  gas  effects,  including  broadening  and  frequency  shift 
%  OUTPUT:  gamma  pressure  broadening  due  to  buffer  gas 

%  shift  frequency  shift  due  to  buffer  gas 

%  dvD  Doppler  broadening 

%  dvL  Lorentzian  homogeneous  broadening  (natural  +  pressure) 

o, 

o _ 

%  REFERENCES: 

%  [6]  A.Andalkar  and  R . B . Warrington 

%  High-resolution  measurement  of  the  pressure  broadening  and  shift 

%  of  the  Cs  D1  and  D2  lines  by  N2  and  He  buffer  gases 

%  Phys .  Rev.  A  65,  032708  (2002) 


[7]  G.  Pitz,  D.  Wertepny,  and  G.  Perram 

Pressure  broadening  and  shift  of  the  cesium  D1  transition  by  the 
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%  noble  gases  and  N2,  H2,  HD,  D2,  CH4,  C2H6,  CF4,  and  3He 

%  Phys .  Rev.  A  80,  062718  (2009) 

o, 

O - 

%  Pressure  broadening  coefficients  are  temperature  dependent  such  that 
%  (gamma  or  shift) *( (T/Tdata) " (1/2) ) ^Pressure 
%  where 

%  gamma  or  shift  are  the  broadening  or  shift  values  (MHz/Torr) 

%  T  is  the  current  temperature 

%  Tdata  is  the  temperature  at  which  buffer  gas  values  were  measured 

o. 

o 

%  Pressure  effect  values  and  Tdata  below  are  taken  from  reference  [2] . 
global  tau  c  kB  MCs  vo; 

gamma  =  (19 . 18e6) * ( (T/313) ~ (1/2) ) *PN2  ... 

+  (20.59e6) * ( (T/313) ~ (1/2) ) *PHe. . . 

+  (25.84e6) *  (  (T/313) ~  (1/2) ) *PCH4; 
shift  =  (-6 . 2e6) * ( (T/313) " (1/2) ) *PN2. . . 

+  (0.69e6) * ( (T/313) ~ (1/2) ) *PHe. . . 

+  (-8 . 8  6e6) *  (  (T/313) ~  (1/2) ) *PCH4; 

dvD  =  (vo  +  shift) * ( (8*kB*T*log (2) ) / (MCs* (c"2) ) ) ~ (1/2) ; 

dvL  =  (1/ (2*pi*tau) )  +  gamma; 

end 


BOLTZMAN  RATIO 


function  [BR]  =  BoltzRatio (F, T) 

%  Calculates  relative  Boltzman  distribution  of  ground  state  energy  levels, 
global  kB  h  dvS12; 

if (F  ==  3) 

BR  =  7 . / ( 7  +  9*exp (- (h*dvS12 ) . / (kB . *T) ) ) ; 
elseif (F  ==  4 ) 

BR  =  ( 9* exp (- (h*dvS12 ) . / (kB . *T) ) ) . / ( 7  +  9*exp (- (h*dvS12 ) . / (kB . *T) ) ) ; 
elseif (F  ^  3  &&  F  ^  4) 

fprintf ( ' \nIMPROPER  INPUTS\n\n') 

BR  =  0; 

end 

end 


%%  VOIGT  LINESHAPE  - 

function  [y]  =  gVoigt (v, dvL, dvD) 

%  A  useful  approximation  for  the  Voigt  convolution  is  given  by  Thompson 

o, 

o _ 

%  REFERENCE: 

%  [8]  William  J.  Thompson 

%  Numerous  neat  algorithms  for  the  Voigt  profile  function. 

%  Computers  in  Physics,  Vol  7,  No.  6,  Nov/Dec  1993. 

o, 

o _ 

%  gV(u,a)  =  Re [Exp ( z ~ 2 ) erf c ( z ) ] 

%  z  =  a  +  iu 

%  a  =  sqrt(ln(2))  vLorentzian/vDoppler 
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u  =  2  sqrt(ln(2))  (v-vo) /vDoppler 

NOTE:  Not  valid  for  Doppler  free  spectra  where  vDoppler  =  0 
NOTE:  erfc(z)  =  1  -  erf(z)  (related  to  the  Faddeeva  function) 

a  =  sqrt (log (2) ) * (dvL/dvD) ; 
u  =  2*sqrt (log (2) ) * (v/dvD)  ; 
z  =  a  +  li*u; 

y  =  2*sqrt (log (2) /pi) * (1/dvD) *real (faddeeva (li*z) ) ; 
end 


gVoigt  subfunctions 


function  w  =  f addeeva ( z ,  N) 

o, _ 

o 


o. 

o 

o, 

o 

o. 

o 

o. _ 

o 

o. 

o 

o. 

o 

o. 

o 

o_ 

o 

o, 

o 

o. 

o 

o. 

o 

o. 

o 

o. 

o 

o. 

o 

o. 

o 

o. 

o 

o. 

o 

o_ 

o 

o_ 

o 

o, 

o 

o_ 

o 

o. 

o 

o. 

o 

o. 

o  _ 

o, 

o 

o, 

o 

o, 

o 

o, 

o 

Q. 

o  _ 

if 

if 


http: //www.mathworks . com/matlabcentral/ f ileexchange/ . . . 

. . . 22207-faddeeva-function-f ft-based 
Posted  by  Kesh  (tikuma@lsu.edu)  -  21  Nov  2008  (Updated  19  Feb  2009) 


FADDEEVA  Faddeeva  function 

W  =  FADDEEVA (Z)  is  the  Faddeeva  function,  aka  the  plasma  dispersion 
function,  for  each  element  of  Z.  The  Faddeeva  function  is  defined  as: 

w(z)  =  exp(-z~2)  *  erfc(-j*z) 

where  erfc(x)  is  the  complex  complementary  error  function. 

W  =  FADDEEVA (Z, N)  can  be  used  to  explicitly  specify  the  number  of  terms 
to  truncate  the  expansion  (see  (13)  in  [1]) .  N  =  16  is  used  as  default. 

Example : 

x  =  linspace (-10, 10,  1001)  ;  [X,Y]  =  meshgrid (x, x) ; 

W  =  faddeeva (complex (X, Y) ) ; 
figure ; 

subplot ( 12 1 ) ;  imagesc (x, x,  real (W) )  ;  axis  xy  square;  caxis ( [-1  1]); 
title ( 're ( faddeeva ( z ) )  ' )  ;  x label ('re(z)  ' ) ;  y label ( ' im ( z )  ' ) ; 
subplot ( 122 )  ;  imagesc (x, x, imag (W) ) ;  axis  xy  square;  caxis ([-1  1]); 
title ( ' im (faddeeva (z) ) ' ) ;  xlabel ( ' re (z) ' ) ;  ylabel ( ' im (z) ' ) ; 


Reference : 

[9]  J.A.C.  Weideman,  "Computation  of  the  Complex  Error  Function,"  SIAM 
J.  Numerical  Analysis,  pp.  1497-1518,  No.  5,  Vol.  31,  Oct.,  1994 
Available  Online:  http://www.jstor.org/stable/2158232 


nargin<2,  N  =  [];  end 
isempty(N),  N  =  16;  end 


w  =  zeros (size (z) ) ;  %  initialize  output 


o,  o,  o.  o_  o, 
o  o  o  o  o 


%  for  purely  imaginary-valued  inputs, 
idx  =  real(z)==0;  % 


use  erf  as  is  if  z  is  real 
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w(idx)  =  exp (-z (idx) . ~2) . *erfc (imag (z (idx) ) ) ; 


if  all (idx),  return;  end 
idx  =  — 'idx; 

o.  o,  o.  o,  o, 
o  o  o  o  o 

%  for  complex-valued  inputs 

%  make  sure  all  points  are  in  the  upper  half-plane  (positive  imag.  values) 
idxl  =  idx  &  imag(z)<0; 
z (idxl)  =  con j (z (idxl) ) ; 

M  =  2  *N; 

M2  =  2  *M; 

k  =  (-M+1:1:M-1) %  M2  =  no.  of  sampling  points. 

L  =  sqrt (N/sqrt ( 2 ) ) ;  %  Optimal  choice  of  L. 

theta  =  k*pi/M; 

t  =  L*tan (theta/2 ) ;  %  Variables  theta  and  t. 
f  =  exp (-t . "2) . * (L'2+t . ~2) ; 
f  =  [0;  f ] ;  %  Function  to  be  transformed. 

a  =  real (fft (fftshift (f ) ) ) /M2;  %  Coefficients  of  transform, 
a  =  flipud (a (2 :N+1) ) ;  %  Reorder  coefficients. 

Z  =  (L+li*z (idx) ) . / (L-li*z (idx) ) ; 
p  =  polyval (a, Z) ;  %  Polynomial  evaluation. 

w(idx)  =  2*p . / (L-li*z (idx) )  . ~2  +  (1/sqrt (pi) )  . / (L-li*z (idx) )  ;  %  Eval  w(z). 

%  convert  the  upper  half-plane  results  to  the  lower  half-plane  if  necesary 
w (idxl)  =  con j (2*exp (-z (idxl) . ~2)  -  w (idxl) ) ; 

end 
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Absorption  Coefficient  Changes  for  Cs  D1 


The  previous  Matlab®  code  example  can  be  easily  modified  for  cesium  Dl  or 
any  other  atomic  species  with  a  single  isotope.  To  change  the  Matlab®  absorption 
coefficient  function  on  page  C  for  cesium  Di,  make  the  following  modifications  at 
the  appropriate  places  within  the  code.  To  modify  use  of  this  code  for  other  atomic 
elements,  be  sure  to  modify  vapor  pressure  and  buffer  gas  effect  functions. 


%  Cs  DI  Data 

vo  =  335 . 116048807el2;  %  Frequency  (Hz) 

tau  =  34.72e-9;  %  Lifetime  (s) 


TransM  =  [7/2  1/2  3  3; 

7/2  1/2  3  4; 

7/2  1/2  4  3; 

7/2  1/2  4  4]; 

%hfsEnergy  Constants... 

AS 12  =  h*2 ,2981579425*1E9;  %(Hz) 

API 2  =  h*291 . 9201*1E6;  %(Hz) 

AP32  =  h*50.28827*lE6;  %(Hz) 

BP32  =  h*-0 . 4934*1E6;  %(Hz) 

CP32  =  h*0 . 560*1E3;  %(Hz) 

%fnCross Sect ion . 

cs  =  (dvnat*pi/2) *  (  .  .  . 


Bolt zRatio (TransM (l,3),T)*sO(l).*... 

( (x+vDl (1) +shift+vo) / (vDl (1) +shift+vo) ).*... 
gVoigt ( (x-vDl ( 1 ) ) , dvL, dvD) + . . . 

Bolt zRatio (TransM (2,3),T)*sO(2).*... 

( (x+vDl (2) +shift+vo) / (vDl (2) +shift+vo) ).*... 
gVoigt ( (x-vDl (2 ) ) , dvL, dvD) + . . . 

Bolt zRatio (TransM (3,3),T)*sO(3).*... 

( (x+vDl (3) +shift+vo) / (vDl (3) +shift+vo) ).*... 
gVoigt ( (x-vDl (3) ) , dvL, dvD) + . . . 

Bolt zRatio (TransM (4,3),T)*sO(4).*... 

( (x+vDl ( 4 ) +shift+vo) / (vDl ( 4 ) +shift+vo) ).*... 
gVoigt ( (x-vDl (4 ) ) , dvL, dvD) ) ; 

%Buf ferGasEf fects .  .  . 

gamma  =  (15. 82e6) * ( (T/318) ~ (1/2) ) *PN2  ... 

+  (24 . 13e6) *  (  (T/323) ~  (1/2) ) *PHe.  .  . 

+  (29.00e6) * ( (T/333) ~ (1/2) ) *PCH4; 

shift  =  (-7. 69e6)  *  (  (T/318C  (1/2)  )  *PN2.  .  . 

+  ( 4 . 24e6 ) *  (  (T/323) ~ (1/2) ) *PHe.  .  . 

+  (-9 . 28e6 ) * ( (T/333) ~ (1/2) ) *PCH4; 
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Cesium  D2  Hyperfine  Energy  Level 


The  following  code  sample  is  a  code  snippet  that  calculates  the  cesium  Do  hyper¬ 
fine  energy. 


function  [dEhfs]  =  hf sEnergy ( Inam, L, J, F) 

%  Calculates  hyperfine  structure  energy  for  Cs  D2.  Given  Inam,  L,  J  and  F 
%  calculates  the  hyperfine  energy  difference. 

o, 

o 

%  REFERENCES: 

%  Arimondo,  E.,  Inguscio,  M.  and  Violino,  P. 

%  Experimental  determinations  of  the  hyperfine  structure  in  the  alkali  atoms 
%  Rev.  Mod.  Phys . ,  American  Physical  Society,  1977,  49,  31-75 

o. 

o 

%  Gerginov,  V.;  Derevianko,  A.  &  Tanner,  C.  E. 

%  Observation  of  the  Nuclear  Magnetic  Octupole  Moment  of  Csl33 
%  Phys.  Rev.  Lett.,  American  Physical  Society,  2003,  91,  072501 

o, 

o 

global  h; 


%Constants 

AS 12  =  h*2 ,2981579425*1E9;  %(Hz) 
API 2  =  h*291 . 9201*1E6;  %(Hz) 
AP32  =  h*50.28827*lE6;  %(Hz) 
BP32  =  h*-0 . 4934*1E6;  %(Hz) 
CP32  =  h*0 . 560*1E3 ;  %(Hz) 


K  =  F*  (F  +  1)  -  Inam* (Inam  +1)  -  J*(J  +  1); 

if (L  ==  0) 

dEhfs  =  . 5*AS1 2  *K; 
elseif (L  ==  1  &&  J  ==  .5) 
dEhfs  =  . 5*AP 12  *K; 
elseif (L  ==  1  &&  J  ==  1.5) 

dEhfs  =  ( . 5*AP32*K)  +  BP32* ( (1 . 5*K* (K+l)  -  ... 

2*Inam* (Inam+1) *J* (J+l) ) / (4*Inam* (2*Inam-l) *J* (2*J-1) ) )  +. . . 

CP 32* (5* (K~2) * ( (K/4 ) +1 )  +  K* ( Inam* ( Inam+1 )  +  J*(J+1)  +  3  -  ... 
3*Inam* (Inam+1) *J* (J+l) )  -. . . 

5*Inam* (Inam+1) *J* (J+l) ) / (Inam* (Inam-1) * (2*Inam-l) *J* ( J-l) * (2*J-1) ) ; 

else 

fprintf ( ' \nIMPROPER  INPUTS\n\n') 
dEhfs  =  0; 

end 


%end  of  function 
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Rubidium  D\  Absorption  Coefficient 


The  absolute  absorption  of  the  rubidium  D  lines  compares  favorably  to  the  theo¬ 
retical  and  experimental  results  reported  by  Siddons  et  a/.(Siddons  et  ah,  2008). 


function  [abs]  =  fnAbsCoef fRbDl (freq, T, PHe, PCH4 ) 

%  Calculates  rubidum  D1  absorption  coefficent  (1/cm)  at  a  given  frequency 
%  at  temperature  T  (K)  and  buffer  gas  pressures  PHe  and  PCH4  (Torr) 

o, _ 

o 

%  Monte  Anderson 

%  Air  Force  Institute  of  Technology,  Department  of  Engineering  Physics 

o, 

o 

%  27  Oct  2009 


%%  CONSTANTS 

global  c  eo  h  hbar  kB  AMU; 
global  AmRb85  MRb85  Inam85  vo85  tau; 
global  AmRb87  MRb87  Inam87  vo87; 
global  TransM  numTrans; 
global  dv85S12  dv87S12  sO  vHF; 


%  Physical  Constants  (NIST) 
c  =  2.997  92458e8 ; 
eo  =  8 . 85418781762e-12; 
h  =  6 . 62  60  68  96e-34 ; 
hbar  =  1 . 054571628e-34; 
kB  =  1 . 3806504e-23; 

AMU  =  1 . 660538782e-27; 


%  vacuum  speed  of  light  (m/s) 
%  vacuum  permittivity  (F/m) 

%  Planck  const  (Js) 

%  Planck  const/2pi  (Js) 

%  Boltzmann  const  (J/K) 

%  Atomic  Mass  Unit  (kg) 


%  Rubidium  Atomic  Data 
AmRb85  =  84.911789732; 
AmRb87  =  86.909180520; 
MRb85  =  AmRb85*AMU; 
MRb8  7  =  AmRb87*AMU; 
Inam85  =  5/2; 

Inam87  =  3/2; 


%  Atomic  mass  (AMU)  of  Rb85 
%  Atomic  mass  (AMU)  of  Rb87 
%  Mass  of  Rb85  atom  (kg) 

%  Mass  of  Rb87  atom  (kg) 

%  Total  nuclear  angular  momentum  Rb85 
%  Total  nuclear  angular  momentum  Rb87 


%  Rb  D1  Data 

vo85  =  377 . 107385690el2; 
vo87  =  377 . 1074633805el2; 
tau  =  27 . 679e-9; 


%  Frequency  (Hz)  Rb85 
%  Frequency  (Hz)  Rb87 
%  Lifetime  (s)  Rb85  and  Rb87 


%  REFERENCES 


%  Transition  Matrix,  contains  allowed  transitions  and  pertenent  quantum 
%  numbers  in  the  format  (Inam,  Jupper,  Flower,  Fupper) 
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%  NOTE:  Jlower 

'  is 

1/2  (ground  state) 

TransM  =  [5/2 

1/2 

2 

2; 

5/2 

1/2 

2 

3; 

5/2 

1/2 

3 

2; 

5/2 

1/2 

3 

3; 

3/2 

1/2 

1 

1; 

3/2 

1/2 

1 

2; 

3/2 

1/2 

2 

1; 

3/2 

1/2 

2 

2] ; 

numTrans  =  length (TransM) ; 

%%  Hyperfine  Structure  -  Transition  Frequencies  (from  TransM) 
vHF  =  zeros (l,numTrans) ; 
for  j  =  1 : numTrans 

vHF ( j )  =  (hf sEnergyRb (TransM ( j , 1 ) , 1 , TransM ( j ,  2 )  ,  TransM  (j,4))/h)-... 

(hf sEnergyRb (TransM (j, 1) ,  0, 1/2, TransM ( j , 3 ) ) /h)  +  dvlso (TransM ( j , 1 ) ) ; 

end 

dv85S12  =  hf sEnergyRb ( Inam85 , 0 , 1 /2 , 3 ) /h  -  hf sEnergyRb ( Inam85 , 0 , 1 /2 , 2 ) /h; 
dv87S12  =  hfsEnergyRb (Inam87, 0, 1/2, 2) /h  -  hf sEnergyRb ( Inam87 , 0 , 1 /2 , 1 ) /h; 

%%  Line  Strength  Factors,  SFF 
SFF  =  zeros ( 1 , numTrans ) ; 
for  n  =  1: numTrans 

SFF (n)  =  S_f un (TransM (n, 3 ) , TransM (n, 4 ) , 0 . 5, TransM (n, 2 ) , TransM (n, 1 ) ) ; 

end 

%%  Cross  Section  at  Line  Center  (sigma_0) 
sO  =  zeros ( 1 , numTrans ) ; 
for  f  =  1: numTrans 

s0(f)  =  sigma_0Rb (SFF, vHF, f ) ; 

end 

%%  main  function  output 

abs  =  NumDensityRb (T) *CrossSect ionRb ( f req,  T,  PHe,  PCH4 )  ; 

end  %end  of  f nAbsCoef f RbDl  main  function - 


%%  f nAbsCoef f RbDl  SUBFUNCTIONS 


%%  HYPERFINE  ENERGY  LEVELS  - 

function  [dEhfs]  =  hfsEnergyRb ( Inam, L, J, F) 

%  Calculates  hyperfine  structure  energy  for  Rb  D1  and  D2 . 

o, 

o 

%  REFERENCES: 

%  Arimondo,  E.,  Inguscio,  M.  and  Violino,  P. 

%  Experimental  determinations  of  the  hyperfine  structure  inthe  alkali  atoms 
%  Rev.  Mod.  Phys . ,  American  Physical  Society,  1977,  49,  31-75 

o, 

o 

global  h  Inam85  Inam87; 
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%Constants  Rb85 

A85S12 

=  h*l . 0119108130e9; 

%  (Hz 

A85P12 

=  h*120.527e6; 

%  (Hz 

A85P32 

=  h*25.0020e6; 

%  (Hz 

B85P32 

=  h*25 . 790e6; 

%  (Hz 

%Constants  Rb87 

A87S12 

=  h*3 . 417341305452145e9; 

%  (Hz 

A87P12 

=  h*407 . 24e6; 

%  (Hz 

A87P32 

=  h  *  8  4 . 7 1 8  5  e  6 ; 

%  (Hz 

B87P32 

=  h  *  1 2 . 4965e6; 

%  (Hz 

K  =  F* (F  +  1 )  -  Inam* (Inam  +1)  -  J* ( J  +  1) ; 

if  (Inam  ==  Inam85)  %  Rb85  hyperfine  energy  levels 

if (L  ==  0) 

dEhfs  =  . 5*A85S12*K; 
elseif  (L  ==  1  &&  J  ==  .5) 
dEhfs  =  . 5*A85P12*K; 
elseif  (L  ==  1  &&  J  ==  1.5) 

dEhfs  =  ( . 5*A85P32*K)  +  B85P32* (  ( 1 . 5*K* (K+l ) 

2*Inam* (Inam+1) *J* (J+l) ) / (4*Inam* (2*Inam-l) *J* (2*J-1) ) ) ; 

else 

fprint f ( 1 \nIMPROPER  INPUT-Hyper f ine  energy  lvl  determinat ion\n\n ' ) 
dEhfs  =  0; 

end 

elseif  (Inam  ==  Inam87)  %  Rb87  hyperfine  energy  levels 

if (L  ==  0) 

dEhfs  =  . 5*A87S12*K; 
elseif (L  ==  1  &&  J  ==  .5) 
dEhfs  =  . 5*A87P12*K; 
elseif (L  ==  1  &&  J  ==  1.5) 

dEhfs  =  ( . 5*A87P32*K)  +  B87P32* ( ( 1 . 5*K* (K+l )  -... 

2*Inam* (Inam+1) *J* (J+l) ) / (4*Inam* (2*Inam-l) *J* (2*J-1) ) ) ; 

else 

fprint f ( ! \nIMPROPER  INPUT-Hyper f ine  energy  lvl  determinat ion\n\n ' ) 

dEhfs  =  0; 

end 

else 

fprint f (' \nIMPROPER  INPUT-Hyperf ine  energy  lvl  determinat ion\n\n ' ) 
dEhfs  =  0; 

end 

end  %end  of  function 

%Hyperfine  subfunction 
function  dviso  =  dvlso(lnam) 

%  Provides  Isotope  shift  for  Rb85  -  all  frequencies  relative  to  Rb87. 
global  vo85  vo87  Inam85; 

dvo  =  vo87  -  vo85;  %  Isotope's  hyperfine  spectrum  rel.  to  same 
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central  frquency  vo  so  that 
voRb85  +  dviso  =  voRb87 


if  (Inam  ==  Inam85) 
dviso  =  dvo; 

else 

dviso  =  0; 

end 

end 


%%  LINE  STRENGTH  FACTORS  - 

function  [S]  =  S_f un (F, Fprime, J, Jprime, I ) 
SI  =  Wigner6 j (J, Jprime, 1 , Fprime, F, I ) ; 

S  =  (2*Fprime+l) * (2*J+1) * (Sl~2)  ; 

end  %end  of  function 


%%  WIGNER  6 j  FACTORS - 

function  Wigner  =  Wigner6 j ( jl, j2, j3, Jl, J2, J3) 

%  Wigner  6j-symbol  calculator.  Written  by  Amita  B  Deb,  Clarendon  Lab.  2007. 
%  Improved  by  Richard  A.  Holt,  Univ.  of  Western  Ontario,  2009. 

%  Further  improved  to  deal  with  large  input  arguments  by  Lee  Harper, 

%  CRL  Oxford,  2009. 

%  Builds  on  the  m-file  Wigner 6 jcoeff  posted  by  Richard  A.  Holt. 

%  Calculates  {  jl,  j2  , j 3 }  using  Racah  formula. 

%  Jl  J2  J3 

%  Finding  Triangular  coefficients 

tril  =  triangle.coef f ( jl, j2, j3) ; 
tri2  =  triangle.coef f ( jl, J2, J3) ; 
tri3  =  triangle.coef f ( Jl, j2, J3) ; 
tri4  =  triangle.coef f ( Jl, J2,  j3) ; 


if  (tril  =  =  0  1 1  tri2  =  =  0  1 1  tri3  =  =  0  1 1  tri4==0) 
Wigner=0 ; 
return 


end 

%  Finding 

the  range  of 

a  ( 1 )  = 

jl 

+ 

j2 

+ 

j3; 

a  (2)  = 

jl 

+ 

J2 

+ 

J3 ; 

a  ( 3 )  = 

Jl 

+ 

j2 

+ 

J3 ; 

a  (4)  = 

Jl 

+ 

J2 

+ 

j3; 

rangei 

=  max (a) 

r 

summation 


in  the  Racah  formula. 
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Ml)  = 

jl 

+ 

j2 

+ 

Jl 

+ 

J2 ; 

M2)  = 

j2 

+ 

j  3 

+ 

J2 

+ 

J3 ; 

M3)  = 

j3 

+ 

jl 

+ 

J3 

+ 

J1; 

rangef 

= 

min (k) 

r 

Wigner 

= 

0; 

for  t=rangei : rangef 

Wigner  =  Wigner  +  ( (-1 ) ) *exp (gammaln (t+2 ) 
fung  (t,  jl,  j2,  j3,  Jl,  J2,  J3)  )  ; 


end 

Wigner  =  (tril*tri2*tri3*tri4 ) " (0 . 5) *Wigner; 


end 

%  -  Wigner6j  subfunctions  - 

function  r  =  fung(t,  jl,  j2,  j3, Jl, J2, J3) 

%  Calculating  the  logarithm  of  the  denominator  in  Racah  Formula,  using 
%  the  gamma  function  in  place  of  the  factorial. 

r  =  sum (gammaln ( [ (t  —  jl— j2  —  j3) ;  (t- j 1- J2- J3 ) ;  (t-Jl-j2-J3) ;  (t- Jl- J2- j  3 ) ;  .  .  . 

( jl+ j2+Jl+J2-t) ; ( j2+ j3+J2+J3-t) ; ( j3+ jl+J3+Jl-t) ]  +  1)); 


end 

o. _ 

o 

function  tri  =  triangle_coef f (a, b, c) 

%  Calculates  triangle  coefficients  for  angular  momenta. 

%  This  version  returns  0  if  the  triangle  inequalities  are  violated.  (RAH) 

if  (a<0  |  |  b<0  |  |  c<0) 

tri=0 ; 
return 

end 

for  xa  =  abs (a-b) : 1 : a+b 
if  c==xa 

tri  =  factorial (a+b-c) *factorial (a-b+c) *.. . 

factorial (-a+b+c) / (factorial (a+b+c+1)  )  ; 
return 

end 

end 
tri=0 ; 

end  %end  of  function 
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%%  CROSS  SECTION  AT  LINE  CENTER  - 

function  [sO]  =  sigma_0Rb (SFF, vD, x) 

%  Calculates  cross  section  at  line  center  for  each  hyperfine  transition 
global  c  vo85  vo87  tau  TransM  Inam85  Inam87; 

if  (TransM (x, 1 ) ==Inam85 ) 
vo  =  vo85; 

elseif  (TransM (x, 1 ) ==Inam87 ) 
vo  =  vo87; 

else 

fprint f ( ' \nERROR  in  sigma_0Rb  -  -  Inam  is  outside  allowed  range.\n\n' )  ; 
return; 

end 

dvnat  =  1/  (2*pi*tau) ; 

sO  =  SFF(x) * (1/tau) * ( ( (c*100) ~2) / (8*pi) ) * (1/ ( (vo+vD(x) ) ~2) ) *. . . 

( (2*TransM (x, 2) +1) / (2*. 5+1) ) * (2/ (pi * dvnat ) ) ; 


end 


%%  NUMBER  DENSITY  - 

function  [ndensity]  =  NumDensityRb (T) 

%  calculates  number  density  (crrU-3)  given  Temperature  in  Kelvin 
global  kB; 

ndensity  =  (RbVaporPressure (T) )/ (kB*T) *101325* (le-6)  ; 
end 


%%  VAPOR  PRESSURE  - 

function  [Pv]  =  RbVaporPressure (T) 

%  Provides  rubidium  vapor  pressure  in  atmospheres  given  T  in  Kelvin 

o. 

o 

%  REFERENCE: 

%  C.B.Alcock,  V.P.Itkin,  and  M.K.Horrigan 

%  Vapour  pressure  equations  for  the  metallic  elements:  298-2500K 
%  Canadian  Metallurgical  Quarterly,  1984,  23,  309-313 

o. 

o 

%  "Precise"  equation  for  Rb  vapour  pressure. 

%  Uses  melting  point  for  Rb  at  312. 45K  (39.30  degC) 

%  Valid  for  temperatures  from  melting  point  to  550K. 

o. 

o 

%  "Practical"  equations  for  RB  are  used  from  298K  to  melting  point. 

%  There  is  a  discontinuity  at  312. 45K.  Yields  +/-5%  or  better  precision. 

%  For  best  precision,  use  temperatures  between  melting  point  and  550K. 

if (T  >  312.45  &  &  T  <  550) 
a  =  8.316; 
b  =  -4275; 
c  =  -1.3102; 
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Pv  =  10."  (a  +  b*(T"-l)  +  c*loglO (T) ) ; 
elseif  (T  >  298  &  &  T  <  312.45) 
a  =  4.857; 
b  =  -4040; 

Pv  —  10. "(a  +  b  * ( T " — 1 ) ) ; 
elseif  (T  <  298  |  |  T  >  550) 

fprint f ( ' \nTemperature  is  outside  of  allowed  range. \n\n') 
return 

end 

end 


%%  CROSS  SECTION  - 

function  [cs]  =  CrossSectionRb (x, T, PHe, PCH4 ) 

%  Calculates  cross  section  as  function  of  frequency,  x 
global  tau  sO  vo85  vo87  TransM  vHF  numTrans  Inam85  Inam87; 


dvnat  =  1/  (2*pi*tau) ; 
csHF  =  zeros ( 1 , numTrans ) ; 


%  natural  linewidth,  a  nu  natural 
%  Cross  section  for  an  individual  HF  transition 


temp  =  0; 

for  n  =  1: numTrans 

Inam  =  TransM(n,l); 

Flwr  =  TransM(n,3); 

[gamma, shift , dvD, dvL]  =  BufferGasRb (T, Inam, PHe, PCH4 ) ; 

if  (Inam  ==  Inam85)  %  Set  central  frequency  based  on  isotope  (Inam) 

vo  -  vo85; 

elseif  (Inam  ==  Inam87) 
vo  -  vo87; 

else 

fprint f (’ \nERR0R  in  CrossSect ionRB-Inam  is  outside  range . \n\n ’) ; 
return; 

end 

csHF  =  RbNatAbundance ( Inam)  *  Bolt zRatioRb ( Inam, Flwr , T)  *  sO (n)  *  ... 

( (x  +  vHF (n)  +  shift  +  vo)  /  (vHF (n)  +  shift  +  vo) )  *  . . . 
gVoigt ( (x  -  vHF (n) ) , dvL, dvD) ; 
temp  =  temp  +  csHF; 

end 


cs  =  (dvnat*pi/2)  *  temp; 
end 


function  eta  =  RbNatAbundance ( Inam) 

%  Return  Rubidium  natural  abundance  based  on  Total  nuclear  angular  momentum 
Eta85  =  0.7217; 

Eta87  =  0.2783; 

if  (Inam  ==  5/2) 
eta  =  Eta85; 
return; 
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elseif  (Inam  ==  3/2) 
eta  =  Eta87; 
return; 

else 

fprint f ( ' \nlnam  is  outside  of  allowed  range  for  Rb85  or  Rb87.\n\n'); 
return; 

end 


end 


%%  BUFFER  GAS  EFFECTS  - 

function  [gamma, shift , dvD, dvL]  =  Buf ferGasRb (T, Inam, PHe, PCH4 ) 

%  Rb  buffer  gas  effects  -  pressure  broadening  and  frequency  shift 

o, 

o 

%  REFERENCE: 

%  Rotondaro,  M.  D.  &  Perram,  G.  P. 

%  Collisional  broadening  and  shift  of  the  rubidium  D1  and  D2  lines 
%  (5~2S_l/2  ->  5~2P_l/2,  5~2S_l/2  ->  5~2P_3/2)  by  rare  gases, 

%  H_2  ,  D_2  ,  N_2 ,  CH_4 ,  and  CF_4 

%  Journal  of  Quan.  Spect .  and  Radiative  Transfer,  1997,  57,  497-507 
global  tau  c  kB  MRb85  MRb87  vo85  vo87  Inam85  Inam87  TransM; 

Tdata  =  394.15;  %  Temperature  (K)  at  which  buffer  gas  values  were  measured 


%  Set  effect  constants 
if  (TransM (1,2) ==1/ 2 ) 
gammaHe  =  18.9e6; 
shiftHe  =  4.71e6; 
gammaCH4  =  29.1e6; 
shiftCH4  =  -7.92e6; 
else  %D2  values 

gammaHe  =  20.0e6; 
shiftHe  =  0.37e6; 
gammaCH4  =  26.2e6; 
shiftCH4  =  -7.00e6; 

end 


based  on  D1  or  D2  line 
%D1  values 
%  MHz/Torr 
%  MHz/Torr 
%  MHz/Torr 
%  MHz/Torr 

%  MHz/Torr 
%  MHz/Torr 
%  MHz/Torr 
%  MHz/Torr 


(J_upper)  . 


%  Set  central  frequency  and  mass  based  on  isotope 
if  (Inam  ==  Inam85) 
vo  =  vo85; 

MRb  =  MRb  8  5 ; 
elseif  (Inam  ==  Inam87) 
vo  =  vo87; 

MRb  =  MRb  8  7 ; 

else 

fprintf ( ' \nERR0R  in  Buf f erGasRb-Inam  is  outside  range . \n\n ') ; 
return ; 

end 

gamma  =  (gammaHe) * ( (T/Tdata) ~ (1/2) ) *PHe  + (gammaCH4) * ( (T/Tdata)  (1/2) ) *PCH4; 
shift  =  (shiftHe) * ( (T/Tdata)  ~  (1/2) ) *PHe  + (shiftCH4) *( (T/Tdata) ~ (1/2) )*PCH4; 


111 


408 

409 

410 

411 

412 

413 

414 

415 

416 

417 

418 

419 

420 

421 

422 

423 

424 

425 

426 

427 

428 

429 

430 

431 

432 

433 

434 

435 

436 

437 

438 

439 

440 

441 

442 

443 

444 

445 

446 

447 

448 

449 

450 

451 

452 

453 

454 

455 

456 

457 

458 

459 


dvD  =  (vo  +  shift) * ( (8*kB*T*log (2) ) / (MRb* (c~2) ) ) ~ (1/2) ; 

dvL  =  (1/ (2*pi*tau) )  +  gamma; 

end 


%%  BOLTZMAN  RATIO  - 

function  [BR]  =  Bolt zRatioRb ( Inam, F, T) 
global  kB  h  dv85S12  dv87S12  Inam85  Inam87; 

if  (Inam  ==  Inam85) 
if (F  ==  2) 

BR  =  5./ (5  +  7 *exp (- (h*dv85S12 ) . / (kB . *T) ) ) ; 
elseif  (F  ==  3) 

BR= ( 7*exp (- (h*dv85S12 )  . / (kB. *T) ) )  J  (5  +  7*exp (- (h*dv85S12)  . /  (kB. *T) ) ) ; 
elseif (F  ^  2  &&  F  ^  3) 

fprintf ( 1 \nIMPROPER  INPUTS\n\n’) 

BR  =  0; 

end 

elseif  (Inam  ==  Inam87) 
if (F  ==  1) 

BR  =  3 . / (3  +  5*exp (- (h*dv87S12) ./ (kB.*T) ) ) ; 
elseif (F  ==  2) 

BR= ( 5*exp (- (h*dv87S12 )  . / (kB. *T) ) )  J  (3+5*exp (- (h*dv87S12)  . /  (kB. *T) ) ) ; 
elseif  (F  ^  1  &&  F  #  2) 

fprintf (’ \nIMPROPER  INPUTS\n\n’) 

BR  =  0 ; 

end 

else 

fprintf (' \nIMPROPER  INPUTS\n\nf) 

BR  =  0; 

end 

end 


%%  VOIGT  LINESHAPE  - 

function  [y]  =  gVoigt (v, dvL, dvD) 

%  A  useful  approximation  for  the  Voigt  convolution  is  given  by  Thompson 

o, 

o 

%  REFERENCE: 

%  William  J.  Thompson 

%  Numerous  neat  algorithms  for  the  Voigt  profile  function. 

%  Computers  in  Physics,  Vol  7,  No.  6,  Nov/Dec  1993. 

Q. 

O 

%  gV(u,a)  =  Re [Exp (z~2) erfc (z) ] 

%  z  =  a  +  iu 

%  a  =  sqrt(ln(2))  vLorentzian/vDoppler 

%  u  =  2  sqrt(ln(2))  (v-vo) /vDoppler 

%  NOTE:  Not  valid  for  Doppler  free  spectra  where  vDoppler  =  0 
%  NOTE:  erfc(z)  =  1  -  erf(z)  (related  to  the  Faddeeva  function) 
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a  =  sqrt (log (2 ) ) * (dvL/dvD) ; 
u  =  2*sqrt (log (2 ) ) * (v/dvD) ; 
z  =  a  +  li*u; 

y  =  2*sqrt (log (2) /pi) * (1/dvD) *real (faddeeva  (li*z) ) ; 
end 


gVoigt  subfunctions 


function  w  =  f addeeva ( z , N) 

%  FADDEEVA  Faddeeva  function 

%  W  =  FADDEEVA (Z)  is  the  Faddeeva  function,  aka  the  plasma  dispersion 

%  function,  for  each  element  of  Z.  The  Faddeeva  function  is  defined  as: 

o, 

o 

%  w(z)  =  exp(-z~2)  *  erfc(-j*z) 

o, 

o 

%  where  erfc(x)  is  the  complex  complementary  error  function. 

o, 

o 

%  W  =  FADDEEVA (Z, N)  can  be  used  to  explicitly  specify  the  number  of  terms 
%  to  truncate  the  expansion  (see  (13)  in  [1]) .  N  =  16  is  used  as  default. 

o, 

o 

%  Example: 

%  x  =  linspace (-10, 10,  1001)  ;  [X,Y]  =  meshgrid (x, x)  ; 

%  W  =  faddeeva (complex (X, Y) ) ; 

%  figure; 

%  subplot (121) ;  imagesc (x, x,  real (W) )  ;  axis  xy  square;  caxis([-l  1]); 

%  title (' re (faddeeva ( z ) )  ' )  ;  xlabel ('re(z)  ' ) ;  y label ( ' im ( z )  ' ) ; 

%  subplot (122 ) ;  imagesc (x, x, imag (W) ) ;  axis  xy  square;  caxis([-l  1]); 

%  title (' im (faddeeva ( z ) )  ' ) ;  xlabel ('re(z)  ' )  ;  y label ( ' im ( z )  ' ) ; 

o. 

o 

%  Reference: 

%  [1]  J.A.C.  Weideman,  "Computation  of  the  Complex  Error  Function,"  SIAM 

%  J.  Numerical  Analysis,  pp.  1497-1518,  No.  5,  Vol .  31,  Oct.,  1994 

%  Available  Online:  http://www.jstor.org/stable/2158232 


if  nargin<2,  N  =  [];  end 
if  isempty(N),  N  =  16;  end 

w  =  zeros ( size ( z )) ;  %  initialize  output 


0,0,0,0  0, 
o  o  o  o  o 

%  for  purely  imaginary-valued  inputs,  use  erf  as  is  if  z  is  real 
idx  =  real(z)==0;  % 

w(idx)  =  exp (-z (idx) . ~2)  . *erfc (imag (z (idx) ) ) ; 

if  all (idx),  return;  end 
idx  =  -ddx; 


o,  o,  o,  o,  o, 
o  o  o  o  o 

%  for  complex-valued  inputs 

%  make  sure  all  points  are  in  the  upper  half-plane  (positive  imag.  values) 
idxl  =  idx  &  imag(z)<0; 


113 


512 

513 

514 

515 

516 

517 

518 

519 

520 

521 

522 

523 

524 

525 

526 

527 

528 

529 

530 

531 

532 

533 


z(idxl)  =  con j (z (idxl) ) ; 

M  =  2  *N; 

M2  =  2  *M; 

k  =  (-M+1:1:M-1)  f;  %  M2  =  no .  of  sampling  points. 

L  =  sqrt (N/sqrt (2 ) ) ;  %  Optimal  choice  of  L. 

theta  =  k*pi/M; 

t  =  L*tan (theta/2 ) ;  %  Variables  theta  and  t. 
f  -  exp (-t . "2)  . *  (L~2+t . ~2) ; 
f  =  [0;  f] ;  %  Function  to  be  transformed. 

a  =  real (f ft (f ftshift (f ) ) ) /M2;  %  Coefficients  of  transform, 
a  =  f lipud (a (2 : N+l ) ) ;  %  Reorder  coefficients. 

Z  =  (L+li*z (idx) )  . / (L-li*z (idx)  )  ; 
p  =  polyval (a, Z ) ;  %  Polynomial  evaluation. 

w(idx)  =  2 *p . / (L-li* z ( idx) ) . ~ 2  +  (1/sqrt (pi) ) . / (L-li*z (idx) ) ; 

%  convert  the  upper  half-plane  results  to  the  lower  half-plane 
w(idxl)  =  con j (2*exp (-z (idxl ) . ~2 )  -  w(idxl)); 

end 


Evaluate  w ( z ) . 
if  necesary 
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Absorption  Coefficient  Changes  for  Rb  D2 


The  previous  Matlab®  code  example  can  be  easily  modified  for  rubidium  D2  or 
any  other  atomic  species  with  multiple  isotopes.  To  change  the  Matlab®  absorption 
coefficient  function  on  page  C  for  rubidium  D2  make  the  following  modifications  at 
the  appropriate  places  within  the  code.  To  modify  use  of  this  code  for  other  atomic 
elements,  be  sure  to  modify  vapor  pressure  and  buffer  gas  effect  functions. 


%  Contstants . . 
%  Rb  D2  Data 
vo85  =  384.230 
vo87  =  384.230 
tau  =  26.2348e 
TransM  =  [5/2 
5/2 
5/2 
5/2 
5/2 
5/2 
3/2 
3/2 
3/2 
3/2 
3/2 
3/2 


4  0  637  3el2 ; 

4844685el2; 

-9; 

3/2  2  1; 

3/2  2  2; 

3/2  2  3; 

3/2  3  2; 

3/2  3  3; 

3/2  3  4; 

3/2  1  0; 

3/2  1  1; 

3/2  1  2; 

3/2  2  1; 

3/2  2  2; 

3/2  2  3] ; 


%  Frequency  (Hz)  Rb85 
%  Frequency  (Hz)  Rb87 
%  Lifetime  (s)  Rb85  and  Rb87 
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Sample  Cesium  Absorption  Coefficient  Code 


The  following  code  sample  illustrates  use  of  MATLAB®  functions  fnAbsCoef  f  CsD2  .m, 
fnAbsCoef  fCsDl  .m  and  fnAbsCoeff  CsD2SingleHFLine  .m  creating  figure  output  as 
shown  in  Figure  46. 


CsD,  Abs  Coeff  @  1 00°C 


Detuning  Frequency  (GHz) 
CsD2  Abs  Coeff  @100°C 


Detuning  Frequency  (GHz) 


Figure  46.  Sample  Absorption  Cesium  Coefficient  Code  figure  output 


%  Cesium  Absorption  Coefficient  D1  and  D2 

o, 

o _ 

%  Monte  D.  Anderson 

%  Department  of  Engineering  Physics,  Air  Force  Institute  of  Technology 
%  6  Aug  2010 

o, 

O - 

%  Calculates  cesium  D2  absorption  coefficent  (1/cm)  at  a  given  frequency 
%  at  temperature  T  (K)  and  buffer  gas  pressures  PN2,  PHe,  and  PCH4  (Torr) 

o. 

O - 

%  REQUIRES: 

%  f nAbsCoef f CsDl 

%  fnAbsCoef fCsD2 

%  fnAbsCoef fCsD2SingleHFLine 

%%  CLEAR 
close  all 
clear  all 
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clc 

tic 

%%  calcuations 

c  = 

2 . 997  92458e8 ; 

%m/  s 

voD2 

=  351725718; 

%MHz 

voDl 

=  335116049; 

%MHz 

nu  = 

linspace  (-8e9, 

8e9, 

T  = 

100  +  273.15; 

PN2 

=  0; 

PHe 

=  0; 

PCH4 

=  0; 

cell_length  =  5; 


absCoeffDl  =  zeros (length (nu) , 1 ) ; 
absCoeffD2  =  zeros (length (nu) , 1) ; 
lambdaD2  =  zeros (length (nu)  ,  1 )  ; 
lambdaDl  =  zeros  (length (nu) , 1 ) ; 

for  i=l : length (nu) 

absCoeffDl (i)  =  fnAbsCoef fCsDl (nu (i ) , T, PN2 , PHe, PCH4 ) ; 
absCoeffD2 (i)  =  fnAbsCoef fCsD2 (nu (i ) , T, PN2 ,  PHe,  PCH4 )  ; 
a32 (i)  =  fnAbsCoef f CsD2SingleHFLine (nu (i) , T, PN2 , PHe, 3, 2 ) ; 
a33 (i)  =  fnAbsCoef f CsD2SingleHFLine (nu (i) , T, PN2 , PHe, 3,  3 ) ; 
a34 (i)  =  fnAbsCoef f CsD2SingleHFLine (nu (i) , T, PN2 , PHe, 3, 4 ) ; 
a43 (i)  =  fnAbsCoef f CsD2SingleHFLine (nu ( i ), T, PN2 , PHe, 4 ,  3 ) ; 
a44 (i)  =  fnAbsCoef f CsD2SingleHFLine (nu (i) , T, PN2, PHe, 4, 4 ) ; 
a45 (i)  =  fnAbsCoef f CsD2SingleHFLine (nu ( i ), T, PN2 , PHe, 4 ,  5 ) ; 
lambdaD2 (i)  =  c. / (nu (i) * (le-9) +voD2*  (le-3) ) ; 
lambdaDl (i)  =  c ./ (nu (i) * (le-9) tvoDl*  (le-3) ) ; 

end 

%%  plot 
figure 

subplot (2,1,1); 

plot  (nu* (le-9) , absCoeffDl,  ' LineWidth ' , 2) ; 
title ([fCs  D_1  Abs  Coeff  @  ' , num2str (T-273 . 15) , ' "oC ' ] ) ; 
xlabel (' Detuning  Frequency  (GHz)1); 
ylabel('Abs  Coeff  (l/cm)’); 
axis  tight 
subplot (2, 1,2) 

plot  (nu* (le-9) , absCoef fD2,  ' LineWidth ' , 2)  ; 
hold  on 

plot  (nu* (le-9) ,a32,  1 LineWidth ' , 1) ; 
plot  (nu* (le-9) ,a33,  1 LineWidth ' , 1) ; 
plot  (nu* (le-9) ,a34,  1 LineWidth' , 1) ; 
plot  (nu* (le-9) ,a43,  ' LineWidth' , 1) ; 
plot  (nu* (le-9) ,a44,  ' LineWidth' , 1) ; 
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plot (nu* (le-9) , a45, 1 LineWidth ' , 1) ; 

title(['Cs  D_2  Abs  Coeff  @  * , num2str (T-273 . 15) ,  ' ~oC 1 ] )  ; 

xlabel (' Detuning  Frequency  (GHz)'); 

ylabel('Abs  Coeff  (1/cm)'); 

axis  tight 

hold  off 


toe 

%  tblSaveDl  =  [nu'*(le-9)  lambdaDl  absCoeffDl] ; 

%  tblSaveD2  =  [nu'*(le-9)  lambdaD2  absCoeffD2]; 

%  save  'C:\SimAbsCoeff-Dl.dat'  tblSaveDl  -ascii  -double; 
%  save  'C:\SimAbsCoeff-D2.dat'  tblSaveD2  -ascii  -double; 
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Sample  Rubidium  Absorption  Coefficient  Code 


The  following  code  sample  illustrates  use  of  MATLAB®  function  fnAbsCoef  f  RbDl  .m, 
creating  figure  output  as  shown  in  Figure  47. 


Figure  47.  Sample  Rubidium  Absorption  Coefficient  Code  figure  output 


%  RbDIAbsCoef fMaker  .m 

o. 

O - 

%  Monte  D.  Anderson 

%  Department  of  Engineering  Physics,  Air  Force  Institute  of  Technology 
%  26  Feb  2010 

o, 

o - 

%  Calculates  absolute  aborption  coefficient  from  temperature  in  Kelvin 
%  and  buffer  gas  pressures 

o, 

o _ _ _ _ 

clear  all;  close  all;  clc 
tic 

startGHz  =  -5.5; 
stopGHz  =  7; 
numPts  =  501; 

aNu  =  (stopGHz  -  startGHz)/ (numPts-1); 

T  =  273.15  +  110;  %temp  in  Kelvin 
PHe  =  0;  %Torr 
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PCH4  =  0;  %Torr 

c  =  2.997  92458e8 ;  %m/s 

vo87  =  377 . 1074633805e6;  %  Frequency 

(MHz) 

Rb87 

alpha  =  zeros  ( 1 , numPts ) ; 
nuref  =  zeros  ( 1 , numPts ) ; 
lambda  =  zeros ( 1 , numPts ) ; 

numberl  =  0; 

for  freq  =  startGHz : aNu : stopGHz  %GHz 

from 

nu  ref 

numberl  =  numberl+1; 

alpha (numberl )  =  fnAbsCoef fRbDl (f req*le9, T, PHe, PCH4 ) ; 
nuref (numberl )  =  freq; 

lambda (numberl )  =  c . / (freq+vo87* (le-3) ) ; 

end 

figure 

plot (nuref, alpha, ' b ' , ' LineWidth ' , 2 ) ; 

xlabel (’ Detuning  (GHz)');  ylabel (' Absorption  Coefficient  (1/cm)'); 
title (['Rb  D_1  @  ' , num2str (T-273 . 15) , ' "oC ’ ] ) ; 
grid  on 

figure 

plot ( lambda, alpha, '  b  '  ,  ' LineWidth '  ,  2 )  ; 

xlabel (' Wavelength  (nm) ' ) ;  ylabel (' Absorption  Coefficient  (1/cm)’); 
title (['Rb  D_1  @  ' , num2str (T-273 . 15) ,' "oC '])  ; 
grid  on 

%  I/Io  - 

%  Calculates  intensity  ratio  I/Io 

%  cellLength  =  5;  %  Cell  length  in  cm 

o, 

o 

%  Iovrlo  =  zeros  (1, 2001) ; 

%  numberl  =  0; 

%  for  freq  =  -10:. 01:10 
%  numberl  =  numberl+1; 

%  Iovrlo (numberl )  =  exp (-fnAbsCoef fRbDl ( freq*le9,  T,  PHe,  PCH4 ) *cellLength)  ; 
%  end 

o. 

o 

%  freq  =  -10 :  . 01  : 10; 

%  figure 
%  hold  on 

%  plot ( freq, Iovrlo, ' r ' , ' LineWidth ' , 2 ) ; 

%  xlabel (' Detuning  (GHz)');  ylabel ('I/Io  (unitless)'); 

%  axis (  [-10  10  0  1] )  ; 

%  grid  on 
%  hold  off 

%  Save  Data  - 

tblSave  =  [lambda'  alpha']; 
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72  %save 


' C : \SimRbDlAbsCoef f-lOOC . txt 


tblSave  -ascii  -double; 
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Kramers-Kronig  Group  Delay  Prediction 


The  following  code  calculates  delay  from  absorption  coefficient  using  the  Kramers- 
Kronig  model.  This  code  can  be  easy  modified  for  any  simulated  absorption  coefficient 
or  used  with  observed  absorption.  Sample  output  is  shown  in  Figure  48  and  Figure  49. 


Cs  D2  Group  Velocity  @  1 00°C 


Figure  48.  Group  velocity  output  from  Kramers-Kronig  delay  calcula¬ 
tion. 


%  DelayPrediction .m 

o, 

o _ _ 

%  Monte  D.  Anderson 

%  Department  of  Engineering  Physics,  Air  Force  Institute  of  Technology 
%  19  Mar  2010 

o. 

o _ _ 

%  Calculates  predicted  delay  given  absorption  coefficient  using  the 
%  Kramers-Kronig  delay  model  outlined  in 

o. 

o 

%  Monte  D.  Anderson,  Glen  P.  Perram 

%  Slow  light  in  cesium  vapor:  pulse  delay  measurements  and  predicted  delay 
%  (Proceedings  Paper)  Published  11  February  2010 
%  Vol .  7612:  Advances  in  Slow  and  Fast  Light  III, 

%  Selim  M.  Shahriar;  Philip  R.  Hemmer,  Editors,  76120A 

o, 

o - 

%  REQUIRES: 
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Cs  D2  Delay  @  100°C 


Figure  49.  Delay  output  from  Kramers-Kronig  calculation. 


%  fnAbsCoef fCsD2  .m 

o, 

o _ 

clear  all; 
close  all; 
clc; 
tic 


%  Constants  - 

%  Physical  Constants  (NIST  CODATA  2006  values) 
c  =  2 . 99792458e8;  %  vacuum  speed  of  light  (m/s) 


%  Cs  D2  Data 

vo  =  351 . 72571850el2;  %  Frequency  (Hz) 


Scenario  Information 


T  =  100; 

PN2  =  0; 

PHe  =  0; 

PCH4  =  0; 
vmin  =  -10100E6; 
vmax  =  10100E6; 
NumPts  =  10000; 


%Temporary  temperature  (C) 

%N2  pressure  (Torr) 

%He  pressure  (Torr) 

%CH4  pressure  (Torr) 

%ref  frequency  in  Hz  from  center  frequency 
%ref  frequency  in  Hz  from  center  frequency 
%Total  number  of  points  will  be  NumPts+1 


A_nu  =  (vmax  -  vmin) /NumPts; 


%discrete  step  size  in  Hz 


%  Absorption  Coefficient  - 

tblSimalpha  =  zeros (NumPts+1 , 2 ) ; 
for  n  =  l:NumPts+l 
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tblSimalpha (n,  1 )  =  vmin  -  A_nu  +  n*A_nu; 

tblSimalpha (n, 2 )  =  fnAbsCoef fCsD2 (tblSimalpha (n, 1 ), 273 . 15+T, .. . 

PN2 , PHe, PCH4 ) ; 

end 

%  figure 

%  plot (tblSimalpha ( : , 1 ) , tblSimalpha ( : , 2 ) ) 

%  Real  Index  of  Refraction  - 

tl  =  zeros ( 1 , NumPts+1 ) ; 
tblnRe  =  zeros (NumPts+1 , 2 ) ; 
for  p  =  OiNumPts 

for  m  =  OiNumPts 
if (m  #  p) 

tl(m+l)  =  (tblSimalpha (m+1 ,  2 )/( le-2 ))/.. . 

( (m*A_nu  +  vmin  +  vo) ~2  -  (p*A_nu  +  vmin  +  vo) ~2) ; 

else 

t 1 (m+1 )  =  0 ; 

end 

end 

tblnRe (p+1 , 1 )  =  p*A_nu+vmin; 

tblnRe  (p+1,2)  =  1+  (  (  (2 . 9  9792458e8)  /  (2*  (pi~2) ) )  *sum(tl)  *A_nu)  ; 

end 

%  figure 

%  plot (tblnRe (i, 1 ), tblnRe , 2 ) ) 

%  dn/dv  - 

tbldndv  =  zeros (NumPts+1 , 2 ) ; 
for  n  =  2:NumPts 

tbldndv (n, 2)  =  (tblnRe (n+1 , 2 )  -  tblnRe (n-1, 2) )/ (2*A_nu) ; 

end 

tbldndv(:fl)  =  tblnRe (:, 1) ; 
tbldndv (1,2)  =  0; 
tbldndv (NumPts+1 , 2 )  =  0; 

%  figure 

%  plot (tbldndv ( : , 1 ) , tbldndv ( : ,  2 ) ) 

%  Group  Velocity  - 

tblVg  =  zeros (length ( 0 : NumPts ) , 2 ) ; 
for  n  =  OrNumPts 

tblVg (n+1,1)  =  n*A_nu  +  vmin; 

tblVg (n+1, 2)  =  1/ (tblnRe (n+1 , 2 )  +  (tblnRe (n+1 , 1 )  +  vo) *tbldndv (n+1 , 2 ) ) ; 

end 

%  figure 

%  plot (tblVg (:,l)*(le-9), tblVg ( : , 2 ) ) 

%  Delay  - 

cell_length  =  5;  %cell  length  in  centimeters 
tblDelay  =  zeros (length ( 0 : NumPts ), 2 ) ; 
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for  n  =  0:NumPts 

tblDelay (n+1 , 2 )  =  cell-length/ (tblVg(n+l,  2) *2 . 99792458el0)  ; 
tblDelay (n+1 ,  1 )  =  n*A_nu+vmin; 

end 

figure 

plot (tblDelay  (:, 1 ) , tblDelay (:, 2 ) ) 

%  Plot/Export  Data  Files  - 

%save  'C:\SimAlpha.dat'  tblSimalpha  -ascii  -double; 

tblSaveVg  =  [tblVg ( :  ,  1 ) * (le-9)  tblVg(:,2)]; 
figure 

plot (tblSaveVg ( : , 1 ) , tblSaveVg (:,2), 'LineWidth',2) ; 
title (['Cs  D_2  Group  Velocity  @  ' , num2str (T) , ' ~ OC ' ] ) ; 
xlabel (' Detuning  Frequency  (GHz)'); 
ylabel (' Group  Velocity  (c) ' ) ; 
axis ( [-10  10  -0.01  0.2] ) ; 

%  save  'C:\SimVg.dat'  tblSaveVgrp  -ascii  -double; 

tblSaveDelay  =  [tblDelay (:, 1) * (le-9)  tblDelay (:, 2 )*( le9 )] ; 
figure 

plot  (tblSaveDelay  (■.,!),  tblSaveDelay  (  :  ,  2  )  ,  '  LineWidth  '  ,  2  )  ; 

title (['Cs  D_2  Delay  @  ' , num2str (T) ,  ' ~oC ' ] )  ; 

xlabel (' Detuning  Frequency  (GHz)'); 

ylabel ( 'Delay  (ns)  ' )  ; 

axis ( [-10  10  -0.5  25] )  ; 

%  save  'C:\SimDelay.dat'  tblSaveDelay  -ascii  -double; 
toe 
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Hole  Burning  Locations 


The  following  code  calculates  hole  burning  locations  for  a  cesium  D2  pump  and 
corresponding  hole  locations  in  cesium  D i,  producing  the  following  output  shown  in 
Figure  50. 


Figure  50.  Hole  location  calculation  output. 


1  %  HoleBurningLocations . m 

2  % _ 

3  %  Monte  D.  Anderson 

4  %  Department  of  Engineering  Physics,  Air  Force  Institute  of  Technology 

5  %  19  Apr  2010 

6  % 

7  %  Given  a  frequency  relative  to  the  Cs  D2  absorption,  calculates  the 

8  %  difference  from  the  D2  hyperfine  transition  frequencies,  and  plots 

9  %  the  same  difference  relative  to  all  the  D1  transitions.  This  should 
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give  relative  positions  of  the  hole  burnt  in  the  D1  absorption  due  to 
the  hole  burnt  in  the  D2  spectrum. 


o, 
o 

o, 
o 

o, 
o 

o, 
o 

%  REQUIRES 

%  fnAbsCoef fCsD2 . m 

%  fnAbsCoef fCsDl  .m 

%  hfsEnergy.m 

o, 

o - 

clear  all; 
close  all; 
clc; 

%  CONSTANTS  :::::::::::::::::::::::::::::::::::::::::::::: 

global  c  eo  h  hbar  kB  AMU; 
global  AmCs  MCs  Inam  voDl  voD2  tau; 

%  Physical  Constants  (NIST  2006  CODATA  recommended  values) 
c  =  2.997  92458e8 ; 
eo  =  8 . 85418781762e-12; 
h  =  6 . 62  60  68  96e-34 ; 
hbar  =  1 . 054571628e-34; 
kB  =  1 ,3806504e-23; 

AMU  =  1 . 660538782e-27; 

%  Cesium  Atomic  Data 
AmCs  =  132.905451931; 

MCs  =  AmCs * AMU; 

Inam  =  7/2; 

%  Cs  D2/D1  Data 
voD2  =  351 . 72571850el2 ; 
voDl  =  335. 11604 9el2; 
tau  =  30.32e-9; 

%Transition  Definitions 
TransMD2  =  [7/2  3/2  3  2; 

7/2  3/2  3  3; 

7/2  3/2  3  4; 

7/2  3/2  4  3; 

7/2  3/2  4  4; 

7/2  3/2  4  5] ; 
numTransD2  =  length (TransMD2) ; 

TransMDl  =  [7/2  1/2  3  3; 

7/2  1/2  3  4; 

7/2  1/2  4  3; 

7/2  1/2  4  4]; 

numTransDl  =  length (TransMDl) ; 

vDl  =  zeros ( 1 , numTransDl ) ; 
for  j  =  1: numTransDl 

vDl ( j )  =  (hf sEnergy (TransMDl ( j , 1 ) , 1 , TransMDl ( j , 2 ) , TransMDl (j,4) ) /h) - . . . 
(hf sEnergy (TransMDl ( j , 1 ) , 0,1/2, TransMDl ( j  ,  3 ) ) /h)  ; 

end 


%  vacuum  speed  of  light  (m/s) 

%  vacuum  permittivity  (F/m) 

%  Planck  const  (Js) 

%  Planck  const/2pi  (Js) 

%  Boltzmann  const  (J/K) 

%  Atomic  Mass  Unit  (kg) 

%  Atomic  mass  (AMU) 

%  Mass  of  Cs  atom  (kg) 

%  Total  nuclear  angular  momentum 

%  Frequency  (Hz) 

%  Frequency  (Hz) 

%  Lifetime  (s) 
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dvS12  =  hf sEnergy ( Inam, 0 , 1/2 ,  4 ) /h  -  hf sEnergy ( Inam, 0 , 1 /2 , 3 ) /h; 

vD2  =  zeros ( 1 , numTransD2 ) ; 
for  j  =  l:numTransD2 

vD2 ( j )  =  (hf sEnergy (TransMD2 ( j , 1 ) , 1 , TransMD2 ( j , 2 ) , TransMD2 (  j  ,  4  ) )/h)-. . . 
(hf sEnergy (TransMD2 (j,l), 0,1/2, TransMD2 ( j ,  3) ) /h) ; 

end 


%  MAIN  . 

nu  =  linspace (-8e9, 8e9,  1000 )  ; 

T  -  50; 

PN2  =  0; 

PHe  -  0; 

PCH4  =  0; 

CellLength  =  5; 


nuProbe  = 

-  3 . 8  4  e  9  ; 

%  11732.1790 

o, 

o 

nuProbe 

=  -3 . 86e9; 

%  11732.1783 

o, 

o 

nuProbe 

=  -3 . 924e9; 

%  11732.1762 

o, 

o 

nuProbe 

=  -  3 . 9  0  e  9 ; 

%  11732.1770 

dnuProbe  =  nuProbe  -  vD2; 

nuabs  =  voD2  +  nuProbe; 
lambda  =  c/nuabs; 
wavenum  =  nuabs/c; 

fprintf('D2  Frequency  - \nf); 

fprint f ( '  %3 . 4 f  nm\n ' , lambda* ( le9 ) )  ; 
fprint f ( '  %5 . 4 f  l/cm\n ' , wavenum* ( le-2 ) ) ; 
strwn  =  sprint f (' %5 . 4 f '  , wavenum* ( le-2 )) ; 
fprint  f  (  ' - \n  '  )  ; 

holesDl  =  zeros (length (vDl ) , length (dnuProbe) ) ; 
holeColorStr  =  zeros (length (vDl ) , length (dnuProbe) ) ; 
for  i=l : length (vDl ) 

for  j  =  1 : length (dnuProbe) 

holesDl(i,j)  =  vDl (i)  -  dnuProbe (j); 

end 

end 

absCoeffDl  =  zeros ( length (nu) , length (T) ) ; 
absCoeffD2  =  zeros (length (nu) , length (T) ) ; 
lambdaDl  =  zeros (length (nu)  ,  length (T) )  ; 

IRatio  =  zeros  (length (nu) , length (T) ) ; 
for  i=l : length (nu) 

absCoeffDl (i, j )  =  fnAbsCoeffCsDl(nu(i),T+273.15fPN2,PHefPCH4); 
absCoef fD2 (i,  j )  =  fnAbsCoef fCsD2 (nu (i) , T+273 . 15f PN2,  PHef PCH4 )  ; 
lambdaDl (if j )  =  c. / (nu (i) * (le-9) +voDl* (le-3) ) ; 

IRatio (i,j)  =  exp (-absCoef fDl (if j ) *CellLength)  ; 

end 
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%  PLOTS  . 

peak  =  max (max (max (absCoef fDl ) ) , max (max (absCoef fD2 ) ) ) ; 
ylgt  =  linspace (0, peak, 250) ; 
ydrk  =  linspace (0, peak,  1000 )  ; 

plotmin  =  -6; 
plotmax  =  -2.2; 

figure 

subplot  (2 , 1 , 1 ) 
hold  all 

plot (nu*  (le-9) , exp (-absCoef fD2* Cel lLength) , 'b ' ) ; 
ax  =  axis; 

axis ( [plotmin  plotmax  0  1]); 

plot ( [nuProbe* ( le-9)  nuProbe*  (le-9) ] ,  [0  1 ] ,  ' c ' ,  1 LineWidth f , 3) ; 
h  =  text (nuProbe* (le-9), ax (4) /2, [  '  Pump  @  ' , strwn, 1  cm"{-l}'],... 
' Vertical Alignment ' , ' middle ' , . . . 

' HorizontalAlignment ' ,  f  center ' ,  .  .  . 

' FontSize ' , 9) ; 
set  (h,  'rotation',  90) 
for  i  =  l:numTransD2 

plot ( [vD2 (i) * (le-9)  vD2 (i) * (le-9) ] , [0  l],'b:') 

end 

xlabel (' Detuning  Frequency  \ nu_{\Delta2 }  (GHz)'); 

ylabel (' Transmittance  (unitless) '); 
hold  off 

subplot (2,1,2) 
hold  all 

plot ( [vDl (1) * (le-9)  vDl (1) * (le-9) ]  ,  [0  l],'r') 
ax  =  axis; 

axis ( [plotmin  plotmax  0  1]); 

plot  (nu*  (le-9) , exp (-absCoef fDl*CellLength)  ,  ' r '  )  ; 
for  i  =  l:numTransDl 

plot ( [vDl (i) * (le-9)  vDl (i) * (le-9) ] , [0  l],'r:') 

end 

for  i=l : length (vDl ) 

for  j  =  1 : length (dnuProbe) 

plot ( [holesDl (i, j) * (le-9)  holesDl (i, j) * (le-9) ] ,  [0  1]  ,  .  .  . 

' c ' , ' LineWidth ' , 2 ) ; 

end 

end 

xlabel (' Detuning  Frequency  \ nu_{\Deltal } ( GHz  )  '  )  ; 
ylabel  (' Transmittance  (unitless)  '); 
hold  off 
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Appendix  D.  Delay- Absorbance  Plots 


Introduction 

Most  of  the  delay  plots  in  this  manuscript  are  given  in  terms  of  delay  as  a  function 
of  detuning  frequency.  Other  graphical  interpretations  are  useful  if  there  is  physically 
meaningful  information  which  can  be  extracted.  In  EIT  slow  light  schemes,  Kasapi  et 
al.( Kasapi  et  al.,  1995)  used  a  particularly  insightful  plot  with  absorbance  presented 
as  a  function  of  delay.  The  slope  of  the  resulting  EIT  delay  was  given  as  the  linewidth 
of  the  transitions  (Kasapi  et  al.,  1995). 

Similar  graphical  representations  of  the  absorbance  as  a  function  of  delay  can 
be  created  for  delays  observed  in  the  linear  dispersion  regime.  The  following  plots 
present  theoretical  absorbance-delay  plots  from  the  linear  dispersion  Kramers-Kronig 
delay  predictions  and  delay  data  from  the  preceding  chapters. 

Temperature  Dependence 

Temperature  dependence  of  the  delay-absorbance  curves  are  displayed  in  Figure  51 
and  Figure  52.  The  initial  slope  of  the  curve  appears  to  follow  a  common  asymptote, 
but  the  curves  depart  at  longer  delays  at  higher  temperatures. 

Pressure  Dependence 

Pressure  dependent  delay-absorbance  curves  are  shown  in  Figure  53  and  Figure  54 
for  several  simulated  helium-broadened  samples  from  0  Torr  to  20  Torr.  Notice  how 
the  curves  no  longer  follow  a  common  asymptotic  limit  compared  to  the  temperature- 
dependent  delay-absorbance  curves  thus  the  initial  slope  is  clearly  a  function  of  pres¬ 
sure. 
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Figure  51.  Delay- Absorbance  plot  showing  temperature  dependence  be¬ 
tween  40  °C  and  90  °C. 


Figure  52.  Delay- Absorbance  plot  showing  temperature  dependence  be¬ 
tween  40  °C  and  90  °C  showing  detail  of  delay  between  0  ns  and  40  ns. 
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Figure  53.  Delay- Absorbance  plot  showing  pressure  dependence  between 

0  Torr  and  20Torr. 


Figure  54.  Delay- Absorbance  plot  showing  pressure  dependence  between 
OTorr  and  20  Torrshowing  detail  of  delay  between  0  ns  and  75  ns. 


132 


Linear  Dispersion  Delay 


Linear  dispersion  data  from  Chapters  III  and  IV  are  used  to  plot  absorbance  in 
Figure  55  as  a  function  of  delay  and  the  predicted  delay-absorbance  curves  (from 
simulation)  are  provided  as  reference. 

Hole-burning  Delay 

Delay-absorbance  curves  in  the  regions  surrounding  the  rapidly  varying  absorption 
created  by  spectral  hole  burning  are  shown  in  Figure  56  and  Figure  57.  Data  showing 
all  three  pump  powers  (0  rnW,  5  rnW  and  16  mW)  are  shown  in  Figure  56.  No¬ 
pump  (0  mW)  delay-absorbance  curve  is  shown  in  Figure  57  with  the  predicted  curve 
for  reference.  Deviation  from  the  predicted  curve  may  be  indicative  of  hole  burning 
effects. 


Figure  55.  Delay- Absorbance  plot  showing  linear  dispersion  delay  as 
shown  in  Figure  7  with  predicted  delay-absorbance  curves. 
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Figure  56.  Delay- Absorbance  plot  showing  hole-burning  data  as  shown 
in  Figure  24. 


Figure  57.  Delay- Absorbance  plot  showing  0  mW  (no  pump)  hole¬ 
burning  data  as  shown  in  Figure  24  with  predicted  delay-absorbance 
curves. 
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